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ABSTRACT 


A  method  o£  fairing  ships'  hulls  and  obtaining  a  surface 
equation  representing  the  ships '  hull  through  the  use  of 
linear  programming  is  presented.  Methods  for  producing 
loft  data  from  this  equation  are  also  presented. 
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Section  I 


INTRODUCTION 


The  mathematical  tool  called  "Linear  Programming"  has  been  in  existence 
for  many  years.  Since  1951,  when  George  Dantzig  first  published  the 
simplex  method  for  solving  linear  programming  problems,  and  through 
increased  availability  of  large-scale  digital  computers,  linear  program¬ 
ming  use  has  become  widespread.  Primarily,  this  use  of  linear  programming 
has  been  confined  to  economic  studies  of  product  mixes,  ship  scheduling, 
plant  locations,  etc.  In  this  report  a  new  use  of  linear  programming  is 
described:  that  of  curve  and  surface  fitting. 

This  use  of  linear  programming  has  enabled  the  formulation  of  a  more 
sophisticated  approach  to  the  problem  of  defining  a  ship's  hull  mathe¬ 
matically.  Specifically,  it  is  possible  to  require  that  the  curve  or 
surface  meet  certain  conditions  of  fairness  and  that  the  curve  come  as 
close  to  the  offsets  as  possible  while  meeting  these  conditions. 

Also  presented  are  some  methods  for  extracting  the  information  required 
to  construct  the  ship  from  the  faired  curve  or  surface. 


T.  C.  Koopmans,  ed.,  Activity  Analysis  of  Production  and  Allocation. 
Chap.  XXI,  New  York,  Wiley,  1951 


1.0. 0-1 


1-1 


Section  II 


CRITERIA  AMD  BACKGROUND 


CRITERIA 

The  title,  "Mathematical  Lofting,"  is  meant  to  apply  to: 

1.  The  overall  process  of  mathematically  fairing  a  ship's  hull 

2.  Providing  from  this  hull,  all  of  the  information  necessary 
to  cut  and  shape  parts  influenced  by  the  hull 

Mathematical  fairing  is  the  basic  key  to  this  process;  however,  mathe¬ 
matical  fairing  becomes  most  efficient  when  integrated  into  a  system 
which  meets  the  following  criteria: 

(1)  The  system  should  accept  as  data  the  preliminary  information 
on  hull  definition  developed  by  the  naval  architect. 

(2)  It  should  include  a  method  for  eliminating  gross  errors  in 
this  input  data. 

(3)  It  should,  from  this  data,  obtain  a  fully  faired  hull  form 
from  which  the  ship  can  be  produced. 

(4)  The  fairing  process  should  operate  with  as  little  human 
intervention  as  possible. 

(5)  The  system  should  develop,  as  much  as  is  practicable,  complete 
loft  information  now  provided  by  the  manual  loft. 

Criterion  (3)  requires  that  the  hull  be  fully  faired.  The  term  "fair" 
has  been  in  use  for  a  long  time  with  many  varied  interpretations .  To 
provide  a  basis  for  a  mathematical  fairing  method,  we  must  reduce  the 
subjective  definition  of  the  word  "fair"  to  a  more  precise  mathematical 
definition.  Curves  which  meet  the  following  mathematical  requirements 
also  meet  the  subjective  requirements  of  experienced  loftsmen: 
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1-a  The  curve  Y(X)  must  be  continuous 

1-b  Its  first  derivative  must  be  continuous 

l-c  Its  second  derivative  muse  be  continuous 

1-d  The  curve  must  be  completely  free  of  extraneous  inflection 
points  while  possessing  those  inherent  in  the  data. 

1-e  Deviation  from  scaled  offsets  must  be  as  small  as  possible, 
proved?*}  above  four  conditions  hold. 

1-f  Curve  must  look  good  to  the  eye. 

Condition  1-e  is  added  for  completeness  when  beginning  with  scaled  offsets; 
it,  and  1-f  are  also  the  only  conditions  possibly  possessing  any  subjecti¬ 
vity  in  their  interpretation.  Condition  1-d  is  imperative  if  the  faired 
curves  are  to  be  used  for  construction  purposes  as  opposed  to  mathematical 
expressions  vjsed  for  hydrodynamic  studies  only. 

A  system  designed  to  meet  the  above  criteria  would  start  with  a  table 
of  roughly  faired  offsets  and  other  information  from  the  naval  architect's 
preliminary  lines  drawing.  It  should  then  provide  a  method  for  eliminating 
any  gross  errors  due.  for  example,  to  misreading  a  scale  or  transcribing  the 
wrong  information. 

From  this  data,  it  should  provide,  with  a  minimum  of  human  intervention 
or  decision,  a  mathematical  description  of  a  fair  hull.  Next,  it  must 
be  able  to  extract  from  the  mathematical  definition  of  the  hull  sufficient 
information  to  produce  the  ship.  A  diagram  of  such  a  system  is  presented 
in  Fig.  1. 


Continuity  of  first  and  second  derivatives  is  equivalent  to  continuous 
slope  and  curvature. 


1.0. 0-1 


II-2 


fMWAL  ARCHITECT* 
PRIUHlNARr  LINE  I 

A 

OTHER  BATA 


6 AT A  PREPARATION 
TOR  OTHER  LOFTINC 
INFORMATION 


1.0. 0-1 


II-3 


BACKGROUND 


An  extensive  study  of  most  methods  proposed  for  mathematical  lofting  was 
undertaken  prior  to  the  initial  formulation  of  the  present  system.  This 
study  revealed  that  all  proposed  systems  fall,  in  one  way  or  another, 
tc  meet  the  criteria  required  of  a  system  oriented  toward  production. 
From  the  system  which  most  fully  met  the  criteria  (method  proposed  by 
F.  Thellheimer*),  was  borrowed  the  class  of  function  with  which  to 
represent  ships-form,  as  well  as  the  notation  for  this  function. 


« 


» 


e 


•fg 

Feodor  Thellheimer  and  William  Starkweather,  The  Fairing  of  Ship  Lines 
on  a  High-Speed  Computer.  David  W.  Taylor  Model  Basin,  Appl.  Math.  Lab., 
Report  No.  1474,  Jan  1961 
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Section  III 


MATHEMATICAL  FAIRING 


GENERAL  DISCUSSION 

The  manual  loft  has  used  the  batten  (graphical-spline)  very  satisfactorily 
for  interpolating  between  points  as  well  as  "fairing"  points.  We  feel 
it  would  be  desirable  to  mathematically-fair  with  the  mathematical 
equivalent  of  the  graphical-spline,  the  "analytic-spline." 

Analytic  Spline 

Thin  beam  theory  (based  on  the  Bernoulli-Euler  Equation)  shows  that  a 
thin  beam  (e.g.,  batten)  supported  at  a  finite  number  of  points,  under 
small  deflections,  can  be  represented  by  a  series  of  cubic  equations 
(a  single* cubic  in  each  interval  defined  by  two  successive  points  of 
support)  with  any  two  neighboring  cubics  having  equal  ordinate  values, 
first  derivatives  and  second  derivatives  at  the  point  of  support  common 
to  both  curves.  The  third  derivatives  at  this  common  point  are  not 
necessarily  equal,  resulting  in  a  discontinuous  third  derivative  for 
the  curve  representing  the  spline.  The  discontinuity  of  the  third 
derivative  at  the  Juncture  point  is  proportional  to  the  force  applied 
to  the  beam  at  this  point. 

This  situation  closely  parallels  the  one  existing  when  ships'  lines  are 
drawn  from  battens,  and  again  when  these  curves  are  faired  full-scale 
on  the  mold  loft  floor. 

Definition:  An  analytic-spline  through  a  set  of  points  is 
a  curve  through  these  same  points  consisting  of  a  cubic  in 
each  interval  defined  by  two  successive  points  with  any 
two  neighboring  cubics  having  equal  ordinates,  firs;: 
derivatives,  and  second  derivatives  at  the  point  common  to 
both  curves. 


Returning  to  the  original  requirements  for  fairness  of  a  line  (Page  II-l) 
it  can  clearly  be  seen  that  all  an  analytic-spline  needs  in  order  for  it 
to  satisfy  all  fairness  criteria,  is  that  it  possess  only  those  inflec¬ 
tion  points  inherent  in  the  data,  while  minimizing  the  deviation  between 
curve  and  points. 

Having  settled  on  the  type  of  curve  with  which  to  represent  ships-form, 
we  must  now  devise  a  procedure  for  minimizing  the  deviation  between 
curve  and  points  while  requiring  the  occurrence  of  inflection  points  to 
be  compatible  with  the  Information  inherent  in  the  data. 

One  point  we  would  like  to  emphasize  is  that  in  curve -fairing,  as  opposed 
to  curve-fitting,  the  main  objective  is  not  to  pass  exactly  through  the 
data  points,  but  that  the  curve  be  pleasing  to  the  eye  and  possess  only 
those  inflection  points  as  indicated  by  the  data.  Only  after  these 
fairness  conditions  have  been  satisfied  can  we  ask  for  the  curve  which 
minimizes  the  deviation.  The  problem  would  be  less  complicated  were  it 
a  curve-fitting  problem  and  not  a  curve-fairing  problem. 

Linear  Programming 

One  computer  routine  which  can  be  very  satisfactorily  adapted  to  both 
curve-fitting  and  curve-fairing  is  entitled  "Linear  Programming."  The 
general  linear  programming  problem  can  be  described  as  follows: 

Given  a  set  of  M  linear  inequalities  or  equations  in  r 
variables,  we  wish  to  find  non-negative  values  of  these  variables 
which  will  satisfy  the  constraints  and  maximize  or  minimize  some 
linear  function  of  the  variables.* 

The  application  of  linear  programming  to  curve  fitting  stems  from  the 


For  a  detailed  definition  of  Linear  Programming  refer  to  G.  Hadley, 
Linear  Programming.  Addison-Wesley  Publishing  Co.,  1962 
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following  use  of  linear  programming : 

Linear  programming  offers  a  method  of  obtaining  an  optimal 
solution  to  a  set  of  linear  equations  which  has  an  infinite 
number  of  solutions  (i.e.,  a  procedure  for  determining  an 
optimal  element  from  a  set  with  an  Infinite  number  of 
e lements ) . 


The  latter  interpretation  of  linear  programming  indicates  that  the 
problem  of  fitting  a  curve  to  a  set  of  points  can  be  stated  as  a 
linear  programming  problem  by  requiring  the  linear  programming  code 
to  select  from  the  infinite  number  of  curves  in  the  neighborhood  of 
the  data  points,  the  one  which  comes  closest  to  these  points  (i.e., 
select  the  curve  which  satisfies  all  other  fairness  requirements  and 
minimizes  the  deviation  between  curve  and  points). 


Point  Smoothing  and  Second  Derivative  Requirement 

The  data-points  which  serve  as  input  to  the  lofting  system  (mathematical 
of  manual)  are  traditional  -  offsets  scaled  from  a  preliminary  design 
drawing  prepared  by  a  naval  architect.  This  drawing  is  approximately 
one-fiftieth  the  size  of  the  actual  ship.  The  smallness  of  this  drawing, 
although  being  beneficial  to  the  naval  architect,  invariably  causes 
small  reading  errors  in  the  scaled  offsets.  There  are  many  opportunities 
for  errors  aside  from  those  described  above  to  enter  into  the  data.  These 
errors  are  usually  quite  large  and  result  from  human  frailty.  Thtese 
errors  could  be  typically: 

(1)  An  error  resulting  from  using  the  wrong  scale  or  incorrectly 
using  the  right  one, 

(2)  An  error  due  to  transcribing  incorrectly, 

(3)  A  writing  error,  such  as  laying  the  point  down  incorrectly  on 
the  mold  loft  floor,  or  correspondingly  mispunching  a  data  card 
for  a  computer. 


1C  is  highly  probable  Chat  a  few  points  will  be  totally  In  error  and 
will  not  contain  any  Information  of  value.  Suparlmpoaed  on  these  and 
the  rest  of  the  points  will  be  small,  random  errors  due  to  reading 
Inaccuracy . 

Our  first  problem  is  then  to  separate  the  points  having  errors  which 
do  not  at  all  reflect  the  naval  architect's  Intention  from  the  rest. 
In  the  mold  loft  this  operation  is  done  by  eye.  Before  the  loftsman 
proceeds  with  putting  fair  lines  through  the  preliminary  offsets,  he 
checks  them  to  see  if  they  seam  to  fit  in  well  with  one  another. 

Those  points  which  seem  to  be  in  error  are  rechecked  with  the  lines 
plan  and  perhaps  even  re-scaled.  However,  this  is  a  liberty  not 
usually  afforded  the  computers  in  the  modern  approach;  therefore,  it 
is  crucial  that  this  process  of  rejection  of  bad  points  be  executed 
if  the  naval  architect's  intentions  are  to  be  preserved  in  a  later 
computation  with  these  points  (e.g.,  fitting  a  fair  surface  to  them). 

A  fair  curve  does  not  have  too  many  or  too  closaly-spacal  inflection 
points.  One  can  easily  verify  for  himself  that  two,  closely-spaced 
inflection  points  in  a  curve  causes  it  to  have  a  "bump"  in  that 
region.  Such  bumps  are  not  considered  fair  by  a  naval  architect.  In 
a  similar  way,  three  closely-spaced  Inflection  points  correspond  to 
two  "bumps."  Thus,  inflection  points  need  to  be  relatively  isolated 
from  one  another  to  prevent  the  occurrence  of  "bumps." 
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Consider  the  set  of  three  points  shown  in  Fig.  2.  Suppose  any 
continuous  curve  which  has  continuous  first  and  second  derivatives 
passes  through  these  points.  Since,  as  shown,  the  middle  point  (2) 
lies  above  the  line  connecting  the  two  other  points,  (1)  and  (3), 
there  must  be  some  portion  of  the  curve  with  negative  curvature. 
This  can  easily  be  deduced  by  repeated  use  of  the  Theorem  of  the 
Mean  Value  and  the  fact  that  the  curve  is  required  to  have  a 
continuous  second  derivative.  One  can  present  this  result  in  a 

) 

slightly  different  way.  The  "second  difference"  of  the  given 
three  points  is  given  by: 


(1)  r2 


VY2  VY1 

X3-X2  X2,X1 


For  convenience  we  refer  to  this  as  the  second  difference 
evaluated  at  point  (2),  even  though  it  obviously  depends  on  points 
(1)  and  (3)  also.  This  second  difference  has  the  property  of 
being  negative  when  (2)  lies  above  the  line  connecting  (1)  and 
(3) 1  and  positive  When  (2)  lies  below  this  line.  In  Fig,  2  the 
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second  difference  evaluated  at  (2)  la  negative  since  (2)  lies  above 

the  line.  Finally,  we  arrive  at  the  well-known  result  for  a  curve 
* 

satisfying  I-c.  There  must  be  a  portion  of  it  in  an  interval  which 
has  curvature  of  the  same  sign  as  the  second  difference  computed 
for  that  interval, 

It  can  be  seen  that  unlike  aigna  of  aecond  differencea  at  two 
neighboring  polnta  lndicatea  an  inflection  point  in  the  Interval 
between  theae  two  pointa.  We  can  now  develop  a  criterion  for 
detecting  pointa  with  large  errora  in  them.  If,  in  any  aet  of 
three  consecutive  pointa,  tha  aecond  difference  changea  aigna 
twice  we  have  a  "bad"  point.  The  "bad  point"  can  now  be  moved 
an  appropriate  amount  requiring  it  to  be  compatible  with  ita 
neighboring  points  (i.e.,  not  cause  any  "bumps"). 

After  the  points  have  been  "smoothed"  (large  errors  removed)  a 
requirement  that  the  second  derivative  and  second  difference  be 
of  the  same  sign  at  all  data  points  will  guarantee  the  analytic- 
spline  to  have  the  correct  amount  as  well  as  correct  spacing  of 
inflection  points. 

The  Coordinate  Axes 

The  coordinate  system  employed  is  indicated  in  Fig.  3,  below. 


Z 


Fig.  3  The  Coordinate  Axes 
'k 

See  Page  II -2 
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The  equation  for  the  surface  of  the  ship  will  be  a  function  Y  *  F(X,Z) 
The  two-dimensional  curves  obtained  by  holding  Z  constant  (Y  *  f(X)) 
will  be  referred  to  as  waterlines.  The  two-dimensional  curves  obtained 
by  holding  X  constant  (Y  -  g(Z))  will  be  referred  to  as  stations. 


Notation 

The  notatio  1  employed  offers  a  means  of  representing  the  analytic- 
spline  by  a  single  continuous  equation,  which  has  continuous  first 
and  second  derivatives  and  a  discontinuous  third  derivative.  The 
discontinuities  of  the  third  derivative  occur  at  pre-selected  points. 

In  essence,  this  equation  consists  of  individual  cubics  between 
points  of  discontinuity  of  the  third  derivative. 

The  analytic  spline  is  represented  by: 

Y(X)  -  Aq  +  AXX  +  A2X2  +  A3X3  +  A4(X-ax)J  +  ...  +  Vl^Vl*#  * 


a^  . . .  ,  a^_^  are  the  abscissae  of  the  points  of  dis¬ 

continuity  of  the  third  derivative  and  should  be  chosen 
to  coincide  with  the  abscissae  of  the  original  data  points. 


In  the  above  representation 

CX-a)3 

Or  equivalently 

(X-.)3 


{(X-a)3  if  X  >  a 
0  if  X  4  a 

*(X-a)3  +  %  |  (X-a)3 


♦This  notation  was  originally  introduced  by  Dr.  Feodor  Theilheimer 


3 

From  the  above  definition  it  is  obvious  that  (X-a will  never  be 
negative 

The  above  notation  is  very  direct. 

As  an  example,  consider  the  following  analytic-spline 

(2')  Y(X)  -  -1  +  3X  -  3X2  +  2X3  -  (X-l)J 

For  X  at  1  (2')  gives  the  cubic 

Yi(X)  -  -1  +  3X  -  3X2  +  2X3 

For  X>  1  (2')  gives  the  cubic 

YU(X)  -  X3 

It  can  immediately  be  seen  that  Y^  and  Y^  have  equal  ordinate 
values,  first  derivatives,  and  second  derivatives  at  X  »  1,  but 
unequal  third  derivatives. 
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B.  TWO-DIMENSIONAL  FAIRING 


As  previously  stated,  a  requirement  that  the  second  derivative  and  second 
difference  be  of  the  same  sign  will  guarantee  compatibility  of  inflection 
points . 

This  constraint  is  expressed  mathematically  as: 

Require  that 

(3)  r^.  •  Y"(Xi)  k  0  i 

Where 

r^  =  second  difference  at  i^1  data 
Y"(X^)  =  second  derivative  at  iC^  data 

Deviation  Requirement 

The  general  analytic-spline  Y(X)  is  given  by  (2)  . 

The  original  offsets  (ordinate  values  of  original  N  points)  are 
denoted  by  i  *  -1  ,  . . .  ,  N 

The  maximum  ordinate-deviation  between  curve  and  offset  is  noted  byX  (we 
take  maximum  to  mean  largest  magnitude  and  not  maximum  in  an  algebraic 
sense.) 

It  readily  follows  that 

(4)  |  Y(Xt)  -  yJS  X  for  i  -  1  ,  . . .  ,  N 

The  above  expression  simply  says  that  the  absolute  values  of  the  ordinate 
deviations  at  the  N  original  points  are  all  less  than  or  equal  to  X , 
(i.e.,  Xis  the  maximum  deviation). 

X  has  a  distinct  value  for  each  specific  analytic-spline  Y(X). 

We  now  ask  the  linear  programming  code  to  give  us  the  analytic-spline 
which  satisfies  fairness  constraint  (3)  and  has  the  smallest  value 
for  X  (i.e.,  minimum-maximum  deviation). 

*Since  the  second  difference  is  not  defined  for  the  first  and  last  points, 
the  second  difference  at  the  second  and  next  to  last  points  will  be  used 
for  these  two  points  respectively. 


=  1 . N 

point* 
po  int 
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Formulation  of  Deviation  and  Fairness  Constraints  as 
Linear  Programming  Problem 


The  curve  fitting  problem  can  now  be  formulated  as  a  linear  program* 
mlng  problem  as  follows! 

Minimize  Xsubject  to  the  constraints 


(3) 

rt  ‘  Y"(X1)i  0 

(4) 

lYtt^  -  Yt|<  X 

Since  linear  programming  codes  handle  linear  constraints  only, 
constraint  (4)  is  equivalently  expressed  by  the  following  pair  of 
linear  inequalities: 

Y(X1)i 

X+  Y(Xt)  i 


"Y, 


i  -  1  .....  N 


There  are  now  three  linear  constraints  per  point  the  analytic- 
spline  must  satisfy.  This  problem  readily  lends  Itself  to  linear 
programming  routines.  The  problem  is  stated  as: 

Minimize X subject  to  the  following  linear  constraints: 


X-  Y(Xt)fc  -Y 

X+  Y(X1)^  Y£ 

rt’  Y"(X1)i  0 

- 


N 


This  formulation  requires  N+3  non-slack  variables  and  3N  con¬ 
straints,  (  N  is  the  number  of  original  data  points). 


1.0. 0-1 
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Actual  Presentation 


The  actual  presentation  of  the  problem  requires  only  slight 
modification.  The  Inequalities  are  transformed  to  equalities  by 
the  subtraction  of  a  positive  quantity  (called  slack  variable) 
from  each  Inequality.  Since  present  linear  programming  codes 
accept  positive  variables  only,  the  variables  (coefficients  In 
equation  (2))  are  represented  by  the  difference  of  two  non¬ 
negative  variables,  l.e., 

A  *  A'  -  A" 

1  11 

A'  A'*  fc  0  V. 

The  formulation  of  the  linear  programming  problem  now  becomes: 

Minimize  Xsubject  to  the  following  constraints: 

X-  (Aq  -  Aq)  -  (A[  -  A,pX1  -  (A£  A”)xJ  -  (A^  -  A«)xJ  -  . . . 

(AN+1  ’  ^d-l^l-  '“n-I*#  "  S1 

X+  (A^  -  Aq)  +  (A^  -  A,')X1  +  (A^  -  A'pX^  +  (A^  -  A'^X^  +  ... 

+  (an+i  “  •an+i^(xi  "  ~  Pi 

*1  '  [2(*J  -  *5)  +  -  jqW,  +  ...  +  «<»-;-  Vl>«l  •  ■» 

"  Q1 


1  ■  1  ,  ...  ,  N 


This  difficulty  can  be  eliminated,  as  is  shown  on  Page  III-23, 
Section  entitled  "Dual  Formulation." 


) 

In  thi*  way  an  equation  of  the  form: 

Y(X)  -  A0  +  AjX  +  A2X2  +  A3X3  -  A4(X-«1)J  +  . . .  +  A^ 

will  be  constructed. 

The  maximum  deviation  between  curve  and  data  points  will  be  given 
by  the  value  of  Xln  the  final  solution. 


k  c. 


(5) 


THREE-DIMENSIONAL  FAIRING 

Extension  of  Notation  to  Thrce,«Oimenslona 


In  attempting  to  fit  a  surface  to  a  three-dimensional  set  of  points, 
the  same  procedure  as  was  used  in  the  two-dimensional  case  is  employed. 


Equation  (2)  is  replaced  by  its  three-dimensional  extension  Y(X,Z)  : 
Y(X,Z)  -  C0(X)  +  C^XjZ  +  c2(x)z2  +  c3(x)z3  +  c4(x)(z-b 

+  -  + 


Where 


C1<X)  -  c±  q  +  clflX+  ct  2X  +  c±  +  ...  +  c1>N+1(*“aN_1) 


i  -  0  ,  ...  ,  Mfl 


a  ■  X-Coordinnte  of  station 


P 

b  ■  Z-Coordinate  of  waterline 

q 

(5)  represents  a  surface  of  N  stations  and  M  waterlines. 


3 

# 


1.0. 0-1 


III-13 


An  expansion  of  Equation  (5)  gives: 


Y(X,Z)  -  I  cQ0  +  c01x  +  C02X  +  C03X  +  •••  +  co,N+l 


+  C10  +  CliX  +  C12X  +  C13X  +  •"  +  C1,N+1 


+  c20  +  C21X+  C22X  +  C23X  +  •“  +  C2,N+1 


[C30  +  C31 


X  +  C32X  +  e33X  +  +  c3,n+1 


+  lc40  +  CA1X  +  C42X  +  C43X  +  +  C4,N+1 


«-*»-!>#] 
<*-%-!>#] 2 
«-*.!>#]  * 
ff-N-i’j] 23 


+  pH-1,0  +  cMfl,lX  +  ^H-l^*  +  Sh-I.S* 

+  "*  +  CMfl,N+l(X-aN-l)#]  (Z-bM-l># 


We  now  have  an  expression  that  consists  of  individual  cubics  in  the 
coordinate  planes  between  any  two  stations  or  waterlines.  This  surface 
equation  will  have  (N+2)  (M+2)  coefficients  for  N  stations  and  M 

water lines. 
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Requirements 

The  deviation  is  again  required  to  be  less  than  or  equal  to 

Y^Zj)  -  X  i  -  1  .....  N 

j  *  1  ,  •••  i  M 

Two  second  difference  constraints  are  now  required  in  order  to 
guarantee  fairness:  one  in  the  X-direction  and  one  in  the  Z-direction. 

FORMULATION  AS  LINEAR  PROGRAMMING  PROBLEM 


( 


( 


The  problem  is  now  formulated  as  a  linear  programming  problem  as 
follows: 


Minimize  X subject  to  the  following  constraints: 

(b) 

X-  Y(X1,Zj)i  -Ytjj 

(bb) 

X+Ki.jJjrt  Yt>j 

i  *  l  ,  ...  ,  N 

(bbb) 

(bbbb) 


> 


*  rw'to 
j  •  rcx^z^xx  2  o 

second  difference  in  Z-direction  at 

second  difference  in  X-direction  at 


J  -  1 . M 

JCk  point  on  i^ 
i**1  point  on 


station 


waterline 
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Y"  (X,  Z.)  indicates  the  second  derivative  of 
i,  J  aa 

to  Z  (X  held  constant)  at  J6*1  point  on  iC^ 


Y(X,Z)  with  respect 
station. 


>Z.)  indicates  the  second  derivative  of 
'  i  j  xx 

to  X  (Z  held  constant)  at  i^  point  on 


Y(X,Z)  with  respect 
waterline. 


This  formulation  requires  4R  constraints 
(R  ■  number  of  original  points.) 


The  maximum  deviation  will,  as  in  the  two-dimensional  case,  be  given 
by  the  value  of^vin  the  final  solution. 
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D.  DOUBLE  SPLINE 


(i) 


On  Page  III-7  the  analytic  spline  is  represented  by 


Y(X)  -  AO  +  AjX  +  A2X2  +  A3X3  +  A4(X- 


a!>3  + 


Vl(X‘ 


N-r 


where  a,  ,  . . .  ,  a„  ,  are.  the  abscissae  of  the  original  input  points. 
With  the  use  of  the  dual  formulation  (Page  III-23)  the  number  of 
coefficients  in  Equation  (1)  becomes  the  number  of  deviation  con¬ 
straints  in  the  dual  problem.  Since  the  size  of  this  number  is  a 
very  critical  factor  .in  the  length  of  time  it. takes  to  solve  the  LP 
problem,  we  may  ask  the  question,  "What  effect  will  the  elimination 
of  some  of-. these  coefficients  have  to  the  final  solution."  In  other 
words,  can  we  effort  to  reduce  the  computer-,  time  at  the  expense  of 
possibly  inducing  a  larger  deviation  between  curve  and  offsets. 

One  way  of  systematically  reducing  the  number  of  coefficients  would 

3 

be  to  add  an  (X-a^^  term  only  at  every  second  data  point, 

(i.e.,  if  a1  ,  ...  ,  a  are  the  original  dat.  points,  add  (X-  .  3 

1  n  V# 

terms  if,  and  only  if  i  ■  even  integer.  Employing  .this  procedure 
■we  would  add  terms  only  at  a^  . a^  ,  . ,  a2^;# 

Motivation  for  this  procedure  can  be  extracted  from  various  policies 
.  the. graphical  or  manual  loft  uses.  Clearly,  the  graphical  loft  does 
not  put  a  weight,  or  "duck,"  at  every  point  and  generally  agrees  the 
curve  with  the  fewest  "ducks"  is  the  most  desirable,  provided  of 
course  that  the  deviation  is  reasonable'. 


It  was  previously  shown  . that  each  coefficient  Ay  for  J  >4  in 
Equation  (i)  corresponds,  to  a  point  where  a  duck  would  be  used  in 
the.  graphical  loft.  By  eliminating  every  second  coefficient  we  have 
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in  effect  eliminated  every  second  duck 


(ID 


In  practice,  it  has  been  shown  that  the  elimination  of  these 
coefficients  not  only  significantly  reduces  the  computing  time,  but 
also  does  not  reduce  the  quality  of  the  curve.  In  many  instances 
the  quality  of  the  curve  was  increased.  A  reasonable  explanation 
of  this  could  be  that  this  procedure  offers  fewer  places  to  have 
inflection  points  or  zero  second  derivatives  (reverse  curvature 
or  flat  splots).  This  procedure  of  joining  three  points  with  a 
single  cubic  is  referred  to  as  "double  splining." 

If  we  were  to  "double  spline,"  Equation  (1)  would  be  replaced  with 
its  counterpart: 

Y(X)  -  AQ  +  AXX  +  A2X2  +  A3X3  +  A4(X-ap#3  +  Afi(X-a4)#3  + 

Where  Aj^  ,  . . .  ,  aR  are  the  original  input  points.  The  only  change 
to  the  original  formulation,  aside  from  the  use  of  this  different 
curve  will  be  that  a  second  derivative  requirement  will  only  be 
needed  at  the  points  where  an  (X-a^)^  term  Is  added.  Second 
derivative  requirements  at  the  other  points  would  be  redundant.  As 
before,  a  deviation  constraint  is  required  at  all  original  points. 
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B.  ALTERNATE  REPRESENTATION  -  SUM  OF  THE  DEVIATIONS 


The  use  of  linear  programming  as  a  curve  fairing  tool  has  suggested 

an  alternate  fit  to  the  one  just  presented.  This  is  the  "sum  of 

the  deviations"  problem.  In  this  formulation,  the  sum  of  the 

absolute  values  of  the  deviations  between  the  surface  function  and 

n 

the  original  data  are  minimized  (i.e.,  minimize  r  |Y(X^)-YjJ  ) 

This  is  made  possible  by  that  characteristic  of  linear  programming 
which  requires  all  variables  to  be  non-negative. 

Assume  a  relationship  between  the  surface  equation  and  the  set  of 
original  offsets  such  that  the  value  of  the  equation  at  each  of  the 
original  offsets  is  less  than,  or  equal  to  the  original  offset. 

That  is, 

1 ' 1 . » 

J  *  1  .....  M 

We  can  change  this  inequality  to  equality  by  adding  a  dummy  variable 
as  follows  (note  that  the  dummy  variable  must  be  non-negative): 

i  -  1 . N 

J  -  1 . M 


y<Vj>  +  si,j '  yi.j 


In  linear  programming,  inequalities  are  always  changed  in  this 
manner  to  obtain  a  full  set  of  simultaneous  linear  equations.  The 
additional  variables  (S,  ,)  are  called  slack  variables.  As  previously 
stated,  in  any  feasible  solution  in  linear  programming,  all  of  the 
variables  in  the  solution  must  be  non-negative.  Therefore  in  any 
optimal,  feasible  solution,  all  of  the  slacks  which  remain  in  the 
solution  will  either  be  at  a  zero  level  or  positive..  By  requiring 
a  very  high  cost  on  the  slacks  when  compared  to  that  on  the 
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coefficients,  the  sun  of  the  slacks  (positive  deviations)  is 
minimized.  This  fairing  problem  is  formulated  as  follows: 


Minimize 


I 


subject  to  the  following 


\ 


Y(Xi,J>  +  Sl,j 


y”(x1  j)  ^  o 


Y"<Xi,j>xx*  0 


/ 


The  deviation  can  be  obtained  at  each  offset  by  reading  the  value 
of  the  slack  at  that  point.  If  the  slack  does  not  enter  into  the 
final  solution,  or  enters  at  a  zero  level,  the  surface  passes 
exactly  through  the  offset. 


In  practice  the  set  of  deviation  slacks,  plus  the  set  of  slacks 
for  changing  the  curvature  constraints  to  equalities,  which  together 
form  an  identity  matrix,  are  used  as  the  first  basic  feasible 
solution  allowing  the  optimization  phase  to  begin  at  once. 


A  modification  of  this  method  has  been  made  utilizing  the  following 
formulation  (this  is  done  so  the  deviations  will  be  unrestricted  in 
sign,  which  in  turn  should  give  a  closer  fit): 

Minimize ^  Si  j  +  ^i  j  8v*Ject  t0  following 


+  *i.j  -  "i.i 


>.  0 

,>  0 


I  i  -  1  ,  ...  N 
/  j  -  1  ,  ...  M 
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Here  the  slacks  S.  .  and  Q,  .  are  both  positive,  however,  only 
one  of  the  pair  may  enter  into  the  solution  since  they  are  linearly 

•ff 

dependent,  (i.e.,  will  cause  a  singular  matrix) 

This  formulation  allows  the  surface  to  pass  over  some  offsets  and 

under  others.  The  first  basic  feasible  solution  again  consists 

of  the  set  of  slacks  S.  ,  plus  those  of  the  curvature  constraints. 

J 


( 


only  one  of 


Positive  deviation 
Negative  deviation 

Total  deviation,  as  stated  above,  in  practice, 

S,  ,  or  Q,  will  enter  into  the  solution. 

J-*!  J-jJ 


/ 

V 


( 


F.  DUAL  FORMULATION 

jc 

It  can  be  shown  that  if  we  write  the  primal  problem  in  the  form: 
r 

ZdiJ*XJ  di  i  -  1 . k 

J-l 

r 

Edi,j  *  xj  -di  i  -  k  + 1 , ... ,  s , 

.  j-i 

Xj  >  0  j  -  1  ,  . . .  ,  r  MAX  Z  -  cm 

then  the  dual  problem  is 


D'w  >  c'  ,  MIN  Z 

■  d'w 

unrestricted 

i  *  1  j  *• «  9  k 

w^  "2.  0 

■  i  ®  lc  ^  1  9  •  •  •  j  8 

That  is,  wherever  iC^  constraint  in  the  primal  is  an  equality,  the 

th  • 

i  dual  variable  is  unrestricted  with  regard  to  sign.  This  property 

has  considerable  importance  in  our  problem,  because  the  coefficients 

of  the  surface  equation  are,  by  nature,  unrestricted.  The  dual 

formulation  offers  a  way  of  eliminating  the  representation  of  each 

coefficient  as  the  difference  of  two  non-negative  variables,  and  this 

greatly  reduces  the  size  of  the  problem 


♦Refer  to  G.  Hadley,  Linear  Programming . p  236  (Addison-Wes ley 
Publishing  Co.,  1962 


(' 
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We  therefore  restate  the  problem  as  follows: 


Minimize  X  subject  to  the  following  constraints 


■  -yi.j 
x+,(xi,i)  ‘  yi,j 
ri,j  *  2  o 

*1.J  •  V"<Xl.)>»  s  0 


\ 


/ 


•  •  • 


,  N 
,  M 


If  we  present  this  formulation  to  the  linear  program  in  the  dual 
form,  the  variables  representing  the  deviation  constraints 
(coefficients  of  the  surface  equation)  will  be  unrestricted  with 
regard  to  sign. 
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Section  IV 


FORMULATION  OF  END  CONDITIONS 


The  previous  formulation  of  the  surface  fitting  problem  can  handle 
only  rectangular  arrays  of  points,  that  is,  every  waterline  must  have 
the  same  number  of  stations;  and  conversely,  every  station  must  have 
the  same  number  of  waterlines.  This  condition,  unfortunately,  does  not 
exist  when  we  consider  the  bow  and  stern  of  a  ship.  The  lower  waterlines 
may  have  fewer  stations  than  the  higher  ones.  The  bow  and  stern  areas 
also  produce  other  problems  which  we  shall  refer  to  collectively  as 
"end  conditions." 

If,  in  the  equation  for  the  surface  of  the  ship  (i.e.,  Y  ■  f(X,Z)),  we 
replace  Y  with  zero  we  obtain  the  expression  for  the  profile  of  the 
ship  (i.e.,  f(x,z)  *  0). 

This  profile  is  specified  in  the  original  lines  drawings  of  the  ship 
and  at  times  may  be  very  complicated.  This  profile  must  also  be  "fair." 
The  present  formulation  offers  no  means  for  guaranteeing  this.  In 
attempting  to  solve  this  problem  we  can  be  guided  by  the  procedure  used 
in  the  manual  loft. 

In  manual  lofting  the  profile  is  "faired"  as  an  independent  unit  of  the 
ship.  The  ship  is  then  required  to  pass  through  this  previously  "faired" 
profile.  In  this  way  the  ship  is  guaranteed  to  have  the  correct  prrofile, 
no  matter  how  complicated  it  may  be. 

We  shall  attempt  to  construct  an  analogous  procedure  in  the  mathematical 
loft.  To  begin  with,  let  us  show  that  it  is  not  sufficient  to  require 
that  the  original  waterlines  pass  through  the  profile,  without  any 


additional  restrictions  on  the  surface  between  the  original  waterlines. 


Fig.  1 


In  Fig.  1,  A,  B,  C,  and  D  are  waterlines  passing  through  the  profile  at 
Points  a,  b,  c,  and  d,  respectively.  The  contour  of  the  profile  desired 
i  shown  in  Fig.  1-a;  but  in  the  surface  fitting  problem,  a  requirement 
that  Waterlines  A,  B,  C,  and  0  pass  through  Points  a,  b,  c,  and  d, 
respectively,  without  any  additional  requirements  on  the  profile  of  the 
surface,  may  very  possibly  result  in  the  Profile  shown  in  Fig.  1-b.  At 
any  rate,  we  cannot  guarantee  that  this  will  not  occur. 

Consider  the  expression  for  the  surface  of  the  ship: 

Y  -  f(X,Z) 

We  wish  to  transform  this  expression  so  that  Y  will  be  equal  to  zero  at 
the  profile  of  the  ship.  Toward  this  end  we  construct 

Y  -  f (X,Z)  •  T(X,Z) 

Where  T(X,Z)  has  the  property  that  it  is  equal  to  zero  if  the  point 
(X,Z)  in  the  Y  ■  0  plane  lies  on  the  profile. 

Let  T(X,Z)  -  |l  -  [l+f(Z)-x]J| 

f(Z)  is  the  "faired"  spline  curve  representation  of  the  profile.  (Clearly, 
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since  Y  ■  0  at  the  profile,  the  profile  will  be  an  expression  In  X 
and  z  only).  Since  X  -  f(z)  Is  the  expression  for  the  profile,  for 

a  given  waterline  -  say  Z  -  ZQ  ,  f(ZQ)  will  be  the  X  value  corresponding 
to  the  intersection  of  the  profile  with  this  waterline. 

When  X  -  f(Z^)  ,  i.e.,  when  Point  (X,Zq)  lies  on  the  profile 

■  [i-  [>  -  £<zo)]#J,  -  l-1  *  0 

SO  Y  -  f(X,Z0)  T(X,Z0  )  -  0 

But,  if  X2f(Z0)  +  1  (i.e.,  X  is  one  unit  away  or  farther  from 

the  profile) 


T(X,Zq) 

"M 

.  +  f(2„)  -  -  l]j]#  - 

Where 

here  we 

define 

such  that 

K*  "  (°  ^  A#:S° 

]A  if  A#>0 

So 

T(X,Zq) 

-  !-A# 

with  A^  <  0  ,  becomes 

X(X,ZQ) 

o 

1 

H 

■ 

-  1  for  X  i  f  (Z0)  +  1 

and 

Y(X,Z) 

-  f (X,Z)  T(X,Z)  -  f (X,Z) 
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This  implies  that  T(X,Z)  is  effective  only  between  the  profile  and  one 
unit  away  from  it,  also 


0  £  T(X,Z)  1  for  all  X,Z 

We  raise  ^1  +  f(Z)  -  xj  to  the  third  power  so  we  may  have  continuous 
first  and  second  derivatives  as  well  as  continuous  curve  at  X  ■  f(Z)  +  1  . 

If  we  consider  stations  instead  of  waterlines,  we  see  that  similar 
results  will  be  obtained. 


The  expression  for  the  surface  now  becomes 

Y(X,Z)  -  f (X,Z)  £l 

f(Z) 


previously  faired  profile. 


If  it  is  desired  to  have  infinite  slope  where  the  waterline  (station) 
intersects  this  profile,  we  merely  replace 


T(X,Z) 


1  -  [l  +  f (Z)  -  x]j 

•  - 


X 


T(X,Z) 


with 


where  0  <  P  <  1 


The  interval  in  which  T(X,Z)  is  effective  does  not  have  to  be  from  the 
profile  to  the  profile  plus  1  unit,  but  can  be  from  the  profile  to  the 
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profile  plus  any  arbitrary  distance 


D  ,  by  using 


v. 


f 


V 


( 


For  X  - 

T(X,Z)  - 


f(Z)  +  D 


,  [P  ±  fZ  -  f(Z)  -  D 
1  D 


T(X,Z)  -  [l]P  -  1 

It  can  very  easily  be  verified  that 

I(X,Z)  «  1  if  X>  f(Z)  +  D 


or 


Y(X,Z)  -  f(X,Z) 

f(X,Z)  is  the  spline-surface  representation  developed  in  a  previous 
section.  Y(X,Z)  is  now  guaranteed  to  have  the  correct  end  profile  and 
becomes  a  regular  spline  curve  D  units  away  from  the  profile. 


In  Fig.  2,  the  cross-hatched  area  represents  the  portion  of  the  surface 
in  which  T(X,Z)  is  effective. 
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Section  V 


LOFT  INFORMATION 


From  the  mathematical  fairing  program  has  come  a  surface  equation  which 
describes  the  desired  hull  form.  It  now  becomes  necessary  to  extract 
from  this  surface  equation  the  data  required  for  parts  definition. 

THE  FINAL  LINES  DRAWING 

One  of  the  first  and  more  important  types  of  information  required 
from  a  loft  is,  of  course,  a  final  body  plan.  The  data  for  the 
preparation  of  such  body  plan  is  easily  obtained  from  the  mathe¬ 
matical  ship  definition.  Substituting  any  specific  frame 
location  X  into  the  hull  surface  equation  Y  ■  f(X,Z),  and 
solving  this  expression  for  various  values  of  Z  ,  will  provide 
offsets  Y  describing  the  frame  of  interest.  The  values  for 
the  coordinates  of  these  points  may  be  obtained  at  any  desired 
scale.  If  an  automatic  plotting  machine  is  available,  such  as 
the  Coordinatograph,  the  body  plan  may  be  scribed  quickly  and 
with  great  accuracy. 

Other  information  for  the  lines  drawing  may  be  obtained  in  the 
same  manner.  Water lines  may  be  described  by  substituting  a 
specific  value  of  Z  ,  into  the  surface  equation  Y  ■  f(X,Z)  , 
and  solving  the  equation  at  various  values  of  X  . 

The  coordinates  of  the  buttock  points  are  not  obtained  in  such 
straight-forward  fashion.  However,  using  a  computer,  a  technique 
such  as  Newton's  method  may  be  employed  to  intersect  the  hull 
surface  with  a  plane  parallel  to  the  center  plane  of  the  ship. 
Heights  of  the  intersection  of  the  plane  and  the  surface  may  be 


) 

obtained  at  many  points  along  the  length  of  the  ship.  Diagonal 
planes  can  be  handled  In  the  same  way,  Interaectlng  the  surface 
with  a  plane  at  a  constant  angle  to  the  hull  surface. 

.  SkEtL  PIATE  DEVEI0PMBKT 

To  develop  the  flat  pattern  to  which  a  plate  must  be  cut  In  order 
to  fit  properly  in  place  on  the  hull  when  formed  requires  several 
types  of  information: 

a.  Arc  lengths  along  intersections  of  coordinate  planes  (station 
or  waterline  planes)  and  the  hull  surface 

b.  Coordinates  of  points  along  the  sight  edges  on  the  hull 
surface 

c.  Arc  length  along  the  sight  edge  on  the  surface 

d.  True  angles  on  the  surface  between  coordinate  planes 

The  mathematics  for  developing  this  information  is  given  below. 


a.  Finding  arc  lengths  on  the  hull  surface  along  station  or 
waterline  planes: 

The  general  formula  for  arc  length  along  a  station  curve, 
X  -  constant  is 


6) 


When  X  is  constant  the  hull  surface  equation  becomes 
7)  Y  ■  A  +  BZ  +  CZ2  +  DZ3  +  . . . 
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Substituting  7)  into  6)  gives 


8) 


S 


(B  +  2CZ  +  3DZ2  + 


h 


dZ 


o 

which  can  be  solved  for  arc  length  between  any  two  heights  Z  . 

Similarly,  the  equation  for  arc  length  along  the  surface  in  a 
waterline  plane  is 


9) 


1  +  (B  +  2C  +  3DX2  + 


dX 


b.  Coordinates  of  any  point  along  sight  edges: 

To  establish  the  plate  sight  edges  mathematically  from  the 
hull  surface  equation  requires  a  function  Z  -  f(X)  which 
describes  the  trace  of  the  sight  edge  on  the  X  -  Z  plane. 
Assume  this  function  is  known  for  the  moment. 

The  equation  describing  the  hull  surface  is  of  the  form 

10)  Y  -  f(X,Z) 

Substituting  Z  ■  f(X)  into  (5) 

11)  Y  -  f(X) 

Given  any  value  of  X  on  Z  ■  f(X)  ,  Z  and  Y  are 
determined  by  equations  6)  and  8) 


( 
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The  method  for  obtaining  the  function  Z  ■  f(X)  le  to  uee 
the  preliminary  offaete  given  by  the  naval  architect.  Theae 
consist  of  the  height  above  the  baseline  for  each  sight  edge 
at  each  station.  The  set  of  heights  (Z.  ,  ...  ,  Z  )  for 

1  ID 

each  sight  edge  describes  the  trace  of  the  sight  edge  on 
the  centerline  plane  of  the  ship. 

To  establish  Z  “  f(X)  two  steps  are  necessary: 

1.  Since  the  heights  (Z^  ,  ...  ,  Z^)  are  given  only  to 
the  accuracy  obtained  from  scaling  the  drawings,  they 
must  be  smoothed  to  prevent  unwanted  inflections. 

2.  The  smoothed  points  can  then  be  fitted  with  a  spline 
curve  using  the  smaller  two-dimensional  linear  program¬ 
ming  formulation. 

It  may  be  possible  to  assume  that  in  the  equation  Z  -  f(X)  , 

Z  is  a  simple  function  of  X  such  as  a  parabola.  This 
would  require  considerable  analysis  and  perhaps  revision  of 
the  naval  architect's  method  of  establishing  sight  edges. 


c.  Arc  lengths  along  the  sight  edges  on  the  surface: 

The  sight  edge  has  now  been  defined  parametrically 

12)  X  -  X 

Z  -  f(X) 

V  -  f(X) 
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Substitution  into  the  general 
integrating  along  the  X  axis 


formula  for  arc  length  and 
give  8 


cor  +  oo) 


dX 


The  true  angle  between  coordinate  planes  on  the  hull  surface 

Substitution  of  the  directional  cosines  of  the  coordinate 
planes  into  the  law  of  cosines  gives  the  angle  between  these 
planes  as: 


LONGITUDINALS 

Many  times  the  position  of  longitudinals  around  the  hull  of  the 
ship  is  given  by  the  arc  length  distance  from  the  keel  or  other 
longitudinal . 


Fig.  4  Establishing  the  Position  of  Longitudinals 

That  is,  (Fig.  4)  given  S  and  Yq  ,  Zq  (which  may  be  the  keel 
line)  find  Y^  and  Z^  which  are  the  coordinates  of  the  position 
of  the  longitudinal. 

The  height  Z^  may  be  found  by  a  trial  and  error  solution  of 
Equation  3)  after  substituting  Zq  and  S  .  Y^  can  then  be 
found  by  substituting  the  station  location  X  ,  and  Z^  into 
the  surface  equation. 
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LONGITUDINALS 

Many  times  the  position  of  longitudinals  around  the  hull  of  the 
ship  is  given  by  the  arc  length  distance  from  the  keel  or  other 
longitudinal. 


Fig.  4  Establishing  the  Position  of  Longitudinals 

That  is,  (Fig.  4)  given  S  and  Y  ,  Z  (which  may  be  the  keel 

o  o 

line)  find  Yj^  and  Z^  which  are  the  coordinates  of  the  position 
of  the  longitudinal. 

The  height  Z^  may  be  found  by  a  trial  and  error  solution  of 
Equation  3)  after  substituting  ZQ  and  S  .  Y^  can  then  be 
found  by  substituting  the  station  location  X  ,  and  Z^  into 
the  surface  equation. 
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Section  I 


INTRODUCTION 


The  Mathematical  Lofting  Syatem  begina  with  offseta  and  data  from  a 
naval  archltect'a  preliminary  lines  plan.  This  data  is  first  checked 
to  be  sure  it  contains  no  extraneous  errors.  It  is  next  faired  by 
fitting  it  with  a  surface  equation  which  represents  the  molded  surface 
of  the  hull  and  which  has  certain  properties  consistent  with  "fairness." 
Then  data  necessary  for  building  the  actual  ship  is  extracted  from  this 
equation. 

Sections  II  through  IX  describe  the  operation  of  this  System.  The 
instructions  for  the  computer  programs  and  program  listings  are  con¬ 
tained  in  the  Appendixes. 

The  coordinate  system  used  throughout  is  shown  in  Fig.  1-1.  All  dimen¬ 
sions  used  in  this  system  are  in  feet  unless  expressly  stated  otherwise. 


f 


Fig.  1-1  Coordinate  System 


A  diagramatic  representation  of  the  Lofting  System  is  shown  in  Fig.  1-2. 

It  is  expected  that  anyone  using  this  System  will  be  familiar  with  the 
SFS  and  FORTRAN  programming  languages  for  the  UM-1620  computer.  Although 
the  programs  included  in  this  report  represent  a  lofting  system,  it 
must  be  remembered  the  purpose  of  the  programs  was  primarily  to  checkout 
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tha  development  work  on  the  mathematical  lofting  method. 

The  organization  of  thia  report  Into  Stapa  followa  tha  numbering  of 
the  boxai  in  Fig.  1-2.  Tha  flow  of  data  through  the  8yat*a  alao 
follows  thia  same  order. 
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Fig.  1-2  Mathematical  Lofting  System 
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Section  II 


PUPAUTION  OF  INPUT  DATA  (SUP  1) 


The  flow  of  deta  through  the  preparation  stage  is  shown  on  Pig.  II-l. 


Fig.  II-l  Flow  Diagram  for  Data  Preparation 


A.  OFF SITS 

The  left-hand  branch  of  Fig.  II-l  outlines  the  procedure  for  pre¬ 
paring  the  table  of  offsets  for  a  ship  to  be  lofted.  This  somewhat 
expanded  table  of  offsets  must  include  values  for  every  intersection 
of  station  and  waterline  shown  in  the  circumscribing  rectangle  of 
Fig.  II-2.  Although  the  values  of  the  offsets  in  the  shaded  portions 
are  not  important  and  usually  are  made  zero,  they  are  included  for 


PPOFILP  FUNCTION 
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Fig.  II-2 


mechanic#!  reason*.  The  offsets  above  the  deck,  or  knuckle  line, 
must  be  consistent  with  the  others  and  faired  to  a  degree  such 
that  the  second  differences  calculated  using  these  points  will  have 
the  proper  sign. 

The  offsets  are  punched  into  cards,  along  with  their  X  and  Z 
locations  on  the  hull.  These  cards  are  read  by  an  HM- 1620  com¬ 
puter  program  which  converts  them  from  feet  and  fractions  to  feet 
and  decimals,  and  "smooths"  them.*  The  principle  is  a  necessary 
part  of  the  lofting  system.  This  procedure  essentially  points  out 
any  offsets  which  are  obviously  in  error,  and  substitutes  a  value 
which  will  allow  a  fair  curve  to  come  near  it,  without  causing  dis¬ 
crepancies  at  other  points  (corrects  the  second  differences). 
Presently,  only  the  offsets  within  the  dotted  rectangle  of  Fig.  II-2 
may  be  smoothed  by  this  routine,  since  it  has  not  been  programmed 
to  accept  the  end  profiles.  The  remaining  points  can  be  examined 
as  noted  in  Step  2 -A.  After  gaining  experience  in  the  fairing 
procedure,  the  smoothing  program  can  often  be  bypassed  and  the 
data  fed  directly  into  the  matrix  generating  routines.  The  listing 
of  the  signs  of  the  second  differences  can  then  be  examined  for 
consistency.  More  information  about  this  can  be  found  in  Step  2-A. 

The  output  from  this  program  is  directly  usable  in  the  second  step 
of  the  lofting  system.  Operating  instructions,  card  formats,  and 
a  program  listing  are  included  in  Appendix  A. 


The  smoothing  process  is  described  in  Mathematical  Ship  lofting  - 
fart  1.  Theory.  Technical  Report  1.0. 0-1 
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b.  nanus 


la  addition  to  thn  smoothed  offsets,  functions  which  describe  the 
fore  end  aft  profiles  must  be  presented  to  the  matrix  generating 
routines  of  Step  2.  This  process  is  represented  bp  the  right-hand 
branch  of  Fig.  II- 1.  The  functions  representing  the  end  profiles 
are  obtained  from  the  end  profile  coordinates  using  the  approach  of 
Step  4. 

The  function  representing  the  profile  is  x  ■  f  (z)  as  shown  in 
*ig.  II-2.  The  origin  of  this  equation  arust  be  the  seme  as  the 
origin  of  the  surface  equation  for  the  hull,  end  must  be  ecaled 
the  seme  (see  Step  2).  The  data  points  for  describing  the  profile 
must  be  on  the  same  waterlines  as  are  used  for  the  surface  grid. 

The  coefficients  of  the  equations  for  the  profiles  are  punched  on 
cards  in  the  format  shown  in  Appendix  B  to  be  entered  directly 
into  the  Step  2  routines. 
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Section  III 


GENEBAIIMG  THE  LUOA1  P10GRAMMIBG  MATE IX  (STEP  2) 


The  matrix  for  the  linear  program  is  produced  by  a  set  of  1620  programs 
called  SMOG  (Surface  Matrix  from  Offsets  Generator).  These  routines 
accept  the  data  assembled  in  Step  1  and  produce  the  specific  matrix 
required  for  a  particular  surface,  including  the  auxiliary  information 
required  by  the  19  program  for  bookkeeping  purposes. 

There  are  four  versions  of  SMOG  (Appendix  B) : 

SMOG-1  produces  the  matrix  required  for  a  surface  problem,  without 
profile,  using  the  dual  formulation  and  >^fit. 

SMOG-2  produces  the  matrix  required  for  a  surface  problem,  with  profile, 
using  the  dual  formulation  and  Xfit. 

SMOG-3  produces  the  matrix  required  for  a  surface  problem,  with  profile, 
using  the  primal  formulation  and  sum  of  the  deviations  fit. 

SMOG-4  produces  a  matrix  for  a  2-dimensional  curve  using  the  primal 
formulation  and  sum  of  the  deviations  fit. 

Most  surface  problems  will  be  run  using  SMOG-1  or  2  .  The  first  three 
surface  versions  produce  their  output  in  the  SHAKE  Standard  UP  format, 
although  the  output  from  SMOG-3  must  be  rearranged  as  shown  in  the 
Appendix. 


The  general  operation  of  the  programs  is: 

(1)  lead  the  data  and,  if  necessary,  convert  the  offsets  from 
feet-fractions  to  feet-in-decimal  form 

(2)  Calculate  and,  if  required,  print  the  intersections,  if  any, 
of  each  station  and  waterline  with  the  profile 

(3)  Calculate  and,  if  desired,  print  values  and  signs  of  the  second 
differences  in  both  the  waterline  and  station  directions 
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(4)  Calculate  and  punch  the  1?  matrix,  consisting  of  deviation 

constraints,  curvature  constraints,  and  the  peripheral  informa¬ 
tion  required  bp  the  linear  program. 


4.  M4904L  SMOOTH IHG 

In  operation  (3)  above,  the  eigne  of  the  second  differences  ere 
found  and  can  be  printed.  The  second  differences  of  the  offsets 
are  calculated  using  a  general  formula,  therefore,  there  is  no 
requirement  that  the  stations  or  waterlines  chosen  to  describe  the 
surface  be  equally  spaced. 

If  the  offsets  have  not  previously  been  run  through  the  smoothing 
routine,  they  may  be  visually  checked  at  this  time.  That  is,  they 
may  be  compared  against  the  preliminary  lines  drawing  to  see  if 
they  actually  describe  a  curve  of  the  same  character  as  shown  on 
the  drawing.  A  second  check  is  to  investigate  the  consistency  of 
the  signs  of  the  second  differences.  For  instance,  if  a  single 
positive  second  difference  turns  up  in  the  middle  of  a  string  of 
negative  values,  the  point  should  be  carefully  checked  for  validity 
and,  if  necessary,  changed  a  small  amount  to  adjust  the  sign. 

There  ere  also  some  cases  when  the  correct  values  for  the  points 
will  cause  the  wrong  sign  to  be  calculated.  This  is  illustrated 
in  Fig.  III-l.  The  example  is  taken  from  Station  14  of  the  DIG -2 6 
class  ships.  Using  points  from  the  base  line,  4'  and  8'  waterlines, 
a  negative  sign  will  be  calculated  for  the  second  difference  at  the 
4'  waterline,  where  the  character  of  the  curve  suggests  a  positive 
value.  Simply  changing  the  sign  within  the  program  before  punching 
the  curvature  constraints  will  correct  this  situation. 


l.Q.Q-2 
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B.  SCALING 

All  of  the  linear  programming  routines  which  have  been  examined 

have  a  restriction  that  the  size  of  a  matrix  element  be  less  than 

the  absolute  value  of  10,000.  Since  some  of  our  matrix  elements 

3  3 

are  the  product  of  two  cubes  (X^  *  Zj  )  their  size  quickly 

becomes  large.  This  has  necessitated  a  limitation  on  the  size  of 
the  X  and  Z  coordinates.  It  can  be  readily  seen  that  this 
limitation  requires  the  product  of  the  greatest  X  and  Z  coordinates 
to  be  less  than  the  cube  root  of  10,000,  or  21.53.  In  practice,  the 
4x5  combination  is  most  common. 

The  coordinates  may  be  scaled  by  dividing  by  some  multiple  of  the 
station  or  waterline  spacing.  A  surface  five  stations  across,  for 
example,  might  have  a  station  spacing  of  25'  for  a  total  length  of 
100'  .  Division  by  the  station  spacing  gives  a  scaled  value  of  4  , 
or,  the  stations  are  one  unit  apart. 

The  most  it  is  presently  practical  to  scale  is  by  four  times  the 
minimum  station  or  waterline  spacing,  since  this  makes  the  minimum 
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spacing  one -fourth  unit.  The  smallest  matrix  element  then,  becomes 
(lA)^  x  (1/4)3  which.  If  the  largest  element  is  greater  than  9000, 
creates  an  extrems  range  in  the  number  else.  This,  of  course,  causes 
round-off  and  other  difficulties  in  solving  the  matrix.  The  maimer 
in  which  scaling  is  handled  by  the  SMOG  routines  is  noted  in 
Appendix  B. 

Another  important  consideration  for  reducing  maximum  nuafcer  size  is 
to  require  that  the  first  waterline  end  first  station  have  zero 
coordinates.  The  SMOG  programs  have  the  ability  to  translate  about 
the  X  or  Z  axis  so  the  origin  will  occur  at  0,0  . 


C.  SUKJACK  UmSWIATUW 

All  three  of  the  surface  matrix  generators  "double  spline"  the  surface 
in  the  X  direction.  That  is,  they  require  a  cubic  to  cover  two 
Intervals  instead  of  one,  as  the  "single  spline"  does.  The  cubic 
is  required  to  have  the  proper  curvature  condition  at  every  other 
station  and  is  required  to  come  as  close  as  possible  to  every  station. 
This,  as  can  be  seen  in  Tig.  III-2,  requires  that  every  surface  have 
an  odd  number  of  stations. 


Solid  lines  -  Double  spline  both  ways 
Combination,  solid  horizontal,  broken  vertical  -  double 
spline  one  way  only 

Tig.  I1I-2  Spaa  of  Cub lcs  on  Surface 
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In  the  Z  direction,  there  ie  e  choice  between  double  epline  end 
•ingle  spline.  The  double  spline  representation  requires  e  smeller 
metrix  then  the  single  spline  representation.  The  double  spline  in 
the  Z  direction  should  be  used  whenever  possible.  It  is,  however, 
sometimes  necessary  to  use  single  spline  in  order  to  gain  the 
necessary  flexibility  for  fitting  surfaces  with  large  variations 
in  curvature  and  slope. 

The  double  spline  representation  requires  the  data  to  contain  an 
odd  number  of  water lines. 
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Section  17 


USE  OF  LItiEAR  FKOGBAM  (STEF  3) 


The  linear  programming  routine  which  we  have  used  most  for  fairing, 

* 

is  called  LF-90.  It  uses  the  SHABE  standard  format  for  input  as  punched 
by  the  SMOG  programs.  Some  information  on  its  use  will  be  found  in 
Appendix  C. 

Another  IF  system  has  been  written  by  Bonner  &  Moore.  This  program 
is  faster  than  LP-90  for  very  dense  problems,  such  as  the  fairing  problem. 
It  is,  however,  somewhat  more  difficult  to  use.  This  system  also  accepts 
SHARE  standard  input  through  the  use  of  a  special  loading  routine. 

The  matrix,  which  was  produced  by  SMOG  on  the  IBM-1620,  must  be  con¬ 
verted  to  card/magnetic  tape  off-line  on  some  machine,  such  as  the 
IBM-1401,  before  its  use  with  either  IF. 

The  output  from  LF-90  is  a  magnetic  tape  that  contains  information  on 
the  problem  written  during  its  solution,  and  the  final  solution,  which 
is  the  set  of  coefficients  for  the  surface  equation.  By  monitoring  the 
output  tape  using  the  printer  on  the  1401,  this  final  solution  may  be 
punched  on  cards.  These  cards  may  then  be  used  as  direct  input  to 
Step  5. 


* 

LF-90  was  written  by  CEIR  Inc.,  Arlington,  Virginia.  It  can  be 
obtained  from  SHARE  by  sending  them  a  magnetic  tape  on  which  they  will 
copy  the  program.  It  is  also  available  at  any  service  bureau  of  CEIR. 

**Bonner  &  Moore  Associates,  Houston,  Texas.  This  system  may  be 
purchased  from  them  or  used  at  any  Service  Bureau  Corp,  Scientific 
Computing  Center. 
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Section  T 


MrormoM  or  longitusiiul  stkjctdui  nnnsscriovs  (snr  4) 


There  Is  one  type  o£  intersection  with  the  hull  surfece  that  requires 
additional  information.  This  is  the  type  which  is  shown  in  Fig.  V-l. 
Examples  are: 

Shell  plate  sight  edges 
Longitudinal  frames  and  stringers 
Seek  edges 
Knuckles 

The  centerplane  projection  of  these  intersections  is  a  fair  spline  curve. 


In  order  to  complete  the  definition  of  these  curves,  an  equation  for 
their  height  at  any  point  along  the  length  of  the  ship  must  be  found 
(a  *  f(x)).  This  function  is  entered  into  the  routines  explained  in 
Step  5,  along  with  the  hull  equation  to  determine  the  height  and  half¬ 
breadth  of  these  intersections  on  the  fair  hull.  The  procedure  for 


obtaining  these  functions  follows 


A.  DATA  PREPARATION 

The  date  for  each  curve  generally  consists  of  the  height  of  the 
curve  at  each  station  used  to  fair  the  hull  surface  (Fig.  V-l). 

This  data  can  be  entered  Into  the  two-dimensional  smoothing  routine, 
Appendix  A,  or  directly  Into  the  two-dimensional  SMOG  (SMOG  4)  which 
is  described  below. 

B.  GENERATION  OF  MATRIX 

The  data  from  above,  plus  certain  additional  information  as  noted 
in  Appendix  B,  are  entered  into  SMOG  4.  This  program  scales  the  data 
as  desired,  calculates  and  prints  the  second  differences,  and  punches 
the  matrix  for  the  1620  linear  program  given  below.  The  matrix  is 
punched  in  the  sum  of  the  deviations  formulation.  SMOG  assumes  the 
curve  goes  through  the  first  point. 

C.  USE  OF  LINEAR  PROGRAM 

The  linear  programming  code  used  for  two-dimensional  curves  on  the 

* 

1620  computer  was  obtained  from  the  1620  Users  Group.  The  code  is 
far  from  ideal,  since  it  contains  many  compromises  in  order  to  fit 
on  the  smallest  1620  available.  It  is,  however,  easily  obtainable 
free  of  charge.  There  is  a  proprietary  code  written  by  Bonner  & 
Moore  Associates,  which  may  be  more  efficient. 

D.  USE  OF  THE  EQUATION 

The  coefficients  obtained  from  the  LP  may  be  directly  entered  into 
the  next  phase  after  punching  on  cards  in  the  proper  format.  These, 
together  with  the  hull  equation,  enable  the  height  and  half-breadths 
of  the  intersections  to  be  obtained. 

*1620  Users  Group  Program  10.0.006 
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Section  VI 


CALCULAIIOH  OF  BASIC  BUi  OF  HULL  POXM  (SHF  5) 


The  determination  of  the  basic  date  necessary  to  check  the  surface 
equation  and  to  construct  the  ship  is  accomplished  by  a  FORTRAN  II 
program  for  the  1XM-1620.  This  program  is  called  GOBACK  (the  program 
reverses  the  fairing  procedure  and  "goes  back"  to  offsets).  GQBACK 
is  documented  in  Appendix  D. 


The  program  solves  the  surface  equation  for  the  following  kinds  of 
date: 

Offsets  along  any  waterline  or  station  (not  necessarily  the 
original  waterlines  or  stations)  at  any  interval  desired 

Slopes  and  curvatures  at  any  given  interval  along  any  water¬ 
line  or  station 

Offsets  along  the  intersection  of  the  hull  surface  and  any 
diagonal  plane 

Offsets  along  any  buttock  line 

The  coefficients  of  a  Thielheimer  equation  for  any  waterline  or 
station  from  the  surface  equation 

The  coefficients  for  standard  cubic  equations  along  any  waterline 
or  station 

Offsets  along  any  waterline  or  station  at  a  variable  interval, 
such  that  a  straight  line  connecting  each  succeeding  pair  of 
points  comes  exactly  at  a  given  tolerance  from  the  actual  curve 

Heights  and  offsets  of  the  intersection  of  the  hull  equation  with 
curves  described  by  equations  developed  in  Step  4,  or  with  any 
cubic,  parabola,  or  straight  line. 


The  data  from  this  program  is  directly  usable  in  the  programs  for  develop¬ 
ing  other  loft  information  which  follow  it. 
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Section  VII 


FRIPARAIION  OF  NUMERICAL  CONTROL  DATA  (SIX?  6) 


In  8tep  5,  one  of  the  output  choices  from  GOBACK  was  the  standard  cubic 
equations  for  each  interval  on  any  waterline  or  station.  With  certain 
alterations,  this  is  the  data  necessary  to  describe  the  hull  contours 
for  the  numerical  control  programs. 


The  necessary  alterations  are: 

The  fairing  process  is  carried  on  entirely  in  feet  and  decimals 
thereof,  while  numerical  control  programs  use  inches.  Therefore, 
the  coefficients  must  be  altered  so  that  the  independent  variable 
can  be  entered  in  Inches,  and  che  dependent  variable  will  be  in 
Inches . 

The  results  of  the  fairing  process  are  nearly  always  scaled  in 
some  amount.  Numerical  control  programs  work  in  full  scale  only. 
Therefore,  the  coefficients  must  be  re-scaled. 

The  origin  of  the  mathematically-lofted  surface  is  often  not  at 
the  point  the  numerical  control  programs  require.  Therefore, 
it  must  be  possible  to  translate  from  one  origin  to  the  other. 


A  program  has  been  written  to  make  these  alterations.  This  program  accepts 
the  output  from  GOBACK  exactly  as  it  is  produced.  In  its  present  form, 
it  produces  data  in  such  a  format  that  it  can  be  inserted  directly  into 
AUTOMAP  programs  for  parts  definition.  Program  details  are  in  Appendix  X. 


*See  Programming  System  for  Numerically  Controlled  Flame  Cutting  of  Shins 
Parts  -  Operating  Manual.  Technical  Report  5.0.0  (Vol.  4  of  Final  Report 
under  Contract  INR  s -4427 ) 


1.0. 0-2 


VII- 1 


Section  VIII 


PLOTTING  THE  LINES  LEANINGS  (STEP  7) 


After  obtaining  the  offsets  from  GOBACK,  It  Is  often  desired  to  make  a 
visual  check  on  the  fairness  of  the  surface.  For  this  purpose,  a 

•jf 

routine  was  written  which  uses  the  1620  computer  and  a  Cal Comp  plotter. 

This  routine  (Appendix  F)  accepts  the  output  of  GOBACK  directly  and 
plots  stations,  waterlines,  or  other  curves  as  desired.  While  this 
program  Is  written  for  a  specific  computer  and  plotter,  the  principles 
of  operation  would  be  the  same  for  nearly  any  other  computer-plotter 
combination. 


California  Computer  Products,  Anaheim,  California 
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Section  IX 


DXTXXXUUI10H  Of  0im  LOR  B AEA  <*Hf  8) 


A.  SHELL  PUH  OXTlLOfMENT 

Using  the  offsets  along  sight  edges  end  vater lines  furnished  by 
G0BACK,  the  shell  piste  development  program  furnishes  the  data  for 
the  flat  plate  patterns  of  the  shell  plate.  This  program  Is 
explained  In  Appendix  6.  The  program  will  do  one  plate  or  an 
entire  strafce  of  plating  at  once.  The  output  can  be  converted  to 
numerical  control  Input,  manually  plotted,  or  plotted  by  the 
computer  using  the  plotting  routine  of  Step  7.  The  data  as  punched 
can  be  entered  directly  Into  the  plotting  program. 

When  using  plating  developed  by  this  program  Intermixed  with  plating 
developed  manually,  allowance  must  be  made  for  a  shift  In  seam  loca¬ 
tion.  This  shift  is  fully  explained  in  Appendix  G  of  this  report. 


1.  STMirant  BSVXLOnSNT 

Stringer  plates  differ  from  shell  plates  in  that  they  have  curvatures 
In  the  longitudinal  direction  only.  That  Is,  plates  whose  cross 
sections  are  straight  lines  but  whose  longitudinal  edges  are  curves. 

These  plates  are  easily  handled  by  the  shell  plate  development  program. 
It  is  necessary  to  provide  a  set  of  data  describing  the  mid line  of 
the  plate  just  as  though  a  she lip late  were  Involved.  In  the  case  of 
the  stringer,  this  information  Is  redundant,  but  it  Is  needed  for 
mechanical  reasons. 
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C.  nu»  AND  PU33  BENDING  OFF  SITS 


A  program  for  computing  the  offsets  for  frame  bending  templates  Is 
presented  In  Appendix  H.  This  program  is  equally  applicable  to 
computing  offsets  of  plate  bending  templates. 

The  program  requires  an  equation  for  the  molded  shape  of  the  frame. 
It  will  furnish  offsets  for  any  portion  of  the  frame  described  by 
the  equation  at  any  desired  interval.  The  program  accepts  the 
equation  in  the  form  of  a  single  standard  cubic  for  —ch  frame 

segment. 


D.  INNER  HULL  OFFSETS 

Inner  Hull  contours  of  naval  combatant  ships,  such  as  the  DIG  class, 
generally  have  the  same  characteristic  as  the  outer  hull  contours. 
That  is,  the  surface,  or  a  large  part  of  the  surface,  must  be  fair. 

That  part  of  the  inner  hull  which  has  a  fairness  requirement  can 
be  faired  by  exactly  the  same  methods  as  are  used  for  the  outer  hull. 
Using  Steps  1  through  8,  the  preliminary  offsets  are  put  through  the 
same  procedure  as  those  for  the  hull  itself. 


E.  INNER  EDGE  OF  NIB  FRAMES 

In  areas  where  the  inner  hull  contours  have  not  been  faired,  many 
web  frames  contain  inner  contours  which  are  related  to  the  outer 
hull  contour.  This  relation  is  often  a  so-called  "linear  taper.” 
The  offset  between  the  two  curves  is  given  at  the  upper  and  lower 
ends  of  the  part,  and  the  intermediate  offsets  are  assumed  to  vary 
in  a  linear  manner. 

A  program  has  been  written  to  produce  this  taper.  It  is  presented 
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in  Appendix  I.  The  program  uses  points  describing  the  outer  hull 
contours  and  the  upper  and  lower  offset.  It  then  finds  points 
describing  the  shape  of  the  inner  hull  contour. 


F.  DECK  CAMBER  OFFSETS 

In  Appendix  J,  a  program  for  defining  decks  is  given.  This 
program  calculates  the  vertical  offsets  from  a  horizontal  plane 
at  a  distance  equal  to  the  molded  depth  above  the  baseline  plane, 
to  the  deck  for  several  configurations.  These  are: 

Straight  line  sheer  and  straight  line  camber 

Straight  line  sheer  and  parabolic  camber 

Parabolic  sheer  and  straight  line  camber 

Parabolic  sheer  and  parabolic  camber 

Offsets  can  be  obtained  at  any  desired  transverse  interval  at  each 
frame. 
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Appendix  A 


SMOOTHING  SOUTINES 


There  ere  two  smoothing  programs.  One  is  for  two-dimensional  curves 
and  the  other  is  for  matrices  of  offsets  describing  a  surface. 

The  three-dimensional  (3-D)  smoothing  routine  is  the  same  as  the  two- 
dimensional  (2-D)  routine  with  an  added  iteration  routine.  Hu  iterating 
routine  smooths  the  points  in  one  direction  (water lines)  making  necessary 
corrections,  then  smooths  the  points  in  the  other  direction  (stations), 
Iterating  between  the  two  directions  until  no  further  changes  are  required 
in  the  points.  Several  iterations  may  be  required,  of  course,  since  a 
change  in  offset  while  smoothing  in  one  direction  may  cause  bad  points  in 
the  other  direction . 

On  any  given  curve,  a  change  in  offset  will  be  made  if,  in  the  middle  of 
a  string  of  second  differences  of  one  sign,  a  single  second  difference  of 
the  opposite  sign  is  found.  If  a  single  second  difference  of  the  opposite 
sign  on  either  end  of  a  curve  is  found,  it  will  not  be  changed.  Since 
the  character  of  the  curve  cannot  be  examined  beyond  this  point,  the  sign 
must  be  assumed  to  be  legitimate,  as  indeed  it  may  well  be.  The  conditions 
just  described  are  illustrated  below: 

O  Only  these  points  will  be  changed 

+ - © - 

++++©++++- 

++++-.._++ 

+  +  +  -  -  +  +  +  ++  - 

+  +  ©+  +  +©  + 

If  the  program  finds  a  case  such  as  follows 

+  +  +  +  +  -  +  -  +  +  +  +  + 

it  will  terminate  execution.  There  are  several  ways  to  change  this  set 
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to  obtain  a  correct  set  of  signs,  and  it  will  be  possible  to  detect  the 
correct  one  only  by  examining  the  lines  drawing. 


3-D  SMOOTHING  ROUT I.® 


This  program  is  written  in  FORTRAN  I  for  an  IBM- 1620  computer 
a.  Input  Data 

Symbols  Used  in  the  FORTRAN  program: 


MA  -  The  number  of  water lines  contained  in  the  data 

NA  •»  The  number  of  stations  contained  in  the  data 

WL(I)  ■  The  vertical  distance  from  the  origin  of 

the  data  to  waterline  i 


STA(J) 

IASF 

IASI 

IASE 
ASET  (I,J) 


The  horizontal  distance  X  from  the  origin  of 
the  data  to  Station  J  ^ 

The  number  of  feet  in  the  offset  (Y,  ,)  when 
entered  in  feet,  inches,  and  eighths^orm 

The  number  of  inches  in  the  offset  when  entered 
in  feet,  inches,  and  eighths  form 

The  number  of  eighth-inch  intervals  in  the  offset 

The  offset  (Y  ,)  when  entered  in  decimal  form 
after  conversioAJfrom  feet,  inches  and  eighths 


Header  Card: 


The  data  on  this  card  consists  of  two  fixed  point  numbers. 

These  numbers  must  be  right  justified  in  the  data  field. 

Variable  Format  Card  Columns 

MA  !  1  -  5 

215 

NA  ^6-10 

The  remaining  data  consists  of  the  locations  of  and  values  for 
the  offsets.  These  offsets  can,  by  an  appropriate  switch  setting, 
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be  read  by  tbe  program  in  either  feet  and  decimals  of  a  foot 
form,  or  in  feet,  inches  and  eighths,  similar  to  the  preliminary 
table  of  offsets.  The  station  and  waterline  coordinates  are 
always  in  feet  and  decimal  form. 

This  data  consists  of  several  subsets  of  data  -  one 
subset  for  each  waterline  in  the  offset  matrix.  The  first 
card  in  each  subset  contains  the  height  of  the  waterline  (Z,j). 
The  remaining  cards  of  the  subset  each  contain  a  station 
location  (X.)  and  the  offset  at  that  point  (Y.  ,).  The  card 
format,  if  the  offset  is  entered  in  feet  and  decimals,  is: 


Tar lab le 

Format 

Card  Columns 

First  oard  of  a  subset: 

WL(I) 

F10.4 

1-10 

Remaining  cards  of  a  subset 

SIA(J) 

2F10.4 

/ 

J  1  -  10 

ASET (I, J) 

S 

|ll  -  20 

If  the  offset  is  entered 

in  feet. 

inches. 

aud  eighths  from  the 

card,  formats  are: 

Tar lab le 

Format 

Card  Columns 

first  card  of  a  subset: 

Wn(I) 

F10.4 

1-10 

Remaining  cards  of  a  subset: 

STA(J) 

F10.4 

1  -  10 

ISA? 

J 

[ll  -  15 

ISAK 

315 

< 

16  -  20 

ISAS 

1 

21  -  25 

w 

b.  putout 

1.  Typewriter  messages 

If  the  program  successfully  completes  the  smoothing  of  the 
points,  the  message,  FAJRIHC  CCKFUCTB,  is  typed. 
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At  any  time  a  bad  point  in  the  data  is  detected,  the  message 

BAD  POINT,  WATERLINE  - _  STATION  - _ 

LIMITS  - _  to  _ ,  NEW  POINT  - _ 

is  typed .  The  limits  denote  the  boundaries  on  the  value  of 
the  bad  offset  for  the  proper  sign  of  second  difference.  The 
new  point  is  the  actual  offset  value  chosen. 

If  the  indeterminate  case  described  earlier  is  found,  the 
following  message  is  typed: 

ADJACENT  BAD  POINTS 

WATERLINE  -  _ STATION  -  _ 


2.  Punched  Cards 

There  is  no  output  from  this  program  if  the  input  is  smooth 
without  changes.  If  one  or  more  offsets  have  been  changed, 
a  complete  set  of  data  is  punched  in  exactly  the  same  format 
as  the  decimal  form  of  the  input  data. 

Sense  Switch  Settings 

SW  1  -  OFF:  Offsets  are  entered  in  feet  and  decimal  form 

SW  1  -  ON  :  Offsets  are  entered  in  feet,  inches,  and  eighths  form 

All  other  switches  are  ignored 


SAMPIE  DATA  FOR  3-0  SMOOTHING 

•§f 

Feet.  Inches,  and  Eighths  Form  (Switch  1  ON) 

£ _ 51  Iflj  iQ j  .'i  dij  Card  ..Co  lumns 

MA7  NA1  ISAF IASI  IASE 


7 

0.0 

25.5 
51.0 

76.5 
102.0 
127.5 


5 

WL(1) 

STA(l) 

STA(2) 


0 

0 

0 

0 

0 


153.0 

1 

0 

178.5 

1 

0 

2.0 

-  WL(2) 

25.5 

-  STA(l) 4 

0 

51.0 

-  STA(2)5 

7 

6 

76.5 

1 

102.0 

8 

10 

127.5 

10 

8 

153.0 

12 

5 

178.5 

13 

8 

4*0 

25.5 

5 

8 

51.0 

7 

10 

76*5 

10 

0 

102.0 

12 

4 

127.5 

14 

7 

153.0 

16 

7 

178.5 

18 

0 

8.0 

25.5 

7 

6 

51.0 

10 

3 

76.5 

13 

3 

102.0 

16 

3 

127.5 

18 

10 

153.0 

21 

0 

178.5 

22 

6 

12.0 

25.5 

8 

5 

51.0 

11 

9 

76.5 

15 

2 

102.0 

18 

4 

127.5 

20 

0 

153.0 

23 

2 

178.5 

24 

9 

0 

0 

0 

0 

0 

0 

0 

3 

6 

3 

0 

0 

1 

7 

0 

0 

3 
2 
6 
1 
2 

0 

7 

5 
0 

4 
0 

6 

1 

3 

0 

1 

7 

7 

2 


8 


Data  for  3-D  smoothing  In  decimal  form  (Switch  1  OFF)  la  exactly  the 
same  as  the  sample  output  shown  on  Pages  A-6,  A- 7 • 


i  - 0.0-  2 
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Tmmim  lxstxmg  vhek  nuxutms  sample  uu 


260000 20000 3K  •*- 


.  "  1  Clear  Memory  Instruction 

ENTER  SUBROUTINES,  PUSH  START 
1620  FORTRAN  SUBR.  AUTO  DIV  9/30/61 
LOAD  DATA 

BAD  POINT,  WATERLINE-  12,0000  '  STATION-  127,5000 
LIMITS-  20.7916  TO  21,5208  NEW  POINT- 


21,1562 


FAIRING  COMPLETE 
STOP 


OUTPUT  OBTAINED  FROM  PROCESSING  SAMPLE  DATA 

7  5 

•0000 

25,5000  1.0000 

51.0000  1.0000 

76.5000  1.0000 

102.0000  1.0000 

127.5000  1.0000 

153,0000  1.0000 

178.5000  1.0000 

2.0000 

25.5000  4.0312 

51,0000  5.5625 

76.5000  7.1145 

102.0000  8.8333 

127.5000  10.6666 

153.0000  12.4270 

178.5000  13.7395 

4.0000 

25.5000  5.6666 

51.0000  7.8333 

76.5000  10.0312 

102.0000  12.3541 

127.5000  14.6458 

153.0000  16.5937 

178.5000  18.0208 
8.0000 

25.5000  7.5000 

5UOOOO  10.3229 

76.5000  13.3020 
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SAMPLE  OUTPUT  FROM  3-D  SHOOTHIMG  IOUTIIB  COWriRUXD: 


102.0000  16.2500 

127.5000  18.8749 
153.0000  21.0000 

178.5000  22.5625 

12.0000 

25.5000  8.4270 

51.0000  11.7812 

76.5000  15.1666 

102.0000  18.3437 

127.5000  21.1562 

153.0000  23.2395 

178.5000  24.7708 


U0»2 


Loor  once  roK  each 

STATION  EACH  ITERATION 
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FORTRAN  LISTING  FOR  3-D  SMOOTHING  ROUTINE 


C  3-D  POINT  FAIRING  PROGRAM  SURFACE 

DIMENSION  ASET( 1 2, 25 ) .WL ( 1 2 ) ,STA( 25 ) .AY (30) 
DIMENSION  RD(25),FSET(30),DEL(25),D(25) 
LAB-1 

READ  1901.MA.NA 
I F (SENSE  SWITCH  1)9200,9110 

9200  DO  9201  1-1. NA 
READ  9300, WL(I) 

DO  9201  J-1 ,MA 

READ  9300,STA(J), ! ASF, IASI , iASE 
ASF- 1  ASF 
AS  I  - 1  AS  I 
ASE-IASE 

9201  ASET(l,J)-ASF+ASI/12,+ASE/96. 

GO  TO  9022' 

9300  F0RMAT(F10.4,3 15) 

9110  DO  9001  1-1,  NA 
READ  1902,WL(I) 

DO  9001  J-1 ,MA 

9001  READ  1902,STA(J) ,ASET( I , J) 

9022  Z-0. 

C  FAIR  EACH  WATERLINE 

ID-1 

LID-1 

N-MA 

9008  DO  3005  J-1,N 

FSET(J)«ASET( ID, J) 

9005  RD( J)-STA(J) 

ADE-WL(SD) 

GO  TO  4000 

9008  DO  9007  J-1 ,N 
9007  ASETC ID, J)-FSET( J) 

ID-ID+1 

IF  (ID-NA)  9008,9008,9009 

9009  LID*>2 

7_n 

LAB-LAB+1 

C  FAIR  EACH  STATION 
JB-1 
N-NA 

9010  DO  9011  1-1, N 
F$ET( I )-ASET( I , JD) 

9011  RD( I )»WL( 1 ) 

ADE-STA(JD) 

GO  TO  4000 

9012  DO  9013  I-1.N 

9013  ASET( I , JD)-FSET( I ) 
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JD-JD+1 

IF  (JO-MA)  9010,9010.9014 
9014  IF  (Z)  9017,9017,9019 
9019  LAB-LAB+1 
GO  TO  9022 

C  NO  CHANGES  FO  A  COMPLETE  CYCLE - FAIRING  COMPLETE 

9017  PRINT  1906 

IF  (LAB-2)  9080,9018,9018 

9018  PUNCH  1901 ,MA,NA 
00  9100  1-1, NA 
PUNCH  1902,WL(I) 

DO  9100  J-1.MA 

9100  PUNCH  1902,STA(J),ASET(I,J) 

9080  STOP 

1901  FORMAT  (1 5, 15. F1 1.6) 

1902  FORMAT  (2F10.4) 

1906  FORMAT  (///16HFAIRING  COMPLETE) 

C  MAIN  PROGRAM 

4000  R-0.' 

4001  BA-O. 

00  4003  1-2, N 

4003  D(l)-RD(I)-RD(I-1/ 

K-N-1 

4004  DO  4005  I-2.K 

DEL( I )-FSET( I )-(FSET(  1+1 )*D(  I  )+FSET( 1-1 )*D( 1+1 ) )/(D( I )+D( 1+1 ) ) 

4005  DEL( I )-SGN(OEL( I )) 

KK-K-1 

DO  4020  l-3,KK 

IF  (ABSF(DEL( I )-DEL( 1-1 )-DEL( 1+1 ) )-2.)  4019,4019,4007 

4007  IF  (BA)  4080,4008,4009 

4008  UPSET-FSET(I)' 

FSET(l)-99876. 

BA-1. 

GO  TO  4020 

4009  SA-(FSET(  1+1  )*D(  1-1  )+FSET(  l-2)*(D(  1+1  )+D(  I )  ) )/(RD(  1+1  )-RD(  1-2)  ) 
SA-SGN(UPSET-SA) 

SC-(FSET( 1+1 )*(D( I )+D(  1-1 ) )+FSET(  l-2)*D( 1+1 ) )/(RD( 1+1 )-RD( 1-2) ) 
SC-SGN(FSET(I)-SC) 

SB-(FSET( l+2)*(D( 1+1  )+D( I ) )+UPSET*0( 1+2) )/(RD( l+2)-RD( 1-1 ) ) 
SB-SGN(FSET(  l+O-SB) 

SD-(FSET( I )*D( I— 2)+FSET( l-3)*(D( I )+D( 1-1 ) ))/(RD( I )-RD( 1-3) ) 
SD-SGN(FSET( I— 2)— SD) 

TA-OELi l-2)+SB-SA 
TC-DEL( 1+1 )+SD-SC 
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R-R+1 . 

IF  (AB$F(TA)-2,)  4014,4014,4015 

4014  FSET(  I-O-UPSET 
DELD-SB 
DEIU+U-SB  ' 

FSEKD-99876. 

BA— 1« 

GO  TO  AO 20 

4015  IF  (ABSF(TC)-2.)  4017,4017,4080 
4017  BA— 1. 

GO  TO  4020 
4080  IDD-l-1 

PRINT  1409 

GO  TO  (Sf01.3102).LID 

3101  PRINT  14O8.A0E.RDCI) 

STOP 

3102  PRINT  1408tRD(l)vADE 
STOP 

1409  FORMAT  C/////19H ADJACENT' BAD  POINTS) 

1408  FORMAT  £ WHWATERLINI-,F8.4,10X,8HSTATI0N-,F9.4) 

4019  BA-O. 

4020  CONTINUE 
DO  4021  I 
IF  (FSETC 

4092  SL-IFSETC 
SA-1F$£TC 
SB«(FSET( 

IF  (ABSFC 

41 19  SA-SB 

4120  FSET(l)«$l+0.5*SA 
XUP-Sl+SA 
XLQW-Sl 

Z-Z+1 , 

60  TO  (3001 .3002) .LID 

3001  PRINT  1003»ADE,RD(I) 

60  TO  30®3 

3002  PRINT  W03.RDOI.ADE 

3003  CONTINUE 

PRINT  1004.XL0M.XUP.FSETO) 

4021  CONTINUE 

60  TO  (9006.9012) .LID  '  „ 

1003  FORMAT  I/21HBAD  POINT-,  WATERLINE-, F8.4.3X,  8HSTATI0N-.F9.4) 

1004  FORMAT  (7HLIMSTS-,F12.4,4H  T0.F12.4,5X,10HNEW  P0INT-.F12.4) 
END 


C 


m ^ 

0-99876.)  4021,4092.4021 
I+0*DCIHFSETCI-|]*D0+11..  .  ,  .  . 

fl-1)*(DC  8>+D(l-1))*“FSET(  0-2)*D(l))/D(  l«  % 

1+1 )*{0( 1+1 )+D( l+2))-FSET( 5+2)*D( l+1))/D( D+2)-SL 
)  4120,4120,4119 


— |)*D(1+1))/CD(t+1)+p( I) j  ^ 
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2-0  SMOOTHING  BOUTINS 


This  program  is  written  in  FORTRAN  I  for  an  IBM-1620  computer.  The 
program  is  exactly  the  same  as  the  3-D  Smoothing,  without  the  iterating 
routine.  Thus  it  is  equipped  to  handle  just  one  curve  at  a  time. 


a.  Input  Data 

Symbols  Used  in  the  FORTRAN  program: 

N  ■  Number  of  data  points  to  define  the  curve 
RD(I)  “  Ordinate  of  the  offset  (X^) 

IFSF  *  The  number  of  feet  in  the  offset  (Y.)  when  entered 
in  feet,  inches,  and  eighths  form 

IFS1  ■  The  number  of  inches  in  the  offset  when  entered  in 
feet,  inches,  and  eighths  form 

IFSB  -  The  number  of  eighths  inches  in  the  offset  when 
entered  in  feet,  inches,  and  eighths  form 

FSKT(I)  ■  The  offset  (Yi>  when  entered  in  decimal  form  or 
after  conversion  from  feet,  inches,  and  eighths 


Header  Card  Format: 

Variable 

Format 

Card  Columns 

N 

15 

1  -  15 

Offset  Data  Card  Format: 


N  cards,  the  set  of  which  must  have  one  of  the  following  formats: 


Variable  Format 


Card  Columns 


Sense  9witch  1  ON: 

ND(I) 

IFSF 

IFSI 

IFSB 

Sense  Switch  1  OFF: 

U>(I) 

FSBT(I) 


F10.4 

315 


2F10.4 


1-10 


-  15 

-  20 
-  25 


-  10 
-  20 
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b.  Output 

1.  If  the  point*  *r«  saooth  without  okaage,  a  halt  with  a 
UUIX  WO  m»  light  will  occur. 

2.  If  any  point  need*  change  it*  location,  tha  Halt*  of  tha 
possible  offset  value*  and  eh*  oeleated  offset  value  will  be 
printed. 

3.  the  indeterminate  case  described  in  tha  3-i  Writeup  will  cause 
the  point  nuabers  of  the  offending  offsets  and  tha  word  ROT 
to  be  printed  on  the  typewriter. 

c.  Sense  ftritche* 

Sane  as  3-D  Smoothing  Trogran 

C 


c 
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TYKWimi  LISTING  WHIN  EXECUTING  THE  2-D  SMOOTHING  PKOGIAM 


260000200003  <<  Instruction  for  Cleerlng  Memory 

LOAD  DATA 

BAD  POINT,  STA*  25,5000 

LIMITS-  20,7916  TO  21.5208  NEW  POINT-  21.1562 
STOP 


SAMPLE  INPUT  DATA 

7 

5.1 

8 

5 

1 

10.2 

11 

9 

3 

15.3 

15 

2 

0 

20.4 

18 

4 

1 

25.5 

20 

0 

7 

30.6 

23 

2 

7 

35.7 

24 

9 

2 
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FLOW  DIAGRAM  FOR  2-D  SMOOTHING  ROUTINE 
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D04O£IX»*KK 


SL=(yUi*  Di  4  Yt-i*  DUt) /(Di  + 1 4  Di) 

SA  RYi-i#(Di4Dl-0-Vi-2*  Dl)/D‘t-i]  -  SL 
SB[(Yi4i*(Di41  +  Dl4l) -Y142*  Di4l)/bl4^-6L 


FORTUM  LI5THG  TOR  2-B  BHWTHDW  RDOTUK 


I 


C  2-0  POINT  FAIRING  PROGRAM 

DIMENSION  RD(20) , FSET(20) ,0(20) , DEL (20) 

READ  1001, N 

! F( SENSE  SWITCH  1)9100,9000 

9100  DO  9101  1-1, N 

HEAD  1005, RD( I), IFSF, IFSI , IFSE 

FSF-IFSF 

FSI-IFSI 

FSE-IFSE 

9101  FSET(l)-FSF+FSI/12.+FSE/96. 

SO  TO  4000 

1005  F0RMAT(F10.4,3 15) 

9000  DO  1  1-1, N 
1  READ  1002,RD( I ) ,FSET( I ) 

4C00  R-0. 

4001  BA-O. 

DO  4003  1-2, N 

4003  D( I )-RD( I )-RD( 1-1 ) 

K-N-1 

4004  DO  4005  1-2. K 

DEL( I )«FSET( I )-(FSET( l+1)*D( I )+FSET( 1-1 )*D( 1+1 ) )/(D( I )+D( 1+1)) 

4005  DEL( I ) — SGN ( DEL ( I ) ) 

KK-K-1 

DO  4020  I-3.KK 

IF  (ABS(DEL ( 1 )-DEL( 1-1 )-DEL( I +1 ) ) —2 . )  4019,4019,4007 

4007  IF  (BA)  4080,4008,4009 

4008  UPSET-FSET(l) 

FSET(  0-99876. 

BA-1. 

GO  TO  4020 

4009  SA-(FSET(I+1)*D(I-1)+FSET(I-2)*(D(I+1)+D(I)))/(RD(I+1)-RD(  1-2)) 
SA-SGN(UPSET-SA) 

SC-(FSET(I+1)*(D( I )+D( 1-1 ) )+FSET( I -2) *D ( I +1 ) ) / ( RD ( I +1 )-RD( 1-2)) 
SC-SGN(FSET( ! )-SC) 

SB-(FSET( l+2)*(D( 1+1 )+D( I ) )+UPSET*D( 1+2) ) / ( RD ( l+2)-RD( 1-1 ) ) 
S3=SGN(FSET( i +1 )-SB) 

SD-(FSET( :}*0(!-2)+FSET(l-3)*(D(l)+D( I —1 ) ) ) / ( RD ( t )-RD( 1-3)) 
SD-SSN(FSET(  I  — 2 >  — SD ) 

TA»DEL( l-2)+SB-SA 
TC-DEL ( I +1 )+SD-SC 
R-fl+1 . 

IF  (A3S(TA)-2.)  4014, 4014, i«015 

4014  FSET( 1-1 )«UPSET 
DEL ( I ) — SU 

OELO+D-SB 
FSET(  0  -99876. 

BA— 1 . 

CO  TO  4020 

4015  IF  (ABS(TC)-2.)  4017,4017,4080 
4017  BA— 1. 
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4080 

1408 

4019 

4020 

4092 

4119 

4120 


4021 

1001 

1002 

1003 

1004 


60  TO  4020 
100-1-1 

PRINT  1408,100,1 
FORMAT  (12,12) 

STOP 

BA-O. 

CONTINUE 
DO  4021  l-3,KK 

IF  (FSET( u-99876.)  4021,4092,4021 

SL-(FSET( 1+1 1*0 ( I )+FSET( l-1)*D( l+1))/(D(l+1)+0(l)) 

SA-1FSET1 1-1 )*(D( I )+0( 1-1 ) )-FSET( l-2)*0( I ))/D( 1-1 )-SL 
SB-1FSET1 1+1 )*(0( 1+1 )+0( 1+2) )-FSET( l+2)*D( 1+1 ) )/D( l+2)-SL 
IF  (ABS(SA)-ABS(SB))  4120,4120,4119 
SA-SB 

FSET(I)-SL+0.5*SA 

XUP-SL+SA 

XLOW-SL 

PRINT  1003,RD(I) 

PRINT  1004,XLOW,XUP,FSET(I) 

CONTINUE 
FORMAT  (15) 

FORMAT  (2F10.4) 

FORMAT  (15HBAD  POINT,  STA-.F10.4) 

FORMAT  (7HUMITS-,F12.4,4M  T0,F12.4,5X,1OHNEW  P0INT-.F12.4) 

STOP 

END 
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Section  I 


INTRODUCTION 


The  purpose  of  the  SMOG  programs  is  to  generate  the  matrices  necessary 
for  input  to  the  linear  programming  routine. 


The  specific  use  of  each  program  presented  in  this  Appendix  is  listed 
below. 


SMOG  1  -  A  one  pass  program  which  generates  a  matrix  for  a 

surface  without  profile  in  the  dual  formulation  and 
>.  fit. 

SMOG  2  -  A  two  pass  program  which  generates  a  matrix  for  a 
surface  with  profile  in  the  dual  formulation  and 
A  fit.  Pass  1  generates  the  deviation  constraints; 
Pass  2  generates  the  curvature  constraints. 

SMOG  3  -  A  two  pass  program  which  generates  a  matrix  for  a 
surface  with  profile  in  the  primal  formulation  and 
sum  of  the  deviations  fit.  Pass  1  generates  the 
deviation  constraints;  Pass  2  generates  the  curvature 
constraints. 

SMOG  A  -  A  single  pass  program  which  generates  a  matrix  for 
fairing  a  two-dimensional  line  with  the  primal 
formulation  and  A  fit,  using  the  1620  linear  program 


There  is  no  requirement  that  the  water  lines  or  stations  used  by  this 
program  be  equally  spaced.  It  is  desirable  because  of  the  variation  in 
sise  of  the  matrix  elements  produced,  however,  to  have  them  as  equally 
spaced  as  possible 
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SECOND  DIFFERENCE  EQUATIONS 


The  equetlon  which  is  used  Co  calculate  the  second  differences  is  given 
below.  The  normalized  values  of  the  coordinates  are  used  in  finding 
these  second  differences. 


For  water lines,  the  equation  is: 


2 _ 

(xJ+l  "  XJ-1? 


xj+i  “  xj 


lLLlZh±± 

*i  '  xJ-i 


For  stations,  the  equation  is: 


SD 


(* 


i+1 


Zi-1> 


zi+i  +  zi 


i-1 


NORMALIZING 


The  procedure  for  normalizing  the  x  and  z  coordinates  is  given 
below.  If  the  program  establishes  the  normalizing  and  translating 
factors  for  the  coordinates,  these  are  found  as  follows: 


S  - 

x 


Tx 


Sz  “  Z1  •  z0  TZ  “  Z0 


Each  x  and  z  coordinate  is  then  altered  as  follows: 


z 


i 


If  it  is  desired  to  scale  and  translate  different  amounts  than  this,  the 
values  Sx  ,  Tx  ,  Sz  ,  and  Tz  may  be  read  in. 
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OUTPUT  DATA 


The  output  data  from  each  of  the  programs  except  SMOG  4  is  in  standard 
SHARE  format.  This  is  the  accepted  linear  programming  format  established 
by  the  IBM  704,  7090  Users  Group.  Examples  of  each  of  the  data  cards 
punched  by  these  programs  is  shown  in  Fig.  B-l.  Figure  B-2  shows  the 
organization  of  the  data  deck  when  ready  for  use  with  the  linear 
program. 
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Fig.  B-l  LP/90  Card  Formats 

(SHARE  Standard  Format) 
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Baa la  (optional) 


Card  Order 


Section  II 


SMOG  1  AMD  2 


The  purpose  of  the  programs,  SMOG  1  end  SMOG  2,  is  to  generate  the 
matrices  required  for  surfaces  with  the  following  characteristics: 

(1)  Dual  formulation 

(2)  X  (minimax  fit) 

(3)  Double  splined  in  x  direction 

(4)  Single  or  double  splined  in  the  z  direction 

In  addition,  SMOG  2  produces  a  matrix  with  the  modifications  necessary 
for  end  profiles.  The  organization  of  each  program  is  given  in  Fig.  B-3. 


The  surface  equation 


Figure  B-4  shows  a  five-station  by  five -water line  surface  and  the 
equation  representing  this  surface.  When  producing  a  matrix  to  describe 
this  problem  to  the  linear  program  we  must  generate  a  requirement  that 
the  surface  come  close  to  each  offset  on  the  grid.  The  requirement 
is  of  the  form 


(•) 

(b) 


X-  Ytx^Zj) 

X+  Ytx^Zj) 


i  -  1 

i  -  i 


...  ,  N 
...  ,  M 


/ 


One  (a)  and  one  (b)  column  (constraints  are  columns  in  the  dual)  are 
generated  for  each  offset.  For  instance,  offset  (1,3)  Fig.  B-4,  would 
have  the  requirements 


(a)  >  -  AQ1x3  -  AQ2x3  -  4Q3x3  -  _y1>3 

ib)  X+  AQjX3  +  *q2x3  +  *03*3  “  Yif3 
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In  Che  Surface 
Program,  Smog  2, 
Pass  I  and  Pass  II 
are  separated  here 


Pig.  B-3  Organisation  of  SMOG  1  and  SMOG  2 


Y  "  A00  +  A01X  +  A02*2  +  A03x3  + 

+  (AHX  +  A12*2  +  A13*3  +  A14(X-X3>#)  * 

+  +  A^x2  +  Ajj^x3  +  A24(x-x3)3|  z2 

+  |a31x  +  A32x2  +  A33x3  +  A34(x-x3)^|  z3 

+  (A41x  +  A42*2  +  A43*3  +  A44(x-x3>  J 

*  ^°^nts  where  the  sign  of  the  x  curvature  is  enforced 
\  Points  where  the  sign  of  the  z  curvature  is  enforced 
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The  constant  coefficient  A  Is  given  the  value  of  the  offset 

oo 

(x^,s.p  and  Is  not  Included  In  the  matrix.  (The  surface  always  goes 
through  this  offset,  therefore  A  Is  known). 

OO 


Similarly,  the  constraints  for  Point  (3,3)  are: 


(a)  X 


-  [Aoix: 

-[ 


3  +  A02x3  +  A03x33 


A11x3  +  A12x32  +  A13X33 


-  j^jX.j  +  A22x3  +  A^Xj 

‘  [A31X3  +  A32X3  +  A33X3 


-  -Y(3,3) 


[ 


(b)  X *  |A0lx3  +  Aq2x3  +  Aq3x3 


+  Jaux3  +  A12x32  +  A13^3 

+  +  ^22*3  + 


A31x 


3  +  A32x3  +  A33x3 


-  Y(3,3) 


There  are 

offsets. 

A 

ri,J  *re 
Appendix 


also  constraints  on  the  sign  of  the  curvature  at  certain  of  the 
These  constraints  have  the  following  form  where  'r*.  and 
the  signs  of  the  second  differences  found  in  Section  1  of  this 


(c) 

"i.J  * 

Y" 

i.  o 

(d) 

Ki  * 

Y" 

i.  o 
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At  Offset  (1,3)  there  is  e  curvature  constraint  in  the  x  direction 
only,  as  shown  in  Fig.  B-4  . 

(c)  *t>.  (AQ2  -  3A03x3)  >  0 

There  are  x  and  z  constraints  both  at  offset  (3,3)  as  follows: 

ri.J  [‘A02  ’  3A03X3  ‘  (A02  +  3A03x3)  *3 

-  (Aqj  +  3aq3x3^  zx  “  ^A02  +  ^03*3^  *3  J  —  0 

ri,J  ["  ^1*3  +  ^2*3  +  ^3*3  * 

-  (A21X3  +  A22x32  +  *23*3  >  3zs]  *  0 

PROVIDING  A  BASIC  FEASIBLE  SOLUTION 

It  is  often  desirable,  when  fairing  a  surface  in  the  dual,  to  provide  the 
linear  program  with  a  basic  feasible  solution.  This  consists  of  an 
identity  matrix  of  slacks  as  shown  in  Fig.  B-5. 

This  identity  matrix  is  produced  by  the  program  given  in  Appendix  K.  The 
program  is  furnished  with  the  number  of  the  last  non-basis  column  and  the 
number  of  rows  in  the  non-basis  matrix.  It  will  then  produce  a  set  of 
basis  headings  and  the  identity  matrix. 

The  set  of  basis  headings  are  inserted  in  the  dual  matrix  immediately 
following  the  row  identifications.  The  identity  matrix  is  placed  just 
ahead  of  the  first  B  card.  This  matrix  can  be  produced  by  SMOG  1,  and 
the  actual  listing  of  the  matrix  for  a  surface  is  included  in  Section  III 
of  this  Appendix. 


in 

i 

PQ 

00 

•H 

Pm 
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Fig.  B-5 


SURFACE  EQUATION  WITH  END  PS0FIUS 


In  order  Co  allow  for  end  profiles  on  the  surfaces,  the  surface  equation 
y  ■  f(x,z)  has  been  modified  by  multiplying  it  by  a  function  T  so 
that  y  ■  f(x,z)  *  T(x,z) 


where 

(e)  T(x,z) 


G(z)  is  the  two-dimensional  Theilheimer  equation  of  the  profile  mentioned 
in  Section  Il-B  of  the  report.  The  equation  G(x)  must  have  the  same 
origin  and  be  scaled  the  same  as  the  surface  equation.  D  is  the 
fraction  of  the  station  spacing  over  which  the  function  is  effective,  end 
P  determines  if  the  slope  of  the  surface  at  the  bow  is  finite  or  infinite. 
The  usual  value  for  P  is  1  for  finite  slope,  1/3  for  infinite  slope. 
The  value  of  r  is  normally  3  . 


For  any  given  value  of  x  and  z  ,  when  generating  a  matrix  for  a 
given  surface  T(x,z)  reduces  to  a  number.  Therefore,  the  deviation 
constraints  can  be  modified  in  the  following  manner. 


(a) 


X-  Y(x^  ,  Zj) 


ii. 


(b)  X+  Y(x.,z  )  -  Yl,j 

J  T^.Zj) 


So  the  only  modification  necessary  for  producing  deviation  constraints  is 
to  divide  the  offset  by  the  value  of  T  .  T  is  given  the  value  one  if 
the  actual  value  is  less  than  a  given  indicator. 
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The  aeeond  derivative  constraint*  modifications  are  more  involved.  The 
constraints  are: 


(1)  for  the  x  direction 


’  d2 

2 

_dx 

/ 

(2)  for  the  y  direction 


[-4  ■  T(xi«v) 


£  0 


Using  the  rule  for  differentiating  the  product  of  two  terms  twice,  these 
become : 

(!)  t1#J  pOii.'j)  ‘  T,,(*i»*j>xx  +  2(T’(xi'*j)x  *  f’<Xi'Xj>x) 

+  T(xi,aj)  *  0 

(2)  *  T”(x1,rj)K  +  2(T'(x1,ej)z  *  f'G^),) 

+  T(xi,zj)  •  fH(xi,zj)MJi  0 

These  constraints  are  produced  for  the  linear  programming  matrix  by  SMOG  2. 
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A.  SMOG  1 


OPERATING  INSTRUCTIONS 

This  program  produces  a  surface  matrix  without  profile  conditions  in 
the  dual  formulation  and  Xfit.  The  surface  will  be  double  splined 
in  the  x  direction  and  may  be  either  single  or  double  splined  in 
the  z  direction.  The  program  is  written  in  FORTRAN  II  for  an 
IBM- 1620  computer. 


Fortran  Symbol  Definitions 


M  -  Total  number  of  waterlines  in  the  surface  to  be 

faired.  This  will  be  an  odd  number  if  the  surface 
is  double  splined  in  the  z  direction. 

N  -  Total  number  of  stations  in  the  surface.  This  will 

always  be  an  odd  number. 


X(J)  -  Distance  along  the  x  axis  from  the  origin  of  the 

surface  to  Station  J  . 

Z(I)  -  Distance  along  the  z  axis  from  the  origin  of  the 

surface  to  waterline  i  . 


Y(I, J)  -  Offset  of  the  surface  at  waterline  i  ,  station  j  , 
in  feet  and  decimals 


IAY 

IBY 

ICY 

NAA 


SX  . 

TX 


Number  of  feet  in  the  offset  if  read  in  feet,  inches 
and  eighths. 

Number  of  inches  corresponding  to  IAY 
Number  of  eighths  corresponding  to  IAY 
Program  indicator 

If  NAA  ■  -1  scale  factors  for  normalizing  and 
translating  X(J)  and  Z(I)  are  read  in. 

If  NAA  •  +1  scaling  and  translating  are  done 

by  the  program  as  given  in  Section  I  of  this  Appendix. 

Scale  factor  for  X(J)  when  these  are  read  in 

Translation  factor  for  X(J)  when  these  are  read  in 


SZ 


Scale  factor  for  Z(I)  when  these  are  read  in 
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TZ 

RHS 


CURV 


NAB 


Translation  factor  for  Z(I)  when  these  are  read  in 

Right  hand  side  -  In  the  dual  this  is  really  the  cost 
on  the  equation  coefficients.  The  ideal  value  for 
RHS  would  be  zero;  however,  this  makes  the  solution 
susceptible  to  cycling  so  that  a  value  of  .0001  is 
usually  used. 

Value  of  the  requirement  vector  for  the  curvature 
constraints,  this  value  actually  goes  in  the  cost 
row  in  the  dual  and  is  always  made  zero. 

Program  indicator 

If  NAB  "  -1  ,  the  surface  will  be  double  splined 
in  the  z  direction  and  M  must  be  an  odd  number. 

If  NAB  "  +1  ,  the  surface  will  be  single  splined 
in  the  z  direction  and  M  is  unrestricted. 


The  input  data  cards  are  described  below.  There  are  certain  input 
variables  on  which  there  are  restrictions  on  the  value  which  can  be 
entered.  These  restrictions  follow. 

Variable  Minimum  No.  Maximum  No. 

N  3  15 

M  3  11 


In  describing  the  various  data  cards,  the  actual  FORTRAN  format  field 
is  used  in  most  cases.  These  fields  come  consecutively  across  the 
card  with  no  gaps  or  blank  columns  between,  except  where  indicated. 

The  field  descriptions  are  the  FORTRAN  F  field,  which  uses  the  FORTRAN 
fixed  point  decimal  number,  and  the  I  field  that  uses  the  FORTRAN 
integer  number,  which  is  always  right  justified.  The  card  numbers 
are  not  punched  on  the  data  cards . 
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Content*  of  Card 


Card  Mo 


Card  Columns 

1 

2-12  13-53 

Contents 

* 

Blank  Any  alphameric  problem 

heading  (This  heading  will 
be  used  for  Identifying  the 
problem  on  the  LP  output) 

Format 

15 

15 

Variable 

NAA 

NAB 

Format 

15 

15 

F11.6 

F11.6 

Variable 

M 

N 

CURV 

RHS 

Format 

F10.4 

Variable 

Z(I) 

Alternate  Cards : 

— 

Format 

F10.4 

F10.4 

Variable 

X(J) 

Y(I.J) 

> 

or: 

Format 

F10.4 

15 

15 

15 

Variable 

XU) 

IAY 

IBY 

ICY 

This  arrangement  can  be  stated  as  follows: 


2 

3 


N  Cards 
per  Group 


Next 

M 

Groups 


The  first  card  of  each  set  gives  the  height  z  of  the  waterline.  The 
following  cards  each  give  a  station  coordinate  and  the  offset  at  the 
intersection  of  the  given  waterline  and  station.  There  is  one  Set 
of  cards  for  each  waterline  in  the  surface. 


Format  F10.4  F10.4  F10.4  F10.4 

Variable  SX  SZ  TX  TZ 

(This  card  only  used  if  NAA  ■  -1 


Last  card 


Output 

The  output  from  this  program  is  designed  to  be  read  directly  by  the 
LP-90  program  after  card  to  magnetic  tape  conversion. 


The  first  card  punched  contains  the  Heading  Read  in  the  first  input 
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card.  The  next  set  of  cards  are  the  row  Identification  cards. 

These  identify  each  row  of  the  matrix  to  be  used  in  this  solution. 

For  this  matrix  all  of  the  rows  ere  always  used  end  identified, 
including  the  cost  row  which  has  the  name  "OFSET." 

After  these  comes  a  card  punched  "MATRIX,"  end  then  the  input  matrix 
to  the  linear  program.  The  matrix  is  in  column  form  with  the  rows  in 
order  within  the  column.  Zero  values  are  not  punched. 


The  card  formats  are  given  in  Fig.  B-  1 .  The  order  of  the  output 
data  is  shown  graphically  in  Fig.  B-  2  .  The  sample  problem  output 
shows  exactly  what  the  output  data  consists  of. 


Sense  Switch  Settings 


Sense  Switch  1 
Sense  Switch  2 

Sense  Switch  2 
Sense  Switch  3 


Not  used 
ON 


OFF 

ON 


Sense  Switch  3  -  OFF 


The  values  of  the  x  and  s  second 
differences  and  the  resulting, signs  are 
printed 

The  above  is  not  printed 

The  offsets  are  read  in  as  feet,  inches 
and  eighths 

The  offsets  are  read  in  as  feet  and 
decimals 


SAMPLE  PROBLEM 

The  sample  problem  is  a  five-station  by  five-waterline  set  of  offsets 
taken  from  the  AG-88  class  icebreaker.  The  stations  used  are  Stations 
13,  30,  32,  34,  and  36  foot  levels.  The  matrix  is  for  an  equation  double 

splined  in  both  the  x  and  t  directions.  .  ■ 

*  if- nj 
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SAM?  12  PROBUM  FOR  SMOG  1 


INPUT 


*  5X5  Z  DOUBLE  SPUNED  AG-88  ICE  BREAKER 


1  -1 

5  5  .OOOOOO  .0001 


28.0000 

137.5000 

30 

3 

0 

150.0000 

29 

2 

2 

162.5000 

27 

6 

4 

175.0000 

25 

2 

4 

187.5 

21 

10 

6 

30.0000 

137.5000 

30 

7 

5 

150.0000 

29 

7 

7 

162.5000 

28 

2 

5 

175.0000 

26 

2 

4 

187.5000 

23 

3 

5 

32.0000 

137.5000 

30 

8 

5 

150.0000 

29 

10 

4 

162.5000 

28 

7 

5 

175.0000 

26 

10 

4 

187.5000 

24 

4 

3 

34.0000 

137.5000 

30 

8 

3 

150.0000 

29 

11 

2 

162.5000 

28 

10 

5 

175.0000 

27 

3 

5 

187.5000 

25 

1 

7 

36.0000 

137.5000 

30 

5 

5 

150.0000 

29 

10 

4 

162.5000 

28 

11 

5 

175.0000 

27 

6 

6 

187.500C 

25 

7 

3 

TYPEWRITER  LISTING  OF  SECOND  DIFFERENCES 
(3aac«  Switch  2  ON) 

5X5  Y.  DOUBLE  SPUNED  AG-88 
2ND  DIFF  -  WL'NES 


Switch  2  OUTPUT 
ICE  BREAKER 


WL  STA  ACT  VALUE  SIGN 


1 

2 

-.003732 

1.00 

1 

3 

-.004400 

1.00 

1 

4 

-.006266 

1.00 

1 

5 

-.006266 

1.00 

2 

2 

-.002933 

1.00 

2 

3 

-.003666 

1.00 

2 

4 

-.005733 

1.00 

2 

5 

-.005733 

1.00 

3 

2 

-.002532 

1.00 

3 

3 

-.003334 

1.00 

3 

4 

-.004799 

1.00 

3 

5 

-.004799 

1.00 

4 

2 

-.001865 

1.00 

4 

3 

-.003400 

1.00 

4 

4 

-.003599 

1.00 

4 

5 

-.003599 

1.00 

5 

2 

-.001999 

1.00 

5 

3 

-.003200 

1.00 

5 

4 

-.003467 

1.00 

5 

5 

-.003467 

1.00 

2ND 

DIFF-5 TAT  IONS 

STA 

WL 

ACT  VALUE 

SIGN 

1 

2 

-.075525 

1.00 

1 

3 

-.026025 

1.00 

1 

/. 

-.052100 

1.00 

1 

*  > 

-.052100 

1.00 

2 

2 

-.062525 

1.00 

2 

3 

-.039050 

1.00 

2 

4 

-.031250 
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I 


IKO.t*0  +  » 


/  rOMCH  \ 
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JTF.TF-I 
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SUBROUTINES 


Subroutine  Pinch 
Offset  Col. No.  Constraint 

.  (YD)  (JFK)  (KFJD) 


Examine  the  sign  of  the  offset 
and  determine  whether  this  is 
constraint  (a)  or  (b),  then 
punch  appropriate  cost  and 
constant  format 


Return 


FORTRAN  LISTING  FOR  SMOG  1 


C 

C 


SMOG1  DUAL  X  OR  X  AND  Z  DBL  SPL  NO  PROFILE 
DIMENSION  RX( 1 1 ,15) *RZ( 15*11) . TRM(150),  IRW(250) 
DIMENSION  X(15),Z(11),Y(11,15),WL( 15,11) ,  IWL(15) 
★  *  I NPUT  *  * 

750  READ  173 

READ  1 72,NAA,NA8 

READ  150,M,N,CURV,RHS 

NL-N+1 

ML-M+1 

NS-N-1 

MS-M-1 

DO  2  1-1, M 

READ  154, Z(l) 

DO  2  J-1 ,N 

I F(SENSE  SWITCH  3)802,803 

802  READ  183,X(J), IAY, IBY, ICY 
AAY-IAY 
BBY-IBY 


CCY-ICY 

Y(l,J)-AAY+(BBY+(CCY/8.))/12. 

GO  TO  2 

803  READ  154,X(J),Y(I,J) 

2  CONTINUE 
IF(NAA)3,3,4 

3  READ  154,SX,SZ,TX,TZ 
GO  TO  5 

4  SX«X(2)-X(1) 

SZ-Z(2 )-Z( 1 ) 

TZ-Z(l! 

TX-X(I) 

*  *  WATERLINE  2ND  DIFFERENCES  *  * 

5  I F ( SENSE  SWITCH  2)6,7 

6  PRINT  173 
PRINT  157 

7  DO  15  1-1, M 


DO  13  J-2,NS 
A-2*/ ( (X( J+1  )-X(  J-1 ) )  ) 
ZA-Y( I ,  J+1 )-Y( I , J) 
A«A*(ZA/(X(J+1)-X(J))- 
.  IF(A)10,9,8 
8  RX(  I ,  J)— 1. 

GO  TO  11 


J+1 )-X( J) )-(Y( I , J)-Y( I , J-1 ) )/(X( J)-X( J-1 ) )) 


9  RX(I.J)-0, 
GO  TO  11 


10  RX( I , J)«1 • 
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( 


( 


c 


c 


c 


DIFFERENCES 


11  I F(SENSE  SWITCH  2)12,13 

12  PRINT  158,I,J,A,RX(I,J) 

13  CONTINUE 
RX(I,N)«RX(I,NS) 

I F ( SENSE  SWITCH  2)14,15 

14  PRINT  158, I , J,A,RX( I ,N) 

15  CONTINUE 

*  *  STATION  2ND 

I FCSENSE  SWITCH  2)16,17 

16  PRINT  159 

17  DO  25  J-1,N 
DO  23  1-2, MS 
A-2./(Z(l+1)-Z(l-1)) 

ZAA-Y( l+1,J)-Y( I ,J) 

A-A*(ZAA/|z( 1+1 )-Z( I ))-(Y( I , J)-Y( l-1,J))/(Z(l )-Z( 1-1 )) ) 
I F( A) 20, 19, 18 

18  RZ(J,I)— 1, 

GO  TO  21 

19  RZ(J, I )-0, 

GO  TO  21 

20  RZ(J, I )-1 • 

21  I F(SENSE  SWITCH  2)22,23 

22  PRINT  158, J, I ,A,RZ( J, I ) 

23  CONTINUE 
RZ(J,M)-RZ(J,MS) 

I FCSENSE  SWITCH  2)24,25 

24  PRINT  158, J, I ,A,RZ(J,M) 

25  CONTINUE 

★  *  NORMALIZING  *  * 

DO  26  J-1 ,N 

26  X(J)-(X(J)-TX)/SX 
DO  27  1-1, M 

27  Z( I )-(Z( I )-TZ)/SZ 

*  ★  ROW  ID  *  * 


PUNCH  173 

NJUL-((N-3)/2)+4 

NhWN 

I F ( NAB ) 719*718,718 
719  JUL-3 
JUL2-2 

NOROW-NJUL*( (M-3)/2+4) 
GO  TO  720 
718  JUL-2 
JUL2-1 

NOROW- (M+2)*NJUL 
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720  PUNCH  160 

DO  28  NR- 1, NOROW 
28  PUNCH  162, NR 

C  *  *  MATRIX 

PUNCH  164 

C  *  *  FIRST 


700 

702 


31 

32 

133 

33 


34 

35 


703 

705 

36 

37 


706 


MLG-(2*NM)-2 

1-2 

DO  37  J-1 ,MLG,2 
Y(1,I)«Y(1,I)-Y(1,1) 

KFJ-1 

CALL  PINCH  (Y(1,I),J,KFJ) 

I  ROW-2 

DO  31  NPR-1.3 
TRM(IROW)-X(l)**NPR 
PUNCH  167, J, IROW,TRM( I ROW) 
IROW-IROW+1 
IF( 1-3)35,35,32 
DO  34  IT-3,1 
IF(I-IT)34,34,133 
I F ( IT-2*( IT/2) |33|34j^3 


TRM( IROW)-(X( I )-X( I T) )**3 
PUNCH  167, J, IROW,TRM( (ROW) 
IROW-IROW+1 
CONTINUE 
JPLS-J+1 
I  ROW- 1  ROW-1 
I WL ( I ) - 1  ROW 
KFJ-2 

CALL  PINCH  ( Y( 1 , I ) , JPLS,KFJ) 
DO  36  IR-2, 1  ROW 
WL( I , IR)-TRM( IR) 

PUNCH  169,JPLS,IR,TRM(IR) 
1-1+1 


*  REMAINING  WL 


DO  64  I Z-2,M 
JL-J 

MG«JL+(2*NM)-2 

IX-1 

DO  63  J-JL,MG,2 
I F— 2 

I RW ( I F)«IR0W+1 

Y ( I Z , IX) — Y ( I Z , I X)-Y(1 , 1 ) 

KFJ-1 

CALL  PINCH  (Y( IZ, I X) , J,KFJ) 


* 


* 


*  * 

WATERL I NE  * 
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708 

41 


42 

142 

729 

730 
144 


145 

146 
709 

711 

43 

44 


244 


45 

46 


IF(J-JL)41,41,44 

DO  42  IZPR-1.3 

TRM(IF)-Z(IZ)**IZPR 

PUNCH  167,J,IRW(IF),TRH(IF) 

IF-IF+1 

IRW(IF)-1RW(IF-1)+NJUL 
I F ( I Z-JUL) 1 46 , 1 46 , 1 42 
DO  145  ITZ-JULJZ 
IF(IZ-ITZ)146.146J29 
GO  TO  (144, 144, 730), JUL 
IF(ITZ-2*(ITZ/2))1M».145,144 
ZTR-(Z(IZ)-Z(ITZ))**3 
TRM(IF)-ZTR 

PUNCH  167,J, IRW( IF),TRM( IF) 
IF-IF+1 

I RW( I F)— I RW( I F-1 )+NJUL 

CONTINUE 

JPLS-J+1 

KFJ-2 

CALL  PINCH  (Y(IZ,IX),JPLS,KFJ) 

JF-IF-1 

DO  43  IF-2.JF 

PUNCH  169,JPLS,IRW(IF),TRM(IF) 
GO  TO  63 

*  *  Z  ,  Z  SQR  ,  Z  CUBE 

IRO  -1 

K-IWL(IX) 

DO  244  1 J-2,K 
IRO  -IRO  +1 

PUNCH  167,J,IR0,WL(IX,U) 

DO  50  IZPR-1,3 
ISTRT-IRW( IFJ+NJUL 
ZTRH-Z( I Z)** I ZPR 
TRH( IF)-ZTRM 

PUNCH  167,J, IRW(  IF) ,TRM( IF) 
IF-IF+1 

IRW(IF)-IRW(IF-1)+1 

DO  45  NPR-1,3 

TRM( I F)-(X( I X)**NPR)*ZTRM 

PUNCH  167, J» IRW( IF) ,TRM( IF) 

IF-IF+1 

IRW(IF)-IRW(IF-1)+1 
IF( IX-3)50,50,46 


DO  49  IT-3. IX 
IF(IX-IT)49,49,47 


ROWS  *  * 
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u 


49 

50 


731 

55 

56 


256 

n 


59 

60 
61 
712 

71% 


262 


62 

63 


IF(IT-2*(IT/2))4at49,48 
TRM(IF)-<(X<IX)-X(ITJ)**3)*ZTRM 
PUNCH  167,J,1RW(IF),TRH(IF) 
IF-IF+1 

IRW(IF)-IRW(IF-1)-H 

CONTINUE 

IRW(IF)-ISTRT 

*  *  THE  I LKE I  HER  ROWS 

IF( IZ-JUL)61 ,61 ,54 
00  60  ITZ-JUL.IZ 
SO  TO  (55,55,731),JUL 
IF( ITZ-2*i  I TZ/2)  >55.60,55 
ISTRT-IRW(IF)+NJUL 
IF(IZ-ITZ)61, 61.56 
ZTRIi»(Z(  IZ)-Z(  ITZ)  )**3 
TRH( IFj-ZTRH  ,  % 

PUNCH  167.J, IRW( IF),TRM( IF) 
IF-IF+1 

IRW( IF)-IRW( IF-1 )+1 
DO  256  NPR-1,3 
trm<if)-<x<i*)**npr)*ztrm 
PUNCH  167, J.  IRW( IF) ,TRM( IF) 
IF-IF4-1 

IRW(IF)-IRW(IF-1H1 

DO  59  IT-3. IX 

IF( IX-IT)59.59.57 

IF( IT-2*l IT/2) j58.59f58 

TRM(IF)— ((X( IX;-Xv IT)5**3)*ZTRH 

PUNCH  167  ,J,IRW( IF).TRH(IF) 

IF-IF4-1 

IRW(IF)-IRW(IF-1)+1 

CONTINUE 

!RW(IF)-ISTRT 

JPLS-J+1 

KFJ-2 

CALL  PINCH  (Y(IZ,IX),JPLS,KFJ) 

IRO  -1 

K-IWL(IX) 

DO  262  IJ-2.K 

IRO  -IRO  4-1  .  .  .. 

PUNCH  169.JFLS.IR0  .WL(IX.IJ) 
JF-IF-1 
DO  62  IS— 2  JF 

PUNCH  169.JPLS. IRW{ IS) ,TRM( IS) 
IX- 1X4-1 


* 
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64  CONTINUE 

C  *  *  X  SEC  DIFFERENCES  *  * 

JX-NS/2 
DO  88  I2:-I,M 
I  “3 

JS-JPLS*1 
JF-JS+.iX-1 
DO  87  J-JS,JF 
PUNCH  1 7 I , J , CURV 
I  ROW-3 
I STRT-3 
WRH-RX(IZ.I) 

501  CALL  SIGN( J, IROW(WRM) 

503  I  ROW-  ’iROW+1 

WRH-?  #*X( I )*RX( I Z,  I ) 

504  CALL  SIGN(J. IR0W,WRM) 

506  I F( '-3  )  75.75,72 

72  i ROW-IROW+1 
DO  74  K-3  *1*2 
l*(K-l)73,75.75  # 

7i  V/RM-3«*(Xt  0— X(K))*RX(  IZ,  I ) 

5T7  CALL  SIGN(J,IR0W,WRM) 

74  IFvOW-IROW+1 

75  IF(IZ-1)87.87.76 

76  DO  79  NPWR-1 ,3 

I S  TRT- 1 STRT+NJUL 
IROW-ISTRT 

WRH«(Z( IZ)**NPWR)*RX( IZ, I) 

509  CALL  SIGN(J,IROW,MRM) 

511  I  ROW- 1 R0W*1 

WRM3.*X( I )*(Z( IZ)**NPWR))*RX( IZ, I ) 

512  CALL  SIGN(J»IROW,WRM) 

514  IROW-IROW+1 
DO  78  K-3, 1.2 
IF(K-I )77|7?*79 

77  WRIV3.*(X( I )-X(K))*(Z( IZ)**NPWR)*RX( IZ, I ) 

515  CALL  SIGN(J,IROW,WRM) 

78  IROW-IROW+1 

S  CONTINUE 

I  F( I Z-JUL) 87,87,81 
81  DO  86  L-JUL. IZ.JUL2 
„  IF(L-IZ)82,87,87 
82  ISTRT-ISTRT+NJUL 
IROW-ISTRT 

WRM-RX(IZ,I)*((Z(IZ)-Z(L))**3) 


l.C.0-2 
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517  CALL  SIGN(J,IROW,WRN) 

519  IROW-IROW+1  ,  , 

WRM3.*X(  I  )*(  (Z(  I  Z)-Z(L)  )**3)*RX(  I Z,  I ) 

520  CALL  SIGN(J,IR0W,WRM) 

522  IR0VMR0W+1 
00  85  K-3» I ,2 

I F(K-I >84,86.86 

84  WR^(3.*(i(l)-X(K))*((Z(IZ)-Z(L))**3))*RX(IZ,l) 

523  CALL  SIGN(J,IR0W,WRM) 

85  I  ROW- 1  ROM-1 

86  CONTINUE 

87  1-1+2 

88  JPLS-JF 

Z  SEC  DIFFERENCES  *  * 

J-J-1 

DO  111  IZ-JUL,M,JUL2 
DO  110  1-1, N 
J-J+1 


PUNCH  171 ,J,CURV 

735  ISTRT-N+6 

736  WRM-RZ(UZ) 

CALL  SIGN(J.ISTRT,WRH) 

562  I F( 1-1 >96,96,91 

91  IROW-ISTRT+1 

Srh2(x(i)J*Wr)*rz(i,iz) 

525  CALL  SIGN(J,IROW,MRM) 

92  IROW-IROW+1 

I F( 1-3)96,96.93 

93  DO  95  K-3.I.S 

94  WRPm! X( I  5-x|k5)**3)*RZ(  I ,  IZ) 

527  CALL  SIGN(J,IROW,WRM) 

529  IROW-IROW+1 

95  CONTINUE 

96  I STRT- I STRT+NJUL 
WRM3  »*Z(  I  Z)*RZ(  I ,  I Z) 

CALL  SIGN(J,ISTRT,WRM) 

532  IROMISTRT+1 

IF( 1-1 >102,102,97 

97  DO  98  NPWR-1,3 

WRH-(X( I )**NPWR)*Z( I Z)*3.*RZ( I , IZ) 

533  CALL  SIGN(J,IROW,WRM> 

98  IROW-IROW+1 

I F( 1-3) 102, 102,99 


99  00  101  K-3,l,2 

IF(K- 1)100, 101. 101 

100  WR^((X(  l)-X(K))**3)*Z(IZ)*3.*RZ(l,IZ) 

535  CALL  SIGN(J, !R0W,WRM) 

537  IR0W-IR0W+1 

101  CONTINUE 

102  I F( I Z— JUL) 110, 110, 103 

103  DO  210  L-JUL, IZ,JUL2 
I F ( L— I Z ) 104,210,210 

104  ISTRT-ISTRT+NJUL 

WRKMZ(  I  Z)-Z(L)  )*RZ(  1 .  1 Z)*3. 

CALL  SIGN(J, ISTRT.WRM) 

540  IROW-ISTRT+1 

IF( 1-1)210,210,105 

105  DO  106  NPWR-1,3 

WRM-(Z(IZ)-Z(L))*(X(  l)**NPWR)*3.*RZ(l,IZ) 

541  CALL  SIGN(J,lR0W,WRM) 

106  I  ROW- 1 R0W+1 

I F( 1-3)210,210,107 

107  DO  109  K-3.l,2 

I F(K— I ) 108, 210. 210 

108  WRH-(Z(  IZ)-Z(L)  )*((X(  I  )-X(K))**3)*3.*RZ(  I ,  IZ) 
543  CALL  SIGN(J,IR0W,WRM) 

109  IR0W-IR0W+1 
210  CONTINUE 

110  CONTINUE 

111  CONTINUE 
PUNCH  174 
PUNCH  175 

738  DO  412  JORG-2, NOROW 
412  PUNCH  1 77, JORG.RHS 
737  PUNCH  176 
PRINT  178 
PAUSE  751 
GO  TO  750 

150  F0RMAT(2I5,F11.6,F11.6) 

154  FORMAT (4F10 ,4) 

157  FORMAT ( 10X , 1 7H2ND  DIFF  -  WLINES// 

C30H  WL  STA  ACT  VALUE  SIGN/) 

158  F0RMAT(2I5,F12.6,F8.2) 

159  F0RMAT(//10X,1 7H2ND  D I  FF— STATIONS// 

C30H  STA  WL  ACT  VALUE  SIGN/) 

160  FORMAT (6HR0W  ID/12X,6H  OFSET) 

162  F0RMAT(12X,2H  R,l4) 


( 
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164  FORMAT (6HMATR IX) 

167  FORMAT ( 6 X,2H  C,I4,2H  R, 1 4, 1H-, F1 1 .5) 

169  FORMAT ( 6 X,2H  C,I4,2H  R, 1 4, FI 2.5) 

171  FORMAT (6X,2H  C,I4#6H  OFSET,F12.5) 

172  FORMATS 2 1 5) 

173  FORMAT (1H*, 5 X,47H 

174  FORMAT? 7HF I RST  B) 

175  FORMAT? 12X.14H  R  1  2.0) 

176  F0RMAT(3HE0F) 

177  F0RMAT(12X,2H  R, 1 4, 1H  , F1 1 .5 ) 

178  F0RMAT(28HENTER  NEW  DATA  -  PRESS  START) 
183  FORMAT (FIO.4,315) 

END 


SUBROUTINE  PINCH  (YD, JFK,KFJD) 

I F(YD) 706, 707,707 

706  YO— YD 

GO  TO  (752  ,  753)  »KFJD 

752  PUNCH  168, JFK, YD, JFK 
GO  TO  754 

753  PUNCH  166, JFK, YD, JFK 

754  YD— YD 
RETURN 

707  60  TO  (755,756),KFJD 

755  PUNCH  166, JFK, YD, JFK 
RETURN 

756  PUNCH  168, JFK, YD, JFK 
RETURN 

166  F0RMAT(6X,2H  C,I4,6H  0FSET,F12.5/6X,2H  C, 

.0) 

168  FORMAT (6X,2H  C,I4,7H  0FSET-,F1 1 .5/6X,2H  C 

1. °)  ,  „ 

183  F0RMAT(F10.5,3I5) 

END 


4,6H  R 
I4,6H  R 


SUBROUTINE  SIGN  (JW,IW,WAM) 
IF(WAM)370,373,371 

370  WAH— WAM 

PUNCH  169,JW,IW,WAM 
RETURN 

371  PUNCH  167.JW, IW.WAM 
373  RETURN 

169  FORMAT (6X,2H  C,I4,2H  R,i4,F12.5) 

167  F0RMAT(6X,2H  C,I4,2H  R, 1 4, 1H-, F1 1 .5) 
END 
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B.  SMOG  2 


This  program  produces  a  surface  matrix  without  profile  conditions 
in  the  dual  formulation  and  X  fit.  The  surface  will  be  double 
splined  in  the  x  direction  and  may  be  either  single  or  double 
splined  in  the  z  direction.  The  program  is  written  in  FORTRAN  II 
for  an  IBM-1620  computer. 


Organization 

SMOG  2  consists  of  two  passes  as  shown  in  Fig.  B-l.  The  input 
data  shown  below  is  read  into  Pass  1  .  At  the  termination  of  the 
execution  of  Pass  1  the  data  necessary  for  Pass  2  is  punched.  This 
data  can  be  read  directly  into  Pass  2  .  If  it  becomes  necessary  to 
change  the  sign  of  a  second  difference  as  mentioned  in  Step  2,  this 
may  be  done  by  altering  the  Pass  2  data.  More  information  on  this 
data  is  found  in  the  description  of  the  Pass  2  input. 
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B.  SMOG  2 


PASS  1 

OPERATING  INSTRUCTIONS 


Pass  1  of  SMOG  2  calculates  the  second  differences  of  the  offsets 
and  the  Intersections  of  the  profile  with  the  waterllnes  and 
stations.  It  then  punches  the  row  Identifications,  deviation 
constraints  and  data  for  Pass  II.  The  program  Is  written  in 
FORTRAN  II  for  the  IBM-1620  computer. 


Fortran  Symbol  Definition 


M  -  Total  number  of  waterlines  in  the  surface  to  be 

faired.  This  will  be  an  odd  number  if  the  surface 
Is  double  splined  in  the  z  direction. 

N  -  Total  number  of  stations  in  the  surface.  This 

will  always  be  an  odd  number. 

X(J)  -  Distance  along  the  x  axis  from  the  origin  of 

the  surface  to  Station  j  . 

Z(I)  -  Distance  along  the  z  axis  from  the  origin  of 

the  surface  to  waterline  i  . 

Y(I,J)  -  Offset  of  the  surface  at  Waterline  i  ,  Station  j  , 
in  feet  and  decimals 


IAY  -  Number  of  feet  in  the  offset  if  read  in  feet, 
inches  and  eighths. 

IBY  -  Number  of  inches  corresponding  to  IAY 

ICY  -  Number  of  eighths  corresponding  to  IAY 

HU  -  Program  indicator 

If  NU  ■  -1,  scale  factors  for  normalizing  and 
translating  X(J)  and  Z(I)  are  read  in. 

If  HU  ■  41,  scaling  and  translating  are  done  by 
the  program  as  given  in  Section  I  of  this  Appendix. 
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SX  -  Scale  factor  for  X(J)  when  these  are  read  in 

TX  -  Translation  factor  for  X(J)  when  these  are  read  in 

SZ  -  Scale  factor  for  Z(I)  when  these  are  read  in 

TZ  -  Translation  factor  for  Z(X)  when  these  are  read  in 

RHS  -  Right  hand  side  -  In  the  dual,  this  is  really  the 

cost  on  the  equation  coefficients.  The  ideal  value 
for  RHS  would  be  zero;  however,  this  makes  the 
solution  susceptible  to  cycling  so  that  a  value 
of  .0001  is  usually  used. 

CURV  -  Value  of  the  requirement  vector  for  the  curvature 
constraints,  this  value  actually  goes  in  the  cost 
row  in  the  dual  and  is  always  made  zero. 

NAB  -  Program  indicator 

If  NAB  *  -1  ,  the  surface  will  be  double  splined 
in  the  z  direction  and  M  must  be  an  odd 
number. 

If  NAB  ■  +1  ,  the  surface  will  be  single  splined 
in  the  z  direction  and  M  is  unrestricted. 

COST  -  The  penalty  assigned  to  X  .  This  will  appear  in 

Row  1  of  the  right-hand  side.  The  value  is  normally 
made  20  . 

EXP  1  -  The  value  of  p  in  the  profile  function,  usually 

made  .3333  . 

EXP  2  -  The  value  of  r  in  the  profile  function,  usually 

3.0  . 

POINT  -  The  value  of  D  in  the  profile  function.  This 

number  is  the  length  of  the  fraction  of  the  station 
spacing  over  which  the  profile  function  is  effective 
in  terms  of  the  full-scale  surface  (the  value  will 
be  translated  and  normalized  exactly  like  the  x 
values  are). 

C(I)  -  Coefficients  of  the  Theilheimer  equation  describing 
the  profile.  There  must  be  (M+2)  coefficients. 
Straight  lines  or  other  curves  may  be  used  by 
making  all  except  the  desired  coefficients  zero. 
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The  input  data  cards  are  described  below.  There  are  certain  input 
variables  on  which  there  are  restrictions  on  the  value  which  can 
be  entered.  These  restrictions  follow. 


Variable  Minimum  No.  Maximum  No. 

N  3  15 

M  3  11 


In  describing  the  various  data  cards,  the  actual  FORTRAN  format 
field  is  used  in  most  cases.  These  fields  come  consecutively 
across  the  card  with  no  gaps  or  blank  columns  between,  except  where 
indicated.  The  field  descriptions  are  the  FORTRAN  F  field,  which 
uses  the  FORTRAN  fixed  point  decimal  number,  and  the  I  field 
that  uses  the  FORTRAN  Integer  number,  which  is  always  right 
justified.  The  card  numbers  are  not  punched  on  the  data  cards. 

Card  No. 

13-53  1 

Any  alphameric  problem 
heading  (This  heading  will 
be  used  for  identifying  the 
problem  on  the  LP  output) 


Format 

15 

15 

F11.6 

F11.6 

2 

Variable 

NAA 

NAB 

POINT 

COST 

* 

Format 

15 

15 

F11.6 

F11.6 

F11.6 

F11.6  3 

Variable 

M 

N 

CURV 

RHS 

EXP1 

EXP2 

Format  F10.4 

Variable  C(I) 


Next 

M+2  cards 


Contents  of  Card 

Card  Columns  1  2-12 

Contents  *  Blank 
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Format  F10.4 

Variable  Z (1) 


Alternate  Cards : 

Format  F10.4  F10.4 

Variable  X(J)  Y(I,J) 


or: 


N  Cards 
per  Group 


Next 

M 

Groups 


Format  F10.4  15  15  15 

Variable  X(J)  LAY  IBY  ICY 


This  arrangement  can  be  stated  as  follows: 


The  first  card  of  each  set  gives  the  height  z  of  the  waterline. 
The  following  cards  each  give  a  station  coordinate  and  the  offset 
at  the  intersection  of  the  given  waterline  and  station.  There  is 
one  set  of  cards  for  each  waterline  in  the  surface. 


Format  F10.4  F10.4  F10.4  F10.4  Last  card 

Variable  SX  SZ  TX  TZ 


(This  card  is  entered  only  if  NAA  ■  -1 


Output 

The  output  from  this  program  is  designed  to  be  read  directly  by 
the  If-90  program  after  card  to  magnetic  tape  conversion. 

The  first  card  punched  contains  the  Loading  Read  in  the  first  input 
card.  The  next  set  of  cards  are  the  row  identification  coordinates. 
These  identify  each  row  of  the  matrix  to  be  used  in  this  solution. 
For  this  matrix  all  of  the  rows  are  always  used  and  identified, 
including  the  cost  row  which  has  the  name  "OFSET." 

Next,  a  card  punched  "MATRIX"  is  produced,  followed  by  the  deviation 
constraint  portion  of  the  matrix.  Finally  the  input  data  for  Pass  2 
is  punched. 
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The  output  formats  are  shown  in  Fig.  B-  1  ,  and  the  order  of  the 
output  deck  in  Fig  B-2  .  The  actual  output  from  PASS  1  is  shown 
in  the  sample  problem. 


Sense  Switch  Settings 

Sense  Switch  1  -  Not  used 

Sense  Switch  2  -  ON  The  values  of  the  x  and  z  second 

differences  and  the  resulting  signs 
are  printed. 


Sense 

Switch  2 

-  OFF 

The  above  is  not  printed 

Sense 

Switch  3 

-  ON 

The  offsets  are  read  in  as 
inches,  and  eighths 

feet, 

Sense 

Switch  3 

-  OFF 

The  offsets  are  read  in  as 

feet  and 

decimals 


SAMPLE  PROBLEM 

The  following  sample  problem  is  a  five-station  by  five-waterline  set 
of  offsets  taken  from  the  DLG-26  class  ships.  The  stations  are 
Numbers  0  ,  1/2  ,  1  ,  1-1/2  ,  and  2  ,  and  the  water lines  are  12,  16, 

24,  and  28  . 

The  surface  contains  the  bow  profile  of  the  DIG.  The  matrix  is 
doifele  splined  in  the  z  direction.  Only  the  output  from  Pass  1 
is  listed  here;  the  output  from  Pass  2  is  included  with  the  instructions 
from  that  program. 
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SAMPLE  PROBLEM  -  SMOG  2,  PASS  1 


*  *  SAMPLE  INPUT  *  * 

*  TODD  DLG  5X5  END  PROF  I  LE  X-Z  DBL  SPL- 

1  -1  1.  20. 

550.  .001  3.  .33333333 

1.52935  1 

-  .5873  | 

End  Profile  Coefficients 


.0000 

.0000  .000000 

12.7500  1.041666 

25.5000  2.458333 

38.2500  3.874999  l  Surface  Offsets 
51.0000  5.333333 
1 .0000 

.0000  .000000 

12.7500  1.239583 

25.5000  2.697916 

38.2500  4.197916 

51.0000  5.812500 
2.0000 

.0000  .078 

12.7500  1.479166 

25.5000  2.979166 

38.2500  4.583333 
51.0000  6.333333 
3.0000 

.0000  .343749 

12.7500  1.791666 

25.5000  3.447916 
38.2500  5.250000 
51.0000  7.166666 
4.0000 

.0000  .708333 

12.7500  2.416666 

25.5000  4.322916 

38.2500  6.416666 
51.0000  8.583333 


l.O. 0-2 


Bb  -91 


i 


TYPEWRITER  LISTING  OF  INTERSECTIONS  AND  SECOND  DIFFERENCES 
(Switch  2  ON  )  OUTPUT 

*  TODD  DLG  5X5  Ft-iD  PROFILE  X-Z  D3L  SPL- 


END  PROFILE  !  IITERSECT I C!!3 


WL  LIES 

1 

1.529350 

- 

5873 

C  .0000 

2 

.S 42050 

— 

5873 

0.0000 

3 

.354750 

- 

5873 

0.0000 

4 

-.232550 

- 

5873 

C.0000 

5 

-.819850 

- 

5873 

0.0000 

STATIONS 

1 

2.604003 

2 

.901367 

3 

0.000000 

4 

0.000000 

5 

0.000000 

2ND  D ! FF  -  WL 

1 NES 

Wl. 

3TA 

NORM  VALUE 

SIGN 

1 

1 

OUT  OF  SURFACE 

1 

2 

OUT  OF  SUR 

FACE 

1 

3 

-5.176766 

1.00 

1 

4 

.041668 

-1  .00 

1 

5 

.041668 

-1 .00 

2 

1 

OUT  OF  SURFACE 

2 

2 

-37.680848 

1.00 

2 

3 

.041667 

-1.00 

2 

4 

.114584 

-1.00 

2 

5 

.114584 

-1.00 

3 

1 

OUT  OF  SURFACE 

3 

2 

-.963248 

1.00 

3 

3 

.104167 

-1  .00 

3 

4 

.145833 

-1.00 

3 

5 

.145833 

-1.00 

4 

1 

-.049093 

1.00 

4 

2 

.208333 

-1 .00 

4 

3 

.145834 

-1 .00 

4 

4 

.114582 

-1 .00 

4 

5 

.114582 

-1.00 
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LA  LA  LA  U>  LA 


1  .927938  -1.00 

2  .197917  -1.00 

3  .187500  -1.00 

4  .072917  -1.00 

5  .072917  -1.00 


2ND  01 FF-STATIONS 

STA  WL  NORM  VALUE  SIGN 

1  2  OUT  OF  SURFACE 

1  3  OUT  GF  SURFACE 

1  4  -.439121  1.00 

1  5  -.439121  1.00 

2  2  -3.216759  1.00 

2  3  .072917  -1.00 

2  4  .312500  -1.00 

2  5  .312500  -1.00 

3  2  .041667  -1.00 

3  3  .187500  -1.00 

3  4  .406250  -1.00 

3  5  .406250  -1.00 

4  2  .062500  -1.00 

4  3  .281250  -1.00 

4  4  .499999  -1.00 

4  5  .499999  -1.00 

5  2  .041666  -1.00 

5  3  .312500  -1.00 

5  4  .583334  -1.00 

5  5  .583334  -1.00 
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TODD  OLG  5X5  END  PROFILE  X-L  DBL  SPL-12/3/62 
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subroutine 


I CM  (YD,  JFK.KFJD) 

,  YD= OFFSET 

12k - -  JFK-  COL.  NO- 

YD  -  -  YD  KFJO-  A  SWITCH 


/\  punch  yo, 

CPJD^p^JFK  FOR 

ROW  2. 


PUNCH  YD 
JFK  For 
Row  i 


SUBROUTINE  PROFL  2.  XOFftEXPOl  T2PDJTZSOt  EXPD2,*f,T  D, 

TPD,  TSD,  NCO 

©THIS  SUBROUTINE  TAKES  THE 

1“  ( ***  derivative  or  the 

EXPRESSION 


N CO  *  2. 
TO  »  I 


TQ.j^POj 


TFfTOr 


T-  I 


POINTS  FZ-X 
POINT 


FOR  THE  PROFILE  CONDITION 
WITH  INFINITE  SLOPE 


TPD*  O 
TSDsO 


D1FD=  EX*>DI  -  I 
PROO‘TD##  01 FD 
TPOi  »expoi*  tzpo*PRoo 


SUBROUTINE  DIF  2.ND  COMPUTES  THE  DIFFERENCE 
OF  A  POINT  Xi,j,  Yi,j  BYJ 

n  -  *•  -  foii IIL  .  Vi 

Xt+,-Xi-,  lXu/-xi  xf-Vi-iy 


ft  Or  O 


The  vjttue  of  the  2-“  Diff* 
it  not  used,  interest  is  only 
in  the  sijn  of  umi  *•  a 
cheek  is  made  on  the  ^ '1919  of 
2**  6 iff.  in  this  subroutine 
and  when  the  sign  is  minus 
a  +1  value  L»  output  and  “I 
when  the  &.•&*■  diff.  1*  positive. 
When  zero  a  terc  is  output/ 


UD=At 


ojliA.  .U 
jxort-xo  XD 


XOMfl 


j 


;  ) 
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SUBROUTINE  DERIVI  (c2D,  C3D >,  C4D,  FXD,  F2PD,F2SD) 


THIS  SUBROUTINE  TAKES  THE  |H  £  DERIVATIVE 
OF  THE  STANDARD  CUBIC 


1.0. 0-2 
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SUBROUTINE  DERIVE (C5D,DIFD,FPD,  FS^FEPgFlS^ 


THIS  SUBROUTINE  TAKES  THE  1*1  4  2*® 
DERIVATIVE  OF  THE  THEJLHEIMER  TERMS 
AND  ADDS  THEM  TO  THE  STANDARD  CUBIC 
1ST  4  2#*  DERIVATIVES 


FCRj  LISTING 


C  SH0G2  PASS  1  DUAL  X  Oft  X  AND  Z  DDL  SPL  WITH  PROFILE 

DIMENSION  RX(11,15),RZ(15,11),TRM(150),  IRW(250),FZ(11),FX(15 

j  DIMENSION  X(15),Z(t1),Y(11,15),WL(15,11),IWL(15),C(13),FPZ(11 

CFSZ (11 ) 

C  *  *  INPUT  *  * 

750  READ  173 
PRINT  173 

READ  150, NAA, NAB, POINT, COST 
READ  150,M,N,CURV,RHS,EXP1,EXP2 
ML-M+2 

DO  751  1-1, ML 

751  READ  154,C(I) 

NCOMP-1 

NL-N+1 

ML-Mtl 

NS-N-1 

MS-M-1 

DO  2  1-1, M 
READ  154, Z( I ) 

DO  2  J-1 ,N 

I F(SENSE  SWITCH  3)802,803 

802  READ  183,X(J), IAY, IBY, ICY 
AAY-IAY 

BBY-IBY 

CCY-ICY 

Y(l,J)«AAY+(BBY+(CCY/8.))/12. 

GO  TO  2 

803  READ  154,X(J),Y(I,J) 

2  CONTINUE 
IF(NAA)3,3,4 

3  READ  154,SX,$Z,TX,TZ 
GO  TO  5 

4  SX-X(2)-X(1 ) 

SZ-Z  2)-Z(1) 

TZ-Z  1) 

TX«X(1) 

5  I F( SENSE  SWITCH  2)509,508 
509  PRINT  165 

C  *  *  NORMALIZING  *  * 

508  DO  26  J-1.N 

26  X(J)-(X(J)-TX)/SX 
DO  27  1-1. M 

27  Z( I )-(Z( I )-TZ)/SZ 
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c 


POINT-(POINT-TX)/SX 

*  *  ENO  PROFILE  INTERSECTIONS  *  * 

DO  694  1-1, M 

FX( I )-C( 1 )+C(2)*Z( l)+C(3)*Z( I )*Z( I )+C(4)*Z( l)**3 
CALL  OERIV1  (C(2) ,C(3) ,C(4) ,Z( I ) ,FPZ( I ) ,FSZ( I ) ) 

I F( 1-3)690,691,691 
691  DO  700  J-3.1 
DIF-Z(I)-Z(J-1) 

FX( I )-FX( I )+C( J+2)*D I F**3 

CALL  0ERIV2  (C(J+2) ,01 F,FPZ( I ) ,FSZ( I ) ,FPZ( I ) ,FSZ( I 
I F( 1-3)690,690.700 
700  CONTINUE 
690  FPZ( l)«FPZ( l)/SZ 

FSZ( I )-FSZ( I )/(SZ*SZ) 

I F(SENSE  SWITCH  2)693,694 

693  PRINT  170, I ,FX( I ) 

PRINT  154,FPZ(I),FSZ(I) 

694  CONTINUE 
J-1 

OLD  -Z(M) 
l-HS 


I F(SENSE  SWITCH  2)204,510 

204  PRINT  179 

510  FZ(J)-FX(I) 

ZZ-Z(I) 

I COMP- 1 

695  IF(ABSF(X(J)-FZ(J))— ,0001 )  66,66,71 

71  60  TO  (65,69). I  COMP 

65  IF(X(J)-FZ(J))  67,66,72 

66  FZ(J)-ZZ 
IF(SENSE  SWITCH  2)  74,80 

74  PRINT  170,J,FZ(J) 

80  J-J-H 
OLD-ZZ 
GO  TO  510 

67  ZZ-ZZ4-(0LD-ZZ)/2, 

I COMP-2 
IJUL-I-H 

68  FZ(J)-C(1)+C(2)*ZZ+C(3)*ZZ*ZZ+C(4)*ZZ**3 
IF( I JUL-3)  695,696,697 

696  FZ(J)-FZ(J)+C(5)*(ZZ-Z(2))**3 
60  TO  695 

697  DO  699  K-3.IJUL 

699  FZ(J)-FZ(J)+C(K+2)*(ZZ-Z(K-1))**3 


)) 
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69 

70 


72 

81 


83 

82 


6 

7 

8 
10 

660 

9 

513 

1 

13 

14 

15 

16 

17 

18 
21 


60  TO  695 

IF(X(J)-FZ(J))  67,66,70 

SEG-OLD-ZZ 

OLD-ZZ 

ZZ-ZZ-SEG/2. 

GO  TO  68 
OLD-Z(I) 

1-1-1 

IF( 1-1)81 ,510,510 
DO  82  l-J.N 
FZ(l)-Z(l) 

I F( SENSE  SWITCH  2)83,82 
PRINT  170, I ,FZ( I ) 

CONTINUE 

*  *  WATERLINE  2ND  DIFFERENCES  * 

I F(SENSE  SWITCH  2)6,7 

PRINT  157 

DO  15  1-1 ,M 

DO  13  J-1.NS 

60  TO  (8,513), NCOMP 

ZY-O. 

I  F(X(J)-FX(  0)10,10,9 
R-20. 

I F(SENSE  SWITCH  2)660,13 
PRINT  151, 1, J 
GO  TO  13 
ZZ-FX(I) 

NCOHP-2 
GO  TO  1 
ZZ-X(J-I) 

ZY-Y( I ,J— 1) 

CALL  0IF2ND  (X(J) ,X(J+1 ) ,ZZ,Y( I , J) ,Y( I ,J+1) ,ZY, I , 
RX( I ,J)«R 
RX( 

II 

PRINT 
NCOMP- 1 

*  *  STATION  2ND 

I F( SENSE  SWITCH  2)16,17 
PRINT  159 
DO  25  J-1 ,N 
DO  23  1-2, MS 
GO  TO  (18, 19), NCOMP 
IF(Z(I)-FZ(J))21,21,20 
R-20. 


tX(l,N)-RX(l,NS) 

I F ( SENSE  SWITCH  2)14,15 
>RINT  158, I ,J,A,RX( I ,N) 


DIFFERENCES 


J,R,A) 
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I F(SENSE  SWITCH  2)22,23 

22  PRINT  151, J, I 
60  TO  23 

20  ZZ-FZ(J) 

NCOHP-2 
GO  TO  514 
19  ZZ-Z(I-I) 

514  CALL  DIF2N0  (Z( I ) ,Z( 1+1 ) ,ZZ,Y( I ,J) ,Y( 1+1 ,J) ,Y( 1-1 , J) , J, I ,R,A) 

23  RZ(J,lj-R 
l  J*M). 


RZ 


-RZ(J,HS) 


i\4iv  •  n/ -f\i-  / 

I  F(SENSE  SWITCH  2)24,2! 


ROW  ID 


»2? 

24  PRINT  158, J, I ,A,RZ(J,H) 

25  NCOHP-1 

C  *  * 

PUNCH  173 
NJUL-((N-3)/2)+4 

NH.N 

IF(NAB)719,718,718 

719  JUL-3 
JUL2-2 

NOROW-NJUL*( (H-3)/2+4)+1 
GO  TO  720 
718  JUL-2 
JUL2-1 

N0R0W-(Mt2)*NJUL+1 

720  PUNCH  160 
DO  28  NR- 1, NOROW 

28  PUNCH  162, NR 
C  *  * 

PUNCH  164 

C  *  * 

MLG-(2*NM) 

1-1 
FUPZ— 1  • 

FUSZ-O. 

DO  37  J-1,MLG,2 
IF(RX(1,I)-20.)250,251,250 
251  K-2 

GO  TO  252 
250  K-1 

200  GO  TO  (201 ,702).NC0MP 

201  CALL  PR0FL2(X( l),FX(1),EXP1 , FUPZ, FUSZ,EXP2, POINT, T, TP, TS, 

IF(T-.0001)703,703,257 
257  ZY-Y(1,I)/T 


MATRIX 


FIRST  WATERLINE 
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703 

702 

90 

252 


253 

31 

32 

133 

33 

254 

255 

34 

35 


259 

705 


260 

36 

37 


F|-|T)34,34,133 
F(IT-2*(lf/2))33, 34.33 
AM(IR0W)-(X(I)-X(IT))**3 


202 

91 

HP) 


GO  TO  90 
NCOMP-2 
ZY-Y(I.I) 

KFJ-1 

CALL  PINCH  (ZY.J.KFJ) 

I  ROW-3 

DO  31  NPR-1.3 
TRM( IROW)-x( I )**NPR 
GO  TO  <253,31) .K 
PUNCH  167,J,IR0W,TRM(IR0W) 

I  ROW- 1 ROW+1 
I F( 1-3)35.35,32 
DO  34  I T-3 , I 
IF(I-IT)34, 

I 

TRH(IROW). 

GO  TO  (254.255). K 
PUNCH  167,J,IR0W.TRM(IR0W) 

I  ROW- 1  ROW+1 

CONTINUE 

JPLS-J+1 

I  ROW- 1  ROW-1 

IWL( I )-IROW 

GO  TO  (259, 705), K 

KFJ-2 

CALL  PINCH  (ZY.JPLS.KFJ) 

DO  36  IR-3, I ROW 
WL(I.IR)-TRH(IR) 

GO  TO  (260. 36). K 
PUNCH  169,JPLS.IR.TRM(IR) 

CONTINUE 

1-1+1 

NCOHP-1 

*  *  REMAINING  WL  *  * 

DO  64  IZ-2.H 
JL-J 

M6*JL+( 2*NM)-2 
IX-1 

DO  63  J-JL.HG.2 
I F— 2 

I  RW(  IF)- 1  ROW+1 

IF(RX( IZ. I X)-20() 202,63 ,202 

00  TO  (91.92).NC0MP 

CALL  PR0FL2(X( I X) ,FX( I Z) ,EXP1 ,FUPZ,FUSZ,EXP2,P0INT,T,TP,TS,NC0 
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258 

93 

92 

95 

708 

41 


42 

142 

729 

730 
144 


145 

146 
709 

711 

43 

44 


244 


IF(TT.0001)93,93,258 

ZY-Y(1Z,IX)/T 

60  TO  95 

NCOKP-2 

ZY-Y(IZ,IX) 

KFJ-1 

CALL  PINCH  (ZY,J,KFJ) 

I F( J-JL)41 ,41 ,44 

00  42  IZPR-1,3 

TRM(IF)-Z(IZ)**IZPR 

PUNCH  167, J, I RW( IF) ,TRM( IF) 

IF-IF+1 

I RW{ I F)-IRW( I F-1 )+NJUL 
IF(IZ-JUL)146,146,142 
DO  145  ITZ-JULjlZ 


,144 


IF( IZ-ITZ)146,146,729 
GO  TO  (144,144.730),JUL 
IF(ITZ-2*(ITZ/2))144.145, 
ZTR-(Z(IZ)-Z(ITZ))**3 
TRH(IF)-ZTR 

PUNCH  167, J, IRW( IF),TRH( IF) 
IF-IF+1 

IRW( IF)-IRW(IF-1)+NJUL 

CONTINUE 

JPLS-J+1 

KFJ-2 

CALL  PINCH  (ZY,JPLS,KFJ) 

JF-IF-1 

DO  43  IF-2.JF 

PUNCH  169,JPLS, IRW( IF) ,TRM( IF) 
60  TO  63 

*  *  Z  ,  Z  SQR  ,  Z  CUBE 

IRO  -2 

K-IWL(IX) 

DO  244  IJ-3,K 
IRO  -IRO  +1 

PUNCH  167,J( IRO,WL( IX, I J) 

DO  50  IZPR-1,3 
ISTRT-IRW( IF;+NJUL 
ZTRH»Z( IZ)**IZPR 
TRM( IF)-ZTRM 

PUNCH  167,J, IRW( IF) ,TRH( IF) 
IF-IF+1 

IRW(IF)-IRW(IF-1)+1 


ROWS  *  * 


1.0. 0-2 


B-146 


( 


45 

46 

ft 


49 

50 


731 

55 

56 


256 

ll 


59 

60 
61 

712 


DO  45  NPR-1,3 
TRM(IF)-(X(IX)**NPR)*ZTRM 
PUNCH  167*J» IRW( IF) ,TRM( IF) 
IP-IF+1 

IRW(IF)-IRW(IF-1)+1 

IF( IX-3)50,50,46 

DO  49  IT-3. IX 

IF(IX-IT)49.4?.47 

IF( IT-2*( IT^2) )4oj49|48 

TRH( I F)-< (X( I X)-X{ IT))**3)*ZTRM 

PUNCH  167,JJ«W(IF),TRH(IF) 

IF-IF+1 

IRW(IF)-IRW(IF-1)+1 

CONTINUE 

IRW(IF)-ISTRT 

*  *  THE  I LHE I HER  ROWS 

IF(IZ-JUL)61,61£54 
DO  60  ITZ-JUL.IZ 
GO  TO  (55,55,731),JUL 
IF(ITZ-2*(ITZ/2))55,60,55 
ISTRT-IRW(IF)+NJUL 
IF(IZ-ITZ)61.61,56 
ZTRM-(Z(IZ)-i(lTZ))**3 
TRH(IF)-ZTRM 

PUNCH  167,J,IRW(IF),TRH(IF) 
IF-IF+1 

IRW(  IF)-IRW(  IF-D+1 

DO  256  NPR-1,3 

TRH( I F)-( X( I X)**NPR)*ZTRM 

PUNCH  167,J,IRW(IF),TRM(IF) 

IF-IF+1 

IRW(  IF)-IRW(  IF-1)+1 

DO  59  IT-3. IX 

IF(IX-IT)59.59.57 

IF(IT-2*(IT/2)j58,59,58 

TRH( I F)-((X( I X)-X( I T) )**3)*ZTRH 

PUNCH  167  ,J, IRW( IF) tTRM( IF) 

IF-IF+1 

IRW(  IF)-IRW(  IF-D+1 

CONTINUE 

IRW(IF)-ISTRT 

JPLS-J+1 

KFJ— 2 

CALL  PINCH  (ZY,JPLS,KFJ) 
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< 


714  IRO  -2 
K-IWL(IX) 

DO  262  IJ-3.K 

mo  -mo 

262  PUNCH  169,JPLS,IR0  ,WL(IX,IJ) 

JF-IF-1 
DO  62  IS-2,JF 

62  PUNCH  l69,JPLS,mW(IS),TRM(IS) 

63  IX-IX+1 

64  NCOMP-1 

PUNCH  180,JPLS,JUL,JUL2,NJUL,M,N,MS,NS,N0R0W,NC0MP,CURV,RHS 
PUNCH  181, EXP2, COST, POINT, EXP1 
DO  87  1-1, N 
PUNCH  181,X(I),FZ(I) 

DO  87  J-2,M 
87  PUNCH  181 ,RZ( I ,J) 

DO  88  1-1, H 

PUNCH  1 81 ,Z( I ),FX( l),FPZ( I) ,FSZ( I ) 

DO  88  J-1 , N 
88  PUNCH  1 81 ,RX( I ,J) 

ML-H+2 

DO  100  J-1 ,  ML 
100  PUNCH  181 ,C( J) 

PRINT  178 
PAUSE  751 
60  TO  750 

150  F0RMAT(2I5,4F11.6) 

151  FORMAT* 2 1 5, 1?H  OUT  OF  SURFACE) 

154  F0RMAT(4F10.4) 

157  F0RHAT(//10X, 1 7H2ND  DIFF  -  WLINES// 

C30H  WL  STA  NORM  VALUE  SIGN/) 

158  F0RMAT(2I5,F12.6,F8.2) 

159  FORMAT ( // 10X , 1 7H2ND  DIFF-STATIONS// 

C30H  STA  WL  NORM  VALUE  SIGN/) 

160  FORMAT 
162  FORMAT 

164  FORMAT 

165  FORMAT 
167  FORMAT 

169  FORMAT 

170  FORMAT 

171  FORMAT 

172  FORMAT 

173  FORMAT 


6HR0W  ID/12X,6H  OFSET) 

12X,2H  R. 14) 

6HMATR I X) 

//6X.25HEND  PROFILE  INTERSECTIONS//15X,6HWLINES/) 
6X,2H  C,I4,2H  R, I4,1H-,F11 ,5) 

6X,2H  C.I4.2H  R, 1 4, FI 2*5) 

5X, I5,F12.o) 

6X^2H  C,I4,6H  OFSET, FI 2,5) 

1H*.5X.47H 
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178  FORMAT ( 28HENTER  NEW  DATA  -  PRESS  START) 

4  M  ti  A  W  i  I  «  I  tj  ftl  1AM  ^  it  W 


IS 

181 


FORMAT 
FORMAT 

FORMAT  _ 

183  F0RMAT(F10.5,3I5) 
END 


/ 1 4X , 8HSTAT I ONS/ ) 
10I5.2F15.9) 
4F15.9) 
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(  ) 

SUBROUTINE  PINCH  (YD, JFK,KFJD) 

I F( YD) 706,707*707 

706  YD- -YD 

GO  TO  (752*753) ,KFJD 

752  PUNCH  168, JFK, YD, JFK, JFK 
GO  TO  754 

753  PUNCH  166, JFK , YD, JFK, JFK 

754  YD— YD 
RETURN 

707  GO  TO  (755*756)  ,KFJD 

755  PUNCH  166, JFK, YD, JFK, JFK 
RETURN 

756  PUNCH  168, JFK, YD, JFK, JFK 
RETURN 

166  F0RMAT(6X,2H  C,l4,6H  0FSET,F12.5/6X,2H  C,l4,6H  R  1,8H  1.0/ 

C6X,2H  C, I 4,6H  R  2,8H-  1.0) 

168  FORMAT(6X,2H  C,l4,7H  0FSET-,F1 1 .5/6X.2H  C,l4,6H  R  1,8H  1.0/ 

C6X,2H  C, 14, 6H  R  2,8H  1.0) 

END 

SUBROUT  I NE  PR0FL2 ( XD , FD , E  XPD 1 , TZPD , TZSD , E XPD2 , XI , TD , TPD , TSD , NCD ) 
TD«(X1+FD-XD)/X1 
I F(TD-. 00001)1, 1.2 

2  DIFD-EXPD1-1 . 

PR0D-TD**D I FD 

TPD 1 -E  XPD 1 *TZPD*PROD 

TSD1-EXPD1*(D I FD*TD**(D I FD-1 . )*TZPD**2+TZSD*PR0D) 

T0«1.-TD**EXPD1 

IF(TD-.00001)7,7,4 

4  D I F0-EXPD2-1 . 

TDEXP-TD**D I FD 
TPD— EXPD2*TPD1*TDEXP 

TSD-EXPD2*(D!FD*TD**(EXPD2-2.)*TPD1*TPD1-TSD1*TDEXP) 

TD-T0**EXPD2 
IF(TD-.00001)7,7*6 
7  TD-O. 

GO  TO  3 
6  NCD-1 
RETURN 
1  NCD-2 
TD-1. 

5  TPD-O. 

TSD-O. 

RETURN 

3  NCD-1 
GO  TO  5 
END 
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SUBROU T I NE  D I F2ND ( XO , XDP 1 , XDM1 , YD , YDP 1 , YDM1 , I D , JD , RD , AD ) 
AD-2./(XDP1-XDM1) 

ZA-YDP1-YD 

AD«AD*(ZA/(XDP1-XD)-(YD-YDM1 )/(XD-XDH1) ) 
:F(AD)800,o09>808 

808  RD— 1. 

GO  TO  tt 

809  RD-O. 

GO  TO  11 
800  RD-1. 

11  I F(SENSE  SWITCH  2)812, 813 

812  PRINT  198, ID, JD.AD.RD 
158  F0RMAT(2I5,F12.6,F8.2) 

813  RETURN 
END 


SUBROUTINE  0ERIV1  (C2D,C3D,C4D,FXD,FZPD,FZSD) 

FZPD-C  2D+2 . *C3D*F  XD+3 • *C  4D*F  XD*  F  XD 

FZSD-2.*C3D+6.*C4D*FXD 

RETURN 

END 


SUBROUTINE  DERIV2  (C5D,DIFD,FPD,FSD,FZPD,FZSD) 
PROD-3. *C5D*D I FD 
FZPD«FPD+PROD*D I FD 
FZSD-FSD+2.*PR0D 

RETURN 

END 
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SMOG  % 
PASS  2 


OPERATING  INSTRUCTIONS 

Past  2  of  the  program  accepts  the  data  punched  in  Pass  1.  It  produces 
the  curvature  constraints  of  the  matrix  and  then  the  right-hand  side  , 
which  for  the  dual,  consists  of  the  costs  on  the  variables. 


Fortran  Input  Symbols 


JPLS  -  The  number  of  the  last  column  produced  in  Pass  1 


JUL 


JUL2 


NJUL 


M 


Program  indicator 

JUL  >3  if  the  surface  is  double  splined  in 
the  z  direction 

JUL  “  2  if  the  surface  is  single  splined  in 
the  z  direction 

This  value  is  equal  to  the  subscript  of  the  first 
z  value  of  a  point  of  discontinuity 

Program  indicator 

JUL2  -  2  if  the  z  direction  is  double  splined 

JUL2  ■  1  if  the  z  direction  is  single  splined 

This  value  is  equal  to  the  number  of  waterline 
intervals  between  points  of  discontinuity 

This  value  +  4)  is  equal  to  the  maximum  number 

of  rows  necessary  to  store  the  x  terms,  each  of 
which  is  multiplied  by  a  specific  z  term 

Total  number  of  waterlines  in  the  surface 


MS  -  M-l 


N 


Total  number  of  stations  in  the  surface 


NS 

NOROW 


N-l 


Total  number  of  rows  in  the  matrix 


For  z  double  spline  NOROW 
For  z  single  spline  NOROW 


M-3 


1+  4  NJUL 


P+2)NJUl] 


+  1 
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NCOKP 


Program  indicator 


COST  -  The  penalty  assigned  to  X  •  This  will  appear  in 

Row  1  of  the  right-hand  side.  The  value  is  normally 
made  20. 


EXP1 

-  The  value 

.3333  . 

of 

P 

in  the  profile  function  usually  made 

EXP2 

-  The  value 

3.0  . 

of 

r 

in  the  profile  function,  usually 

POINT 

-  The  value 

of 

D 

in  the  profile  function.  This 

number  is  the  length  of  fraction  of  the  station 
spacing  over  which  the  profile  function  is  effective 
in  terms  of  the  full  scale  surface  (the  value  has 
been  translated  and  normalised  exactly  like  the  x 
values  are). 

C(l)  -  Coefficients  of  the  Theilheimer  equation  describing 
the  profile.  There  must  be  (Mf2)  coefficients. 
Straight  lines  or  other  curves  may  be  used  by  making 
all  except  the  desired  coefficients  zero. 

RHS  -  Right-hand  side  -  In  the  dual  this  is  really  the 

cost  on  the  equation  coefficients.  The  ideal 
value  for  RHS  would  be  zero;  however,  this  makes 
the  solution  susceptible  to  cycling  so  that  a  value 
of  .0001  is  usually  used. 

CURV  -  Value  of  the  requirement  vector  for  the  curvature 

constraints.  This  value  actually  goes  in  the  cost 
row  in  the  dual  and  is  always  made  zero. 

X(J)  -  Distance  along  the  x  axis  from  the  origin  of  the 

surface  to  Station  j  . 

Z(I)  -  Distance  along  the  z  axis  from  the  origin  of  the 

surface  to  Hater line  1  . 

FX(I)  -  The  x  coordinate  of  the  Intersection  of  Waterline 

i  with  the  profile 

FZ(J)  -  .The  z  coordinate  of  the  intersection  of  Station 

j  with  the  profile 

FPZ(I)  -  Tha  first  derivative,  with  respect  to  z  ,  of  the 

profile  equation  at  Waterline  i  . 

FSZ(I)  -  The  second  derivative  with  respect  to  z  of  the 

profile  equation  at  Waterline  i  . 
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&Z(I,J)  -  The  sign  of  the  second  difference  of  the  Station 

Curve  j  at  Waterline  i  .  This  will  have  the 
opposite  sign  from  the  actual  second  difference. 

RX(I,J)  -  The  sign  of  the  second  difference  of  Waterline 

Curve  i  at  Station  j  .  This  will  be  the  opposite 
of  the  actual  sign. 


Input  Data 


The  following  input  formats  are  listed  as  a  matter  of  record  in  case 
modifications  of  the  data  cards  are  to  be  made,  normally  this  data 
will  be  obtained  from  executing  Pass 1  in  a  form  which  can  be  entered 
directly  into  Pass  2. 


an tents  of  Card 


Card  No. 


Format 

15 

15 

15 

15  15 

15  15  15 

Variable 

JPLS 

JUL 

JUL2  NJUL  M 

N  MS  NS  1 

Format  (continuation 

of  Card  1) 

15 

15 

F15.5  F15.5 

Variable 

MORON 

NCCMP 

CURV  RHS 

Format 

F15.9 

F15.9 

P15.9 

F15.9 

2 

Variable 

EXP2 

COST 

POINT 

EXP1 

Format 

F15.9 

F15.9 

Variable 

X(I) 

FZ(I) 

1st  card/set 

y  Next 

Format 

F15.9 

N  Sets 

Variable 

RZ(I,J) 

M  Cards/set 

Format 

F15.9 

F15.9 

F15.9 

F15.9 

n 

Variable 

Z(I) 

FX(I) 

FPZ(I) 

FSZ(I) 

1st  card/set  i  Next 

Format 

F15.9 

M  Sets 

Variable 

RX(I,J) 

N  Cards /set 

J 

Format 

F15.9 

Next 

Variable 

C(J) 

M+2  Car 
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Output 


The  curvature  constraint  elements  are  produced  first.  These  are 
punched  with  the  formats  listed  In  Fig.  B-  ,  as  are  the  following 
cards.  The  elements  are  punched  In  column  order  with  the  rows 
ordered  within  the  columns.  Next,  a  card  punched  "FIRST  B"  Is 
produced.  Afterwards,  the  right-hand  side  is  punched  in  row  order. 
Finally,  the  last  card  in  the  deck  is  punched  "EOF." 

After  the  input  data  for  Pass  2  is  removed  from  the  Pass  1  output, 
the  Pass  2  output  is  placed  directly  behind  the  Pass  1  output.  The 
deck  is  now  complete  and  ready  for  input  to  LP-90  . 


Sense  Switches 


No  sense  switches  are  used. 
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0.00000000 


SAMPLE  PROOBUEM  -  Pats  2  Input  Punched  by  Past  1 


50  3  2 

.00100000 
.33333333 
0.00000000 
20.00000000 
20.00000000 
1.00000000 
1.00000000 
1 .00000000 
1.00000000 
-1 .00000000 
-1 .00000000 
-1.00000000 
2.00000000 
-1 .00000000 
-1 .00000000 
-1.00000000 
-1 .00000000 
3.00000000 
-1 .00000000 
-1.00000000 
-1 .00000000 
-1 .00000000 
4.00000000 
-1 .00000000 
-1 .00000000 
-1 .00000000 
-1 .00000000 
0.00000000 
20.00000000 
20.00000000 
1 .00000000 
-1 .00000000 
-1.00000000 
1 .00000000 
20.00000000 
1.00000000 
-1 .00000000 
-1.00000000 
-1 .00000000 
2.00000000 
20.00000000 
1 .00000000 
-1 .00000000 
-1 .00000000 
-1 .00000000 
3.00000000 


*  *  SAMPLE  INPUT  *  * 

5  5  5  4  4  26 

20.00000000  .07843137 

2.60400380 


.90136719 


0.00000000 


0.00000000 


0.00000000 


1.52935000  -.58730000 


,94205000  -.58730000 


.35475000  -.58730000 


-.23255000  -.58730000 


1 

3.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 
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1 .00000000 
-1 .00000000 
-1 .00000000 
-1 .00000000 
-1 .00000000 

4.00000000  -.81985000  -.58730000 

-1 .00000000 
-1 .00000000 
-1 .00000000 
-1 .00000000 
-1 .00000000 
1.52935000 
-.58730000 
0.00000000 
0.00000000 
0.00000000 
0.00000000 
0.00000000 


0.00000000 
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PASS  2  *  *  SAMPlc  OUTPUT 


( 


* 


oooooooooooooooooooooooooooooooooooooo 

OOOOOOOOOOOOOOOOOOOOQOOOOOOOOOOOOOOOOO 

oooooooooooooooooooooooooooooooooooooo 

oooooooooooooooooooooooooooooooooooooo 

oooooooooooooooooooooooooooooooooooooo 

O—  -OO  —  fM-OO  —  O—  -O'—  -O  —  -OO  —  —  CVI-C— cg-CO"—  -OCSJlNJ:».3’®X> 

—  —  —  —  —  —  (\J  -» 

UJ  ULI  UJ  —  —  ~  rg  jj  —  CNJlM'JJ  —  JM 

(/)</)</>  </>  u> 

u.  u.  u.  u.  u. 

Oo£a:oacot«OKaCa:aacocaco£00faccoc«ate:a£oeBCeCQCOaeQCa'a:ceorofa: 

<M(M(MiOfOPo>nininvninifii^iri/>ir>'0'0-o>0'0-0'0-oo-o<5'0'Oxiao®ooooco<r)ODoo 

oooooooooooooooooooooooooooooooooouooo 


S-/*7 


( 


o  o  o 
o  o  o 

o  o  o 
o  o  o 
o  o  o 


O  —  CM 


oooooooooooooooooo 

oooooooooooooooooo 

oooooooooooooooooo 

oooooooooooooooooo 

oooooooooooooooooo 


■OCMJMJ®4CQOaO-’«1>h-'*0- 
<N  —  -3  (M  O-  -*  CM 


o  o  o  o  o  o  o 
oo  o  o  o  o  o 
oo  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 


OIOBO 

— 


o 

o 

o 

o 

o 


o  o  o 
o  oo 
ooo 
o  o  o 
ooo 


o  o  — 


ooooooooo 

ooooooooo 

ooooooooo 

ooooooooo 

ooooooooo 


—  rO  —  O  if.  cm 


h-  in 

aJ 

m 

a. 

O  OC  DC 


—  —  —  —  —  CM  CM  UJ  —  <—  cm  at  —  «-  <—  —  CM  CM  CM  Ui  — 

mm  m 

u.  u.  u. 

acaacacacacacacacacoacQeacacacoacc'ocetKocacacaeacooeaeacaCQCeeacaeccor 


0-0'»l>»0'OOl>0000'000000’-'-'---,*'"’-’"'',",“tMlMN!MlMCMN(MCMtMC\ 

ininminminLninininininin<3-o<3'0>o<30-oooo-0'0'00'0'0'00'0'00-0'0-o-0'Oo 


OUOOOOUOOUOOOOOOOOOOOOOOOOOOOOOlJOOOOOUOOU 


J3  1 

I 

I 

*  I 


/.*.***• 


I 

I 


I 

I 

i 

I 


I 

I 

I 

I 


I  ( 


ooooooooooooooooooooooooooooooooooooooooo 

ooooooooooooooooooooooooooooooooooooooooo 

ooooooooooooooooooooooooooooooooooooooooo 

ooooooooooooooooooooooooooooooooooooooooo 

Q  Q  Q  Q  Q  Q  Q  O  Q  Q  O  Q  Q  Q  ooqoooooo  ooooooo  o  o  o  oooooooo 


»N-N'00-S'04®0--04»'0«44®050’-N'04®4'ON'Odn«ffl-OflDO-(M 
«M  O  —  —  O  CM  —  fr  -O  00  Jt  —  a  CM  —  o  >  -o  >o  oo  o  a  — 

**1  **  W  »-*  K  IO 


I 


I  CM  CM  CM  CM  CM  OJ  '•"•NU1  —  r-cvevtsiuj  —  CMCNCMCMCMUJ'—  ~ 

1  c/>  «/»  VO  t/> 

I  _ _  _  £  t  Ik  u. 


eiceocceacoacecacocacoececeeectieccacacacacoaeccaeeexacatacaconcieBCKctcicoacx 

"OOOO'OO-OO-OO-OO^O-O'O'O-O'O-O'O-O-O-O-O-OO-O-O'O'OOO-O'OO'O'O-O-O'OO 


OOcJOOUOUOOOUO<JOOU<JUOk>VCJc-MJUOOOOOOUOUU>k>OaOO 


I 

I 


I 


( 


3.IU 


I 


pgt_cg«)j}-tr>r-®(>o»-c\j»njfvnf>-coO‘OCMr»ij»'1nH-c\i'o^in^aoo-0(— cNjfO^riDf'-oooo 
—  —  CMt'JCM(NCVJUJ>— —  —  cmuj— -  —  —  --CSJ 

VJ  </>  t/»  VO 

u.  u.  a.  u. 

ccoecacacecacececacoxcccceceeacecotecacacccoeixoeaca'acacecoocecoiicxecececc* 

t^ascocoajcoaaaocDfloo-O'O'OO'O-o-oooO'O-oooooooooo  — — 
-o-O'O'O-O'O'O'O'O'C'OO'OO'O-O'O'O'O'O'O'O'Os-N-f'f-f'-N-r'-f-r-Kr-f-N-r'-f'-r^f'-h- 


U  O  UUOOOOUOOOOOO<JOOOOUOOOOOUUOOOOOOO<JOUO(JO 


t,0.o  *  »- 


0-/4#. 


71  R  22  6.00000 

71  R  23  12.00000 

71  R  24  24.00000 

71  R  25  48.00000 

72  OFSET  0.00000 


oooooooooooooooooooooooooooooooooooo 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOQOOOOOOO 

oooooooooooooooooooooooooooooooooooo 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

oooooooooooooooooooooooooooooooooooo 


cm  —mo  cm  — mocM—  — mo  — O  CM  O  00 

—  —  ro  —  OO 


(MMtfifl’ONffloo-  »-cMmjrmoN-ooO"0  —  (M'Ojfmoh-cgivt^moN.coO'O 
—  —  —  ——  ——  —  CMCMUJ  —  —  —  —  —  ——  —  CMCMCMCMCMCNCMUJ  —  —  ——  ——  —  —  CM 

«/»  m 

u.  u. 

ecacaca:ecoceC(XacocoecacacctecacecoCocacc£ccaCBC;a:oececo(:aca:o£a:cctx 


OOOOOUOOOUOOOOOUUOOlJOOOOOOUOOOUOOOUOUOOOO 
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oooooooooooo 

oooooooooooo 
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—  -o  - 


'-(-<\|NiJ,iri'0»^cooo  —  cMWjfvri'O 
<MOJ  —  —  PJCMPgfS|P<|(NfM 

l/) 

u. 
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96 .00000 
384.00000 
48.00000 
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ooooooooooooooooooooo 

ooooooooooooooooooooo 
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ececocaeaeoiarotaeactta'aracococococoeacesasoc 
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00100 

00100 
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B— 171 
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For  Flow  Diagram  of  Profile  2  Subroutine,  See  SMOG  2,  Pass  1 


SUBROUTINE  SEDIF l(TD/TPD/TSD<1  XD, JROWD,  JD,PWRl*IROW^ 


THIS  SUBROUTINE  COMPUTES 
THE  2fiS  DERIVATIVE  OF 
THE  STANDARD  CUBIC  USING 
SEDIF4>  SEDIFSy  ^  SEDIF 6 
SUBROUTINES 
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SUBROUTINE  SED1F2(dIFD,TD/TPD,TSD,PWRD,IROWO, 


WRM-(3KDtFD*TD+3#TPD*WFD#*24£-SD*l>IFD**3)/2)  *  PWRD 

^QOTOSIQN^V  WRM  IS  THE  DERIVATIVE 
^SUBROUTINES  THE  THEILHEIMER  TERM 

rETMN 


subroutine  sedips^xo.fxp.prd.irowd.jsd) 


WRM«(kxD*FXD**-(lXD-l)-FXD**IXD)*  pro 


WRM  IS  THE  Zt&  DERIVATIVE 
FOR  POINTS  OKI  THE 
PROFILE 


1. 0.0-2 
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SUBROUTINE  SEDIF4-(tSD,PWRD,  JD,IR0WD) 


SUBROUTINE  SEOIFS  (TPD,  TSD,  XD„  PWRD,  JO,  IROWO) 
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SUBROUTINE  SEDIft^TDJPD.TSD.XD,  JD,  PWRD,IROWD) 


SUBROUTINE  SIGN  ( JW,  3  W,  WAM^ 


THIS  SUBROUTINE  PUNCHES 
THE  Z"  DERIVATIVE  FROM 
THE  SEDIF  SUBROUTINES, 
ALONG  WITH  THE  identifying 
SUBSCRIPT  NO'S  I  $  J,  IN 
PROPER  FORMAT 
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FORTRAN  LISTING 


C  SMOG 2  PASS  2  DUAL  X  OR  X  AND  Z  DBL  SPL  WITH  PROFILE 

DIMENSION  RX(11,15),RZ(15,11),FZ(11),FX(15),X(15),Z(11),Y(11 
15) 

DIMENSION  FPZ(11),FSZ(11),C(13) 

750  READ  180,JPLS, JUL, JUL2,NJUL,M,N,MS,NS,N0R0W,NC0MP,CURV,RHS 
READ  181,EXP2,C0ST,P0INT,EXP1 
DO  32  l-f ,N 
READ  181.X(I),FZ(I) 

DO  32  J-2.M 

32  READ  1 81 ,RZ( I ,  J) 

DO  33  I-1.M 

READ  l8l.Z(l),FX(l),FPZ(l),FSZ(l) 

DO  33  J-1.N 

33  READ  181 fRX( I ,J) 

ML-M+2 

DO  100  J-1.ML 
100  READ  181, C(J) 

C  *  *  X  SEC  DIFFERENCES  *  * 

JX-NS/2 
DO  88  IZ-1,M 
1-1 


JS-JPLS+t 

JF-JS+JX 

DO  87  J-JS.JF 

IF(RX(IZ, I )-20.)501 ,87,501 

501  PUNCH  171.J.CURV 
I  ROW- 2 

ISTRT-2 

60  TO  (503,504) ,NCOMP 

503  FUPZ— 1  • 

FUSZ-O. 

CALL  PR0FL2  (X( I ) ,FX( IZ) ,EXP1 ,FUPZ,FUSZ,EXP2, POINT, T, TP, TS 
,NCOMP) 

504  PWR-RX(IZ.I) 

CALL  SEDIF1  (T,TP,TS, X( I ) , IROW,J,PWR, IROW) 

IROW-IROW+1 

CALL  SEDIF2  (X( I ) ,T,TP,TS,PWR, IROW,J) 

I  F(  1-3)75,75,502 

502  IROW-IROW+1 

DO  505  K-3,l,2 
IF(K-I)73, 75,75 
73  DIF-X(I)-X(K) 

CALL  SEDIF2  (DIF,T,TP,TS,PWR, IROW,J) 

505  IROW-IROW+1 

75  IF(IZ-1)87,87,76 


» 
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76  00  79  NPWR-1,3 

I STRT-I STRT+NJUL 
IROW-ISTRT 

PWR-RX( I Z, I )*Z( 1 Z)**NPWR 

CALL  SED t FI  (T,TP,TS,X( I ) , IROW,J,PWR, IROW) 

IROW-IROW+1 

CALL  SEDIF2  (X( I ) ,T,TP,TS,PWR, IROW, J) 
IROW-IROW+1 
DO  78  K-3 ,1,2 
I F(K- I >77.79.79 

77  DIF-X( I )-X(K) 

CALL  SE0IF2  (01 F,T,TP,TS,PWR, IROW,J) 

78  IROW-IROW+1 

79  CONTINUE 
IF(IZ-JUL)87,87,511 

511  DO  86  L-JUL, IZ.JUL2 
IF(L-IZ)512,87,87 

512  I STRT-I STRT+NJUL 
IROW-ISTRT 

PWR-RX( IZ.I )*(Z( IZ)-Z(L) )**3 

CALL  SED I  FI  (T,TP,TS,X( l) ,IROW,J,PWR, IROW) 

IROW-IROW+1 

CALL  SE0IF2  (X( I ) ,T,TP,TS,PWR, IROW.J) 
IROW-IROW+1 
DO  85  K-3  *1.2 
IF(K-I >84.86.86 

84  DIF— X( I )-X(K> 

CALL  SEDIF2  (D I F,T,TP,TS,PWR, IROW,J) 

85  IROW-IROW+1 

86  CONTINUE 

87  1-1+2 
NCOMP-1 

88  JPLS-JF 

Z  SEC  DIFFERENCES  *  * 


DO  110  1-1  ,N 
DO  111  IZ-JUL.M.JUL2 
GO  TO  (550.551) .NCOMP 
550  IF(RX(IZ,l)-20.)556, 111,556 
556  PWR-RZ(I.IZ) 

FZD-Z(IZ) 

FOZ-FX(IZ) 

FUPZ-FPZ(IZ) 

FUSZ-FSZ(IZ) 

562  CALL  PR0FL2(X(I),F0Z,EXP1,FUPZ,FUSZ,EXP2,P0INT,T,TP,TS, 
NCOMP) 
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GO  TO  563 
551  PWR-RZ(I.IZ) 

FZD-Z(IZ) 

563  PUNCH  171.J.CURV 

GO  TO  (564,565) ,NCOMP 

565  ISTRT-2+NJUL 
GO  TO  566 

564  ISTRT-2 

CALL  SEDIF4  (TS,PWR,J, ISTRT) 

567  I F( 1-1)569,568,569 

569  IROW-ISTRT+1 

DO  570  IXPWR-1,3 
PROD«PWR*X( I )**l XPWR 
CALL  SEDIF4(TS,PR0D,J, (ROW) 

570  IROW-IROW+1 

I F( 1-3)568,568,571 

571  DO  572  K-3,1.2 
IF(K-I)604,568,604 

604  PR0D«PWR*(X(I)-X(K))**3 
CALL  SED I F4( TS , PROD , J , 1  ROW) 

572  IROW-IROW+1 

568  ISTRT-ISTRT+NJUL 

CALL  SEDIF5(TP,TS,FZD,PWR,J, ISTRT) 

573  IF( 1-1)575, 566, 57$ 

575  IROW-ISTRT+1 

DO  576  IXPWR-1,3 

PROD-PWR*X ( I ) ** I XPWR 

CALL  SEDIF5  (TP,TS,FZD,PROD,J, IROW) 

576  IROW-IROW+1 

IF( 1-3)566,566,577 

577  DO  578  K-3,1.2 

I F ( K— I ) 605,5o6,605 

605  PR0D-PWR*(X(I)-X(K))**3 

CALL  SEDIF5  ( TP, TS,FZD, PROD, J, IROW) 

578  IROW-IROW+1 

566  ISTRT-ISTRT+NJUL 

CALL  SEDIF6  (T,TP,TS,FZO, J.PWR, I STRT) 

579  IF( 1-1)581,580,581 

581  IROW-ISTRT+1 

DO  582  IXPWR-1,3 

PROD-PWR*X(l)**IXPWR 

CALL  SEDIF6  (T,TP,TS,FZD, J, PROD, IROW) 

582  IROW-IROW+1 

I F( 1-3)580,580,583 

583  DO  584  K-3,1.2 
IF(K-I)606,580,606 
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606  PR0D-PWR*(X( I )-X(K))**3 

CALL  SE0IF6  (T,TP,TS,FZD, J,PR0D, IROW) 

584  IROW-IROW+1 

580  I STRT- I STRT+NJUL 

CALL  SEDIF2  (FZD.T,TP,TS,PWR, ISTRT,J) 

585  IF( 1—1)586,590,506 

586  IROW-ISTRT+1 

DO  587  IXPWR-1,3 

PROD-PWR*X( I )**l XPWR 

CALL  SE0IF2  (FZD, T, TP, TS, PROD, IROW, J) 

587  IROW-IROW+1 

I F( 1-3)590,590,588 

588  DO  589  K-3,l,2 
IF(K-I)597,590,597 

597  PROD-PWR*(X( l)-X(K))**3 

CALL  SEDIF2  (FZD, T, TP, TS, PROD, IROW, J) 

589  IROW-IROW+1 

590  IF( IZ-JUL)591 ,1 11 *591 

591  I STRT- I STRT+NJUL 

DO  592  K-JUL, IZ,JUL2 
IF(K-IZ)593. 111,593 
593  D I F-FZD— Z(K) 

CALL  SE0IF2  (DIF.T,TP,TS,PWR, ISTRT, J) 
1  I F( 1-1)607.592,607 

607  IROW-ISTRT+1 

DO  608  I XPWR-1 ,3 

PROD-PWR*X( I  )**  I XPWR 

CALL  SEDIF2  (DIF, T, TP, TS, PROD, IROW, J) 

608  IROW-IROW+1 

I F( 1-3)592.592.598 

598  DO  599  KJ-3.1.2 
IF(KJ-I ) 600. 592, 600 

600  PR0D«PWR*(X(I)-X(KJ))**3 

CALL  SEDIF2  (DIF,T,TP,TS,PROD, IROW,J) 

599  IROW-IROW+1 

592  I STRT- I STRT+NJUL 
111  J-J+1 


110  NCQMP-1 
PUNCH  174 
PUNCH  175, COST 
DO  412  I -2, NOROW 
412  PUNCH  177, I ,RHS 
PUNCH  176 
PRINT  178 
PAUSE  751 
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171 

j73 

174 

175 

176 

B 

181 


GO  TO  750 

F0RHAT|6X^2H^C^|4,6H  0FSET.F12.5) 


FORMAT  ( 


FORMAT (7HF I RST  B) 

FORMATi  12X,  6H  R  1,F12.5) 
F0RMATC3HE0F) 

FORMATt 12X,2H  R,I4,1H  ,F11.5) 

FORMAT  28HENTER  NEW  DATA  -  PRESS  START) 
FORMAT  (10I5,2F15.9) 

FORMAT  (4F15.9) 

END 
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SUBROUT  I NE  PROFL 2  ( XD ,  FD ,  E  XPD  1 , TZPD , TZSD , E XPD2 ,  XI ,  TO ,  TPD ,  TSO ,  NCD ) 
TD«( X1+FD-XD)/X1 
IF(TD-. 00001)1, 1,2 

2  0IFD-EXPD1-1. 

PROO=TD**D I FD 
TPD1«EXP01*TZPD*PR0D 

TSi)1»£XPD1*(D  I  FD*TD**(D I  FD-t  ,)*TZPD**2+TZSD*PR0D) 

TD«1.-TD**EXPD1 
I Ff  TD-, 00001 ) 7, 7,4 

4  DIFD-EXPD2-1 . 

TDE:<P«TD**DiFD 
TPD— EXPD2*TPD1*TDEXP 

TSD«EXPD2*(D I FD*TD**(EXPD2-2.)*TPD1*TPD1-TSD1*TDEXP) 

T[)»*TO**E  XPD  2 
IF(TD-.00001)7,7,6 
7  TD-O. 

GO  TO  3 
6  NCD«*1 
RETURN 
1  NCD-2 
TD-1. 

5  TPO-O. 

TSD-O. 

RETURN 

3  NCO-1 
GO  TO  5 
END 


SUBROUTINE  SEDIF1  (TD,TPD,TSD,XD,JROWD, JD,PWRD, IROWD) 
IROWD-JROWD 

CALL  SEDIF4  (TSD.PWRD, JD, IROWD) 

IROWO-IROWD+1 

CALL  SEDIF5  (TPD, TSD, XD,PWRD,JD, IROWD) 

IR0WD=IR0WD+1 

CALL  *;E0IF6  (TO, TPD,TSD,XD, JD.PWRD,  IROWD) 

RETURN 

END 


SUBROUTINE  SEDIF2  (DIFD,TD,TPD,TSD,PWRD, IROWD, JD) 
WRM-( 3 .*D I FD*TD+3.*TPD*D I FD**2+(TSD*D I FD**3 ) / 2 . ) *PWRD 
CALL  SIGN  (JD, IROWD ,WRM) 

RETURN 

END 
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SUBROUTINE  SEDIF3  ( IXD.FXD, PRD, IROWD, JSD) 
EXO-IXO 

WRM-(EXD*FXD**( I XD-t )-FXD** I XD)*PRD 
CALL  SIGN  (JSD, IROWD, WRM) 

RETURN 

END 


SUBROUTINE  SEDIF4  (TSD,PWRD,JD, IROWD) 
WRM«(TSD*PWRD)/2. 

CALL  SIGN  (JD, IROWD, WRM) 

RETURN 

END 


SUBROUTINE  SEDIF5  (TPD,TSD,XD,PWRD,JD, IROWD) 

WRM-((TSD*XD)/2.+TPD)*PWRD 

CALL  SIGN  (JD, IROWD, WRM) 

RETURN 

END 


SUBROUTINE  SEDIF6  (TD,TPD,TSD,XD, JD,PWRD, IROWD) 
WRM-TD 

WRM-WRM+( TSD*XD**2 ) /2 . 

WRM-WRM+2 . *TPD*  XD 

WRM-WRM*PWRD 

CALL  SIGN  (JD, IROWD, WRM) 

RETURN 

END 


SUBROUTINE  SIGN  (JW,IW,WAM) 
IF(WAM)370,373,37t 
370  WAM— WAM 

PUNCH  169,JW,IW,WAM 
RETURN 


371  PUNCH  167,JW,IW,WAM 
373  RETURN 
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Section  III 


SMOG  3 


The  purpose  of  the  program  named  SMOG  3  is  to  generate  a  matrix  with 
the  following  characteristics: 

(1)  Primal  formulation 

(2)  Sum  of  the  deviations  fit 

(3)  Double  splined  in  the  x  direction 

(4)  Single  or  double  splined  in  the  z  direction 

(5)  Containing  profile  conditions 

Figure  B-4  (Section  II  of  this  Appendix)  shows  a  five-station  by  five- 
waterline  surface  and  an  equation  representing  this  surface.  The  matrix 
which  represents  this  equation  consists  of  two  parts.  The  first  part 
consists  of  the  deviation  constraints.  The  second  part  is  the  curvature 
constraints . 

The  deviation  constraints  have  the  form: 


(a) 

Y<V*J)  +  S1.J  -  V)  ■  Yi.j 

thus  offset 

(1,3)  Fig.  B-4  would  have  the  constraint 

(a) 

(a01  "  b01‘)x3  +  (a02  ‘  b02Jx32 

+  (a03  "  b03)x3  +  §1 

*01  •  b01  ■  *01 

a02  '  b02  "  A02 

a03  "  b03  “  A03 

1.0. 0-2 

B-193 

,3  ’ 


-  Y 


1,3 


The  curvature  constraints  are  the  same  as  those  of  SNOG  1  and  SMOG  2 


(b) 

f"(x. ,  z. ) 
v  i’  J  xx 

>  0 

(O 

f"  (x^  j  z  j  )zz 

^  o 

A  tableau  of  t’r>. 

matrix 

is  given  in  Fig.  B-6.  The  positive  identity 

matrix  correspo  iding  to  S  .  .  is  always  used  to  provide  a  basic  feasible 

*  >  J 

solution,  and  SMOG  3  provides  the  basis  headings  for  this  purpose. 


Data  Assembly  fnstiuctlons 


Before  the  data  punched  by  Passes  1  and  2  of  SMOG  3  can  be  used,  it 
must  be  rearranged  according  to  the  instructions  which  follow: 


Pass  1 

(1)  Be  sore  Pass  2  input  data  has  been  removed  from  behind  the 
output  deck 

(2)  Remove  the  first  four  cards  from  the  front  of  the  output 
decl  .  These  are  punched  with  the  heading,  ROW  ID,  OFSET, 
MATRIX. 

(3)  Sort  the  output  deck  on  Column  16  with  the  edit  off.  The 
output  deck  will  fall  in  the  reject  hopper  and  the  cost  row 
cards  will  fall  in  the  Column  2  hopper. 

Set  the  cost  row  elements  aside. 

Pass  2 

(1)  Remove  the  last  two  cards  from  the  output  deck.  These  are 
punched  FIRST  B  and  EOF 

(2)  Place  the  Pass  2  output  deck  directly  behind  the  Pass  1 
output  deck. 

Etff  l  JP*  1 

(1)  Sort  the  combined  deck  on  Column  8  with  the  edit  off.  The 
right-hand  eide  carde  and  row  idantification  will  fall  in 
the  reject  hopper.  Sat  thaae  aside. 

(2)  Sort  the  combined  deck  on  Column  23  with  the  edit  off.  The 
basis  heading  carda  will  fall  in  the  reject  hopper. 

(3)  Keypunch  a  card  with  BASIS  in  Columns  1*5  and  place  in  front 
of  the  basis  heeding  deck. 
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Fig.  B-6  Matrix  Tableau 


(4)  Sort  the  row  identification  -  right-hand  side  deck  on  Column  25 
with  the  edit  off.  The  row  identification  cards  will  fall  in 
the  reject  hopper. 

(5)  Place  the  following  cards  in  the  fron  of  a  card  box  in  the 
following  order: 

Heading  card 

Row  ID 

OFSET 

The  row  identification  deck 
The  basis  heading  deck 
MATRIX 

(6)  Sort  the  combined  output  decks  on  Columns  12,  11,  10.  This 
will  place  them  in  Column  order.  Put  the  resulting  deck  behind 
the  MATRIX  card. 

(7)  Put  the  Pirst  B  card  ahead  of  the  right-hand  side  deck  and 
the  EOF  card  behind  it.  Then  put  this  deck  behind  the  combined 
output  decks. 


A,  fully  assembled  output  deck  combining  the  sample  problem  of  Pfas  1  and 
Pass  2  is  given  immediately  following  the  Pass  2  sample  output  listing. 
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SMOG  3  PASS  1 


OPERATING  INSTRUCTIONS 

Pass  1  of  SMOG  3  calculates  the  second  differences  of  the  offsets 
and  the  Intersections  of  the  profile  with  the  water lines  and  stations. 
It  then  punches  the  deviation  constraints  for  the  matrix.  The 
program  is  written  in  FORTRAN  II  for  the  IBM-1620  computer. 


Fortran  Symbol  Definitions 


M^' 

N 

X(J) 

Z(I) 

Y&J) 

IAY 

IBY 

ICY 

NAA 


SX 

TX 

sz 

TZ 


Total  number  of  waterlines  in  the  surface  to  be 
faired.  This  will  be  an  odd  number  if  the  surface 
is  double  splined  in  the  z  direction. 

Total  number  of  stations  in  the  surface.  This  will 
always  be  an  odd  number. 

Distance  along  the  x  axis  from  the  origin  of  the 
surface  to  station  j  . 

Distance  along  the  z  axis  from  the  origin  of 
the  surface  to  Hater line  i  . 

Offset  of  the  surface  at  waterline  i  ,  station  j  , 
in  feet  and  decimals 

Number  of  feet  in  the  offset  if  read  in  feet,  inches 
and  eighths. 

Number  of  Inches  corresponding  to  IAY 
Number  of  eighths  corresponding  to  IAY 
Program  indicator 

If  NAA  ■  -1,  scale  factors  for  normalizing  and 

translating  x(J)  and  z(I)  are  read  in 

If  NAA  “  +1,  scaling  and  translating  are  done  by 
the  program  as  given  in  Section  II  . 

Scale  factor  for  X(J)  when  these  are  read  in 

Translation  factor  for  X(J)  when  these  are  read  in 

Scale  factor  for  Z(I)  when  these  are  read  in 

Translation  factor  for  Z(I)  when  these  are  read  in 
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RHS 

CURV 

NAB 


COST 

EXP1 

EXF2 

POINT 


C(I) 


The  cosC  on  the  slack  variables  associated  with 
the  deviation  constraints.  Usually  made  equal  to 
1000  . 

The  cost  on  the  non-slack  variables.  Usually  made 

.001  . 

Program  Indicator 

If  NAB  ■  -1  ,  the  surface  will  be  double  splined 
In  the  z  direction,  and  M  must 
be  an  odd  number 

If  NAB  ■  +1  ,  the  surface  will  be  single  splined 
In  the  z  direction  and  M  is 
unrestricted 

Used  in  Pass  2.  A  zero  value  must  be  given 

The  value  of  p  In  the  profile  function,  usually 
made  .3333 

The  value  of  r  In  the  profile  function,  usually 
3.0  . 

The  value  of  D  in  the  profile  function.  This 
nuaber  is  the  length  of  the  fraction  of  the  station 
spacing  over  which  the  profile  function  is  effective 
in  terms  of  the  full-scale  surface  (the  value  will 
be  translated  and  normalized  exactly  like  the  x 
values  are). 

Coefficients  of  the  Theilheimer  equation  describing 
the  profile.  There  must  be  (Mt2)  coefficients. 
Straight  lines  or  other  curves  may  be  used  by  making 
all  except  the  desired  coefficients  zero. 


The  input  data  cards  are  described  below.  There  are  certain  input 
variables  on  which  there  are  restrictions  on  the  value  which  can  be 
entered.  These  restrictions  follow. 

Variable  Minimum  Maximum 

N  3  15 

M  3  11 

-  In  describing  the  various  data  cards,  the  actual  FORTRAN  format 
field  description  is  used  in  most  cases.  These  fields  come  con¬ 
secutively  across  the  card  with  no  gaps  or  blank  columns  between. 


( 


except  where  lndiceted.  The  field  deecriptlone  ere  the  FORTRAN  V 
field,  which  ueee  the  FORTRAN  fixed  point  decimal  number,  end  the 
I  field  that  usee  the  FORTRAN  integer  number,  which  is  elweya 
right  justified.  The  card  numbers  ere  not  punched  on  the  date 
cards . 


Contents  of  Card 


Card  So. 


Card  Columns  1  2-12  13-53 

Contents  *  Blank  Any  alphameric  problem 

heading  (This  heading 
will  be  used  for  identify- 


ing  the  problem 
output) 

on 

the  IP 

Format 

15 

15  F11.6 

F11.6 

2 

Variable 

NAA 

NAB  POINT 

COST 

Format 

15 

15  F11.6 

F11.6  F11.6 

F11.6 

3 

Variable 

M 

N  CURV 

RHS  BXF1 

XXP2 

Format 

F10.4 

Next 

Variable 

Format 

C(l) 

F10.4 

M+2  cards 

Variable 

Z(I) 

Alternate  Cards: 

Naxt 

If 

Format 

F10.4 

F10.4 

1 

Variable 

X(J) 

or 

U  1.J) 

N  cards 

H 

Groups 

per 

group 

Format 

F10.4 

15  15 

15 

Variable 

X(J) 

IAY  IBY 

ICY 

This  arrangement  can  be  stated  as  follows:  The  first  card  of  each 
set  gives  the  height  z  of  the  waterline.  The  following  cards 
each  give  a  station  coordinate  and  the  offset  at  the  intersection  of 
the  given  waterline  and  station.  There  is  one  set  of  cards  for 
each  waterline  in  the  surface. 


Format  F10.4  F10.4  F10.4  F10.4  Last  Card 

Variable  SX  SZ  TX  TZ 

(This  card  is  entered  only  if  NAA  ■  -1) 


( 
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fia&Bas. 


All  of  the  date  necessary  for  the  deviation  constreinte  In  the 
matrix  ere  punched  by  this  program.  The  first  cards  punched 
contain  the  heeding,  "ROW  ID,"  "OFSET"  (this  is  the  identification 
of  the  cost  row)  end  "MATRIX." 


The  matrix  elements  are  punched  in  row  order  with  the  columns  in 
order  within  the  rows.  Bach  row  also  contains  in  order  of 
appearance : 

(1)  The  row  identification  for  that  row 

(2)  The  right-hand  side  element  for  that  row 

(3)  The  costs  of  the  two  slack  columns  associated  with 
the  row 

(4)  The  slack  elements 

(5)  The  basis  heading  card 

This  output  must  be  sorted,  together  with  the  output  from  Pass  2, 
before  being  used  with  IP-90  .  The  sorting  instructions  are  in 
the  section  called  "Output  Assembly." 


Sense  Switch  Settings 


Sense  9witch  1 
Sense  9wltch  2 

Sense  Switch  2 
Sense  Switch  3 

Sense  9witch  3 
Sense  Switch  4 
Sense  Switch  4 


Not  Used 

ON  The  values  of  the  x  and  z  second 

differences  and  the  resulting  signs  are 
printed 

OFF  The  above  is  not  printed 

ON  The  offsets  are  read  in  as  feet,  inches 
and  eighths 

OFF  The  offsets  are  read  in  as  feet  and  decimals 

ON  The  data  for  Pass  2  is  punched 

OFF  The  data  for  Pass  2  is  not  punched  but 
remains  in  storage  for  later  use  with 
Pass  2  . 


SAM?  12  PB0B1XM 


The  following  sample  problem  Is  a  five-station  by  five-waterline 
set  of  offsets  taken  from  the  D LG-26  class  ships.  The  statione  are 
Numbers  0  ,  1/2  ,  1  ,  1-1/2  ,  2,  and  the  waterlines  are  12,  16,  20, 

24,  end  28.  The  surface  contains  the  bow  profile  of  the  DIG.  The 
matrix  is  double  splined  in  the  »  direction.  Only  the  output  from 
Peas  1  in  the  order  in  which  it  was  punched  is  listed  here.  The  output 
from  Pass  2  is  included  with  the  instructions  from  that  program.  The 
sorted,  assembled,  and  combined  outputs  of  Pass  1  and  2  is  listed  after 
the  Pass  2  output. 
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SAMPLE  PROBLEM  FOR  SMOG  3,  PASS  1 


*  *  SAMPLE  INPUT  *  * 


* 


I 


TODD 


1  -1 
5  5 

1.52935 
■  .5873 


0. 

0. 

0. 

0. 

0. 


.0000 

.0000 

12.7500 

25.5000 

38.2500 

51.0000 

1.0000 

.0000 

12.7500 

25.5000 

38.2500 

51.0000 

2.0000 

.0000 

12.7500 

25.5000 

38.2500 

51.0000 

3.0000 

.0000 

12.7500 

25.5000 

38.2500 

51.0000 

4.0000 

.0000 

12.7500 

25.5000 

38.2500 

51.0000 


OLG  5X5  END  PROFILE  X-Z  DBL  SPL-1 2/3/62 

1.  0. 

0.001  1000.  3.  .33333333 


.000000 

1.041666 

2.458333 

3.874999 

5.333333 


.000000 

1.239583 

2.697916 

4.197916 

5.812500 

.078  _ 
1.479166 
2.979166 
4.583333 
6.333333 


.343749 

1.791666 

3.447916 

5.250000 

7.166666 


.708333 

2.416666 
4.322916 

6.41 6666 
8.583333 


TOOD  DLG  5X5  END  PROFILE  X-Z  DBL  SPL-1 2/3/62 


* 


END  PROFILE  INTERSECTIONS 
WLINES 

1  1.529350 

-.5873  0.0000 

2  .942050 

-.5873  0,0000 

3  .354750 

-.5873  0.0000 

4  -.232550 

-.5873  0.0000 

5  -.819850 

-.5873  0.0000 

STATIONS 

1  2. 604)03 

2  .901367 

3  0.000000 

4  0.000000 

5  0.000000 


2ND  DIFF  -  WLINES 
WL  STA  NORM  VALUE  SIGN 


2 

2 

2 

2 

2 

3 

3 

3 

3 


4 

4 

4 

4 


1 

2 

3 

4 

5 
1 
2 


5 

1 

2 

3 

4 

5 
1 
2 


1 

5 


OUT  OF  SURFACE 
OUT  OF  SURFACE 
-5.176766  1.00 

.041668  -1.00 

.041668  -1.00 

OUT  OF  SURFACE 
-37.680848  1.00 

.041667  -1.00 

.114584  -1.00 

.114584  -1.00 

OUT  OF  SURFACE 
-.963248  1.00 

.104167  -1.00 

.145833  -1.00 

.145839  -1.00 

-.049093  1.00 

.208333  -1.00 

.145834  -1.00 

.114582  -1,00 

.114582  -1.00 
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uurvimmn 


o 


t  .927938  -1.00 

2  .197917  -1.00 

3  .187500  -1.00 

4  .072917  -1.00 

5  .072917  -1.00 


2ND  DIFF-STATIONS 
STA  WL  NORM  VALUE  SIGN 


1  2  OUT  OF  SURFACE 

1  3  OUT  OF  SURFACE 


1 

4 

-.439121 

1.00 

1 

5 

-.439121 

1.00 

2 

2 

-3.216759 

1.00 

2 

3 

.072917 

-1.00 

2 

4 

.312500 

-1.00 

2 

5 

.312500 

-1.00 

3 

2 

.041667 

-1.00 

3 

3 

.187500 

-1.00 

3 

4 

.406250 

-1.00 

3 

5 

.406250 

-1.00 

4 

2 

.062500 

-1.00 

4 

3 

.281250 

-1.00 

4 

4 

.499999 

-1.00 

4 

5 

.499999 

-1.00 

5 

2 

.041666 

-1.00 

5 

3 

.312500 

-1.00 

5 

4 

.583334 

-1.00 

5 

5 

.583334 

-1.00 
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FLOW  DIAGRAM  -  SMOG  3,  PASS  1  (See  SMOG  2,  Past  1  for  flow  diagrams 

of  tha  subroutines) 
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DO  HOT  MODIFY 
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DO  M 

IT.  3,1)1 


1*  >  IT 


PUNCH  ROW  III  Ml  NT! 
OF  TH1  FORM  ' 

TRH.(l(t*)-  *(?T))** 
ZTKM 
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FORTRAN  LISTING  -  SMOG  3,  PASS  1 


c 

) 

), 


SMOG3  PASS  1  PRIMAL  SUM  OF  DEV  X  OR  X  AND  Z  DDL  SPL  Wl  PRU 
DIMENSION  RX(11,15),RZM5,11),TRM(150),  IRW(250)  ,FZ(  11 )  ,FX(  15 

DIMENSION  X(15),Z(11),Y(11,15),WL(15,11),IWL(15),C(13),FPZ(11 

CFSZ(II) 

COMMON  JPLS,JUL,JUL2,NJUL,M,N,MS,NS,  NUCOL , NCOMP , CURV , COST 
COMMON  EXP2,COST,PO INT.EXP1 , IWW, ICOL 
COMMON  X,FZ,RZ,Z,FX,FPZ,FSZ,RX,C 

*  *  INPUT  *  * 

750  READ  173 
PRINT  173 

READ  150,NAA,NAB,POINT,COST 
READ  150, M,N, CURV, RHS,EXP1,EXP2 
ML-M+2 

DO  751  1-1  .ML 

751  READ  154, C*  I) 

NCOMP- 1 

NL-N+1 

ML-M+1 

NS-N-1 

MS-M-1 

DO  2  1-1 .M 
READ  154, Z(l) 

DO  2  J-1,N 

IF(SENSE  SWITCH  3)802, 803 

802  READ  183,X( J) , I  AY, IBY, ICY 
AAY-IAY 

BBY-IBY 

CCY-ICY 

Y( I ,J)«AAY+(BBY+(CCY/8.))/12, 

GO  TO  2 

803  READ  154,X(J),Y(I,J) 

2  CONTINUE 
!F(NAA)3,3,4 

3  READ  154.SX,SZ,TX,TZ 
GO  TO  508 

4  SX-X(2)-X( 1 ) 

SZ-Z(2)-Z(1) 

TZ-Z(I) 

TX-X(1) 

★  *  NORMALIZING  *  * 

508  DO  26  J-1.N 

26  X(J)-(X(J)-TX)/SX 
DO  27  1-1, M 

27  Z( I )-(Z( I )-TZ)/SZ 
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8  ) 


691 


700 

690 


693 

694 


204 

510 


695 

71 

65 

66 

S 


67 

68 
696 
697 


PO I NT-(PO I NT-TX)/SX 

★  *  END  PROFILE  INTERSECTIONS  *  * 

I F( SENSE  SWITCH  2)509.820 
PRINT  165 
DO  694  1-1, M 

FX( I )-C(1)+C(2)*Z( I )+C(3)*Z( I )*Z( I )+C(4)*Z( I )**3 
CALL  DERI VI  (C(2),C(3)»C(4),Z( I ),FPZ( I ) , FSZ ( I )) 

I F( 1-3)690,691,691 
DO  700  J-3.1 
DIF-Z(I)-Z(J-1) 

FX( I )-FX( I )+C( J+2)*D I F**3 

CALL  DERIV2  (C( J+2) ,DI F,FPZ( I ) ,FSZ( I ) ,FPZ( I ) ,FSZ( I ) ) 

I F( 1-3)690,690, 700 

CONTINUE 

FPZ( I }-FPZ( I )/SZ 
FSZ(I)-FSZ(I)/(SZ*SZ) 

I F( SENSE  SWITCH  2)693.694 
PRINT  170, I .FX( I) 

PRINT  154,FPZ( I ) , FSZ ( I ) 

CONTINUE 

J-1 

OLD  -Z(M) 

I -MS 

I F(SENSE  SWITCH  2)204,510 

PRINT  179 

FZ(J)-FX(I) 

ZZ-Z(I) 

I  COMP- 1 

IF(ABSF(X( J)-FZ(J) )-,0001 )  66,66,71 
GO  TO  (65,69), ICOMP 
I F(X( J)-FZ(J) 5  67,66,72 
FZ(J)-ZZ 

I F(SENSE  SWITCH  2)  74,80 
PRINT  170,J,FZ(J) 

J-J+1 

OLD-ZZ 

GO  TO  510 

ZZ-ZZ+(0L0-ZZ)/2. 

I  COMP-2 
I JUL-l+1 

FZ(J)-C(1)+C(2)*ZZ+C(3)*ZZ*ZZ+C(4)*ZZlWr3 
IF( I JUL-3)  695,696,697 
FZ(J)-FZ(J)+C(5)*(ZZ-Z(2))**3 
GO  TO  695 
DO  699  K-3.IJUL 


1.0. 0-2 
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c 


c 


699  FZ(J)-FZ(J)+C(K+2)*(ZZ-Z(K-1))**3 
60  TO  695 

69  IF(X(J)-FZ(J))  67,66,70 

70  SE6-0LD-ZZ 
OLD-ZZ 
ZZ-ZZ-SEG/2. 

GO  TO  68 

72  OLD«Z(!) 

1-1-1 


I F( 1-1 >61,510,510 

81  DO  82  l-J.N 

FZ(l)-Z(lJ  . 

I F(SENSE  SWITCH  2)83,82 
83  PRINT  170,I,FZ(I) 

82  CONTINUE 

*  *  WATERLINE  2ND  DIFFERENCES  * 

805  IF(SENSE  SWITCH  2)6,7 

6  PRINT  157 

7  DO  15  1-1 ,M 
DO  13  J-1.NS 

GO  TO  (8,513), NCOMP 

8  ZY-O, 

I F(X\J)-FX( I )) 10,10,9 
10  R-20. 

I F(SENSE  SWITCH  2)660,13 
660  PRINT  151, 1, J 
GO  TO  13 

9  ZZ-FX(I) 

NCOMP-2 
GO  TO  1 

513  ZZ-X(J-I)  % 

jymyf |  jM | j 

1  CALL  d!f2ND  (X(J),X(J+1),ZZ,Y(I,J),Y(I,J+1),ZY,I, 

13  RX( I ,J j-R 

RX( I ,N)-RX( I ,NS) 

I F ( SENSE  SWITCH  2)14,15 

14  PRINT  158, I ,J,A,RX(I ,N) 

15  NCOMP-1 

*  *  STATION  2ND  DIFFERENCES  *  * 

I F(SENSE  SWITCH  2)16,17 

16  PRINT  159 

17  DO  25  J-1 ,N 
DO  23  1-2, MS 

GO  TO  (18.19),NC0HP 

18  IF(Z(I)-FZ(J))21,21,20 


J,R,A) 
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21  R-20. 

I F (SENSE  SWITCH  2)22,23 

22  PRINT  151, J, I 
GO  TO  23 

20  ZZ-FZ(J) 

NCOMP-2 
GO  TO  514 

514  call*d7f2no  (Z(I),Z(I+1),ZZ,Y(I,J),Y(I+1,J),Y(I-1,J),J,I,R,A) 

23  RZ(J,I)-R 
(J.H)-F 


RZ 

IF 


j,n/>RZ(J,HS) 

SENSE  SWITCH  2)24,2' 


_ >,25 

24  PRINT  158,J,I,A,RZ(J,H) 

25  NCOMP-1 

*  * 

PUNCH  173 

NJUL-(((N-3)/2)+4)*2 

NP^N 

IF(NAB)719, 718,718 

719  JUL-3 
JUL2-2 

NUC0L-((N+5)*(Hf5))/2 
GO  TO  720 
718  JUL-2 

JUL2-1  ,  % 

NUC0L-2*(M*2)*(((N-3)/2)+4) 

720  PUNCH  160 

I COL-NUCOL-1 

*  ★ 

PUNCH  164 

* 


ROW  ID 


MATRIX 


FIRST  WATERLINE 


MLG-NN 

1-1 

FUPZ— 1. 

FUSZ-O. 

DO  37  J-1.MLG' 

810  IF(RX(1,I)-20«)250,251,250 
251  K-2 

GO  TO  252 
250  K-1 

200  GO  TO  (201 ,702).NC0MP 

201  CALL  PROFLl(X( |J,FX(1),EXP1 ,FUPZ,FUSZ,EXP2,P0INT,T,TP,TS 
,NCOMP) 
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( 


257 

703 

702 

90 

252 


253 

31 

32 


133 

33 


254 

255 

34 

35 
435 


463 


36 

37 


I F(T-.0001 >703,703,257 

ZY-Y(1.I)/T 

GO  TO  90 

NCOMP-2 

ZY-Y(I.I) 

PUNCH  162, J 
PUNCH  177,J,ZY,J,J 
I  ROW-3 

00  31  NPR-1,3 

KR0W-IR0W+1 

TRM( I R0W)-X( I )**NPR 

TRM(KROW)-TRM(IROW) 

GO  TO  <253.31),K 

PUNCH  167, IROW, J,TRM( IROW), KROW,. 

IROW-IROW+2 

IF( 1-3)35.35,32 

00  34  IT-3.1 

KROW-IROW+1 

IF(l-IT)34,34f133 

IF(lT-2*(IT/2))33,34,33 

TRM(IR0W)-(X(I)-XUT))**3 

TRM(KROW)-TRM( I  ROW) 

GO  TO  (254. 255), K 

PUNCH  167, IROW,J,TRM( IROW) ,KROW,« 

IROW-IROW+2 

CONTINUE 

GO  T0(435,463),K 
I COL- I COL+2 
KCOL-I C0L+1 

PUNCH  171, ICOl,RHS.KCOL,RHS 

PUNCH  152, ICOL,J,K*OL,J, ICOL,J 

JPLS-J+1 

I  ROW- 1  ROW-1 

IWL( I )-IROW 

DO  36  JR-3, IROW 

WL( I . JR)-TRM(JR) 

1-1+1 

NCOMP-1 

K3-1 

*  *  REMAINING  WL  * 

DO  64  IZ-2,M 
JL-J 

M6-JL+NM-1 
I X— 1 

DO  63  J-JL.MG 


,TRM( IROW) 


,TRM( IROW) 
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IF-2 

I RW( I F)» I ROW+1  ' 

812  IF(RX(IZ,IX)-20.)202,  63,202 
202  GO  TO  (91,92) .NCOMP 

91  CALL  PROFL2(X(IX),FX(IZ),EXP1,FUPZ,FUSZ,EXP2,POINT,T,TP,T$ 
, NCOMP) 

IF(T-.0001)93,93,258 
258  ZY-Y(IZ,IX)/T 
GO  TO  95 
93  NCOMP-2 

92  ZY-Y(IZ.IX) 

95  PUNCH  1b2,J 

PUNCH  177.J,ZY,J#J 
708  IF(J-JL)41,41,44 

41  DO  42  IZPR-1,3 
TRM(IF)-Z(IZ)**IZPR 
KR0W-IRW(IF)+1 

PUNCH  167,IRW(IF),J,TRM(IF),KR0W,J,TRM(IF) 

IF-IF+1 

42  IRW(IF)-IRW(IF-1)+NJUL 
IF< I Z-JUL) 146, 146,142 

142  DO  145  ITZ-JUL.IZ 
IF(IZ-ITZ)146, 146,729 

729  GO  TO  ( 1 4h, 1 44, 730) , JUL 

730  IF(ITZ-2*(ITZ/2))144,145,144 

144  ZTR-(Z(IZ)-Z(ITZ))**3 
TRM(IF)-ZTR 

KR0W-IRW(IF)+1  ,  . 

PUNCH  167,IRW(IF),J,TRM(IF),KR0W,J,TRM(IF) 

IF-IF+1 

IRW(IF)-IRW(IF-1)+NJUL 

145  CONTINUE 

146  JPLS-J+1 
GO  TO  259 

C  *  *  Z  .  Z  SQR  ,  Z  CUBE  ROWS  *  * 

44  IRO-1 
K-IWL(IX) 

DO  244  I J-3,K,2 
IRO  -IRO  +2 
KR0W-IR0+1 

244  PUNCH  167,IR0,J,WI.(IX,U),KR0W,J,WL(IX,IJ) 

DO  50  IZPR-1.3 
ISTRT-IRW(IF)+NJUL 
ZTRM-Z(IZ)**IZPR 
TRM(IF)-ZTRM 
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45 

46 

S 


49 

50 

53 

54 

731 

55 

56 


ITj)**3)*ZTRM 


256 

!  I 


KR0W-IRW(IF)+1 

PUNCH  167, IRW(IF) ,J,TRM( I F) ,KR0W,  J 
IF-IF+1 

IRW(IF)-IRW(IF-1)+2 

00  45  NPR-1,3 

TRH( I F)-{ X( I X)**NPR)*ZTRM 

KR0W-IRW(IF)+1 

PUNCH  1 67, I RW( I F) , J , TRM( I F) ,KR0W, J 
IF-IF+1 

IRW(IF)-IRW(IF-1)+2 
I F( I X— 3)50,50,46 
DO  49  IT-3, IX 
IF(IX-IT)49,49.47 

IF(lT-2*(IT/2))48,49.48 

TRM(IF)-((X(IX)-X(IT] 

KROW-IRW( IFj+1 
PUNCH  167, IRW( I F) , J,TRH( IF) ,KROW,J 
IF-IF+1 

IRW( IF)-IRW( IF— 1)+2 

CONTINUE 

IRW( IFJ-ISTRT 

*  *  THE  I LHE I MER  ROWS 

IF( IZ-JUL)61 .61,54 
DO  60  ITZ-JUL.IZ 
GO  TO  (55,55, 731 ),JUL 
IF(ITZ-2*(ItS/2))55,60,55 
ISTRT-IRW(IF)+NJUL 
IF(IZ-ITZ)61 ,61,56 
ZTRM-(Z(IZ)-Z(ITZ))**3 
TRH(IF)-ZTRH 
KROW-IRW( IF)+1 

PUNCH  167, I RW( I F) , J ,TRM( I F) ,KROW,. 
IF-IF+1 

IRW(IF)-IRW(IF-1)+2 

DO  256  NPR-1,3 

TRM( I F)-( X( I X)**NPR)*ZTRM 

KROW-IRW( I F)+1 

PUNCH  1 67, I RW( I F) , J , TRM( I F) ,KROW,, 
IF-IF+1 

IRW(IF)-IRW(IF-1)+2 
DO  59  IT-3. IX 
IF(IX-IT)59.59.57 
I Fl IT-2*( IT/2) 158,59,58  % 

TRM( I F)-( ( X( I X)-X( I T) )**3)*ZTRM 
KROW-IRWC I F)+1 

PUNCH  167, I RW( I F) ,J,TRM( I F) ,KROW,, 


,TRM( IF) 


,TRH( IF) 


,TRM( IF) 

*  * 


,TRH( IF) 


,TRM(IF) 


,TRH( IF) 
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IF-IF+1 

IRW(IF)-IRW(IF-1)+2 

59  CONTINUE 
IWW-KROW 

60  IRW(IF)-ISTRT 

61  JPLS-J+1 

GO  TO  (259.63) ,K3 
259  ICOL-ICOL+2 
KC0L-IC0L+1 

PUNCH  171,IC0L,RHS.KC0L,RHS 
PUNCH  152,ICOL,J,KCOL,J,ICOL,J 

63  IX-IX+1 

64  NCOMP-1 

DO  62  1-1, I WW 

62  PUNCH  171.1 .CURV 

I F(SENSE  SWITCH  3)400,401 

400  CONTINUE 

PUNCH  180,JPLS,JUL,JUL2,NJUL,H,N,HS,NS,mdcol  .NCOMP, 
PUNCH  1 81 ,EXP2,C0ST,P0INT,EXPt, IWW, ICOL 
DO  87  1-1 In 
PUNCH  181,X(I),FZ(I) 

DO  87  J-2,M 
87  PUNCH  181 ,RZ( I ,J) 

DO  88  1-1 ,M 

PUNCH  181 ,Z( I ) ,FX( I ) ,FPZ( I ) ,FSZ( I ) 

DO  88  J-1 ,N 
88  PUNCH  181,RX(I,J) 

ML-H+2 

DO  100  J-1.ML 
100  PUNCH  181,C(J) 

401  PRINT  178 
PAUSE  751 
GO  TO  750 

150  F0RMAT(2I5,4F1 1.6) 

151  F0RMAT(2I5,17H  OUT  OF  SURFACE) 

152  F0RMAT(7X,1HC,I4,2H  R,I4,8H  1 ,0/7X,1HC, I4,2H  R, 

C  7X.1HC, I4,2H  R, 14) 

154  FORMAT (4r 10,4) 

157  F0RMAT(//10X,17H2ND  DIFF  -  WLINES// 

C30H  WL  STA  NORM  VALUE  SIGN/) 

158  F0RMAT(2I5.F12,6.F8.2) 

159  F0RMAT(//10X,17H2MD  DIFF-STATIONS// 

C30H  STA  WL  NORM  VALUE  SIGN/) 

160  FORMAT (6HR0W  ID/12X,6H  OFSET) 


COST, CO ST 


14, 8H- 
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162 

164 

165 
167 

169 

170 

171 

172 

173 

177 

( 

178 

18 

181 

182 

"ill 


F0RMAT(1 2X,2H  R.I4) 

FORMAT  6HHATR IX) 

FORMAT  //6X.25HEND  PROFILE  I NTERSECT I ONS// 1 5  X , 6HWL I NES/ ) 
FORMAT  6X,2H  C,I4,2H  R, I4,1H-,F11.5/6X,2H  C,I4,2M  R,I4,F12.5) 
FORMAT  6 X,2H  C.I4,2H  R, 14, FI 2.5) 

FORMAT  5X, 1 5, FI 2.6) 

FORMAT  6X.2H  C,l4,6H  0FSET,1H  ,F11.5) 

FORMAT  215) 

FORMAT \1H*, 5 X,47H 

F0RMAT(12X,2H  R,I4.1H  ,F1 1 ,5/7X,6HC  1  ,1HR, 14, 1H-,4X,3H1 .0/ 

:  7X,6HC  2  ,1HR.I4,5X,3H1.0) 

FORMAT ( 28HENTER  NEW  DATA  -  PRESS  START) 

FORMAT 
FORMAT 
FORMAT 
FORMAT ( 


/14XV8HSTATI0NS/) 

10I5',2F15.9) 

4F15.9,I4,I4) 

X,2H  C,I4,6H  0FSET,1H  ,F11.5/6Xf2H  C,I4,6H  0FSET.1H  F 


F0RMAT(F10.5,3I5) 

END 


Mote:  See  SMOG  2,  PASS  1  for  listings  of  the  subroutines 
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SMOG  3  PASS  2 


OPERATING  INSTRUCTIONS 


Pass  2  of  the  program  accepts  the  data  punched  in  Pass  1.  It  produces 
the  curvature  constraints  for  the  matrix. 


Fortran  Input  Symbols 


JPLS 

JUL 


JUL2 


NJUL 


M 

MS 

N 

NS 

NUCOL 


NCOMP 

COST 


The  number  of  the  last  column  produced  in  Pass  1. 
Program  indicator 

JUL  ■  3  if  the  surface  is  double  splined 
in  the  z  direction 

JUL  ■  2  if  the  surface  is  single  splined 
in  the  z  direction 

This  value  is  equal  to  the  subscript  of  the  first 
z  value  of  the  point  of  discontinuity. 

Program  indicator 

JUI2  ■  2  if  the  z  direction  is  double  splined 

JU12  *  1  if  the  z  direction  is  single  splined. 

This  value  is  equal  to  the  number  of  waterline 

intervals  between  points  of  discontinuity. 

N“3 

This  value  2  (^-1  4)  is  equal  to  the  maximum 

number  of  rows  necessary  to  store  the  x  terms, 
each  of  which  is  multiplied  by  a  specific  z  term 

Total  number  of  water lines  in  the  surface 

M-l 

Total  number  of  stations  in  the  surface 
N-l 

Total  number  of  columns  in  the  non-slack  matrix 
For  z  double  spline,  NUCOL  - 

For  z  single  spline,  NUCOL  “  2(M+2)NJUL 
Program  Indicator 

Becomes  equal  to  curvature.  Not  used  (See  RHS 
and  CURV) 


1.0. 0-2 


B-257 


EXP1 

EXP2 

POINT 


C(I) 

CURV 

RHS 

X(J) 

Z(I) 

FX(I) 

FZ(J) 

FPZ(I) 

FSZ(I) 

RZ(X,J) 

KX(I,J) 

ICOL 

IWW 


The  value  of  p  In  the  profile  functions,  usually 
made  .3333 

The  value  of  r  in  the  profile  function,  usually 
made  3.0  . 

The  value  of  D  in  the  profile  function.  This 
number  is  the  length  of  the  fration  of  the  station 
spacing  over  which  the  profile  function  is  effective 
in  terms  of  the  full-scale  surface  (the  value  has 
been  translated  and  normalized  exactly  like  the 
x  values  are). 

Coefficients  of  the  Theilheimer  equation  describing 
the  profile.  There  must  be  (W-2)  coefficients. 
Straight  lines  or  other  curves  may  be  used  by 
making  all  except  the  desired  coefficients  zero. 

Is  equal  to  the  value  of  COST  in  Pass  1.  This  is 
the  value  of  the  cost  on  the  slacks  corresponding 
to  the  curvature  constraints.  Usually  equal  to 
zero . 

Equal  to  the  value  of  COST  from  Pass  1  .  This  is 
the  value  of  the  right-hand  side  for  the  curvature 
constraints. 

Distance  along  the  x  axis  from  the  origin  of 
the  surface  to  Station  j  . 

Distance  along  the  z  axis  from  the  origin  of 
the  surface  to  Waterline  i  . 

The  x  coordinate  of  the  intersection  of  Waterline 
i  with  the  profile. 

The  z  coordinate  of  the  intersection  of  Station 
j  with  the  profile 

The  first  derivative  with  respect  to  z  of  the 
profile  equation  at  Waterline  i  . 

The  second  derivative  with  respect  to  z  of  the 
profile  equation  at  Waterline  i  . 

The  sign  of  the  second  difference  of  the  Station 
Curve  j  at  Waterline  i  .  This  will  have  the 
opposite  sign  from  the  actual  second  difference. 

The  sign  of  the  second  difference  of  Waterline 
Curve  i  at  Station  j  .  This  will  be  the 
opposite  of  the  actual  sign 

The  column  number  of  the  first  curvature  slack. 

Not  used  in  Pass  2  . 
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The  following  input  formats  are  listed  as  a  matter  of  record  in  case 
modifications  of  the  data  cards  are  to  be  made.  Normally,  this  data 
will  be  obtained  from  executing  Pass  1  in  a  form  which  can  be 
entered  directly  into  Pass  2. 

Contents  of  Card  Card  NO. 


Format 

15  15  15 

15 

15 

15 

15  15 

Variable  JPLS  JUL  JUL2 

NJUL 

M 

N 

MS  NS 

Format  (Continuation  of 

15 

15 

F15.5 

F15.5  1 

Card  1) 

Variable 

NUCOL 

NCOMP 

RHS 

CURV 

Format 

F15.9 

F15.9 

F15.9 

F15.9 

14 

14  2 

Variable 

EXP2 

COST 

POINT 

EXP1 

IWW 

ICOL 

Format 

F15.9 

F15.9 

Variable 

X(I) 

FZ(I) 

1st  card/set 

Format 

F15.9 

Next 

w 

Variable 

Format 

Variable 

R2(I,J) 

F15.9 

Z(I) 

F15.9 

FX(I) 

F15.9 

FPZ(I) 

F15.9 

FSZ(I) 

M-l  Cards /Set 

Sets 

1st  Card/Set 

Next 

Format 

F15.9 

M 

Variable 

RX(I, J) 

N-l  Cards/Set 

Sets 

Format 

F15.9 

Next 

Variable 

C(J) 

M+2 

Cards 

Output 

This  program  produces  the  curvature  constraints  in  row  order  with  the 
columns  ordered  within  the  rows.  It  also  produces  the  row  identifications, 
slack  cost  entries,  slack  elements,  basis  headings,  and  right-hand  side 
entries  for  these  rows.  The  output  must  be  sorted  together  with  the 
output  from  Pass  1  before  use. 


Sense  Switches : 

Sense  Switch  3  ON  -  The  data  punched  in  Pass  1  is  read 

Sense  Switch  3  OFF-  The  data  is  picked  up  from  the  common  storage 
locations . 


SAMPLE  PROBLEM  -  SMOG  3,  PASS  2 


*  *  SAMPLE  INPUT  *  * 


26  3  2 

1000,00000000 
.33333333 
91 

0.00000000 
20.00000000 
20.00000000 
1 .00000000 
1 .00000000 
1 .00000000 
1 .00000000 
-1 .00000000 
-1 .00000000 
-1 .00000000 
2.00000000 
-1 .00000000 
-1 .00000000 
-1 .00000000 
-i.nonooooo 
3.0 0000000 
-1.00000000 
-1.00000000 
-1 .00000000 
-1.00000000 
4.00000000 
-1 .00000000 
-1.00000000 
-1.00000000 
-1 .00000000 
0.00000000 
20.00000000 
20.00000000 
1.00000000 
-1 .00000000 
-1 .00000000 
1 .00000000 
20.00000000 
1.00000000 
-1 .00000000 
-1 .00000000 
-1 .00000000 
2.00000000 
20.00000000 
1.00000000 
-1 .00000000 
-1 .00000000 
-1 .00000000 


10  5  5 

o.ocoooooo 

2.60400380 

.90136719 

0.00000000 

0.00000000 

0.00000000 

1.52935000 

.94205000 

.35475000 


4  4 

.07843137 


-.58730000 


-.58730000 


-.58730000 


1  0.00000000 
3.00000000  50 


0.00000000 


0.00000000 


0.00000000 
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3.00000000 
1 .00000000 
-1 .00000000 
-1 .00000000 
-1 .00000000 
-1.00000000 
4.00000000 
-1 .00000000 
-1.00000000 
-1.00000000 
-1 .ooocoooo 
-1 .00000000 
1.52935000 
-.58730000 
C. 00000000 
0.00000000 
0.00000000 
0.00000000 
0.00000000 


.23255000 


-.81985000 


-.58730000 


-.58730000 


0.00000000 


0.00000000 
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SAMPLE  OUTPUT 


( 


* 


o  o  o  o  o  o 
o  o  o  o  o  o 
o  o  o  o  o  o 
o  o  o  o  o  o 

*  o  o  o  o  o  o 

<<•••• 

*  o  —  —  -o  -o  o 

o 

o 


VJ  VJ  V/  V  J  VJ  W 

oooooooo 

OOOOQOOO 

oooooooo 

oooooooo 

•  «  1 

O  —  —  CMfM'O'OO 

-  O 

O 


oooooooooooooooooo 

oooooooooooooooooo 

oooooooooooooooooo 

oooooooooooooooooo 


o - —  oo~ 


—  -O  «  ~  —  -000 

o 

o 


II  III  I  I  I  I  I  I  1  I 

N.N.(v-t>-N.(v.»-N-!^®XoO®aOCD30®>-®®OOOOOOOOOOOOOOOOOOt-0 
(M(MCMP>jrslf'JuJP>ICNC<jfMCMfMf'J<N<S»f'Jyufs|<MPO,Oryi(«nNirorO'<1'<lwirorofr)rriPOi>0(nfOL'JfO 
l/l  </>  VO 

LL  LL  LL 

ac&acacazacoacecxacacacacacacacoac&acacacacacacacacixacacacccacacacccxoac 

i/\  «  r- co  «i  ki  ni«M0>o4i4  i/iOf-cDin-Of'-ooiri'Or'-aovi-cf-coinin 

0-0  0-  >—  O  O-  O-  —  —  •—  CMtMtMCMrO«0«1N0>O 

OOOOOOO  OOOUOOOOO  OOOOOOOUOOOOOOOOOO 


/.D-a 


( 


o 

o  o 

o 

O 

o 

o 

o 

o  o  o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o  o 

o 

O 

o 

o 

o 

o  o  o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

a 

o 

o 

o 

o  o 

o 

o 

o 

o 

o 

o  o  o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o  o 

o 

o 

o 

o 

O 

o  o  o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

a 

o 

o 

o  o 

o 

o 

o 

o 

o 

o  o  o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

• 

• 

•  • 

• 

• 

« 

* 

• 

•  •  • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

« 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

# 

• 

• 
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FORTRAN  LISTING  -  SMOG  3,  PASS  2 


c 

15) 


SMOG3  PASS  2  PRIMAL  SUM  OF  DEV  X  OR  X  AND  Z  DBL  SPL  Ml  PRFL 
DIMENSION  RX(11 ,15),RZ(15,11),FZ(11),FX(15),X(15),Z(11),Y(11 , 


DIMENSION  FPZ(11),FSZ(11),C(13) 

COMMON  JPLS , JUL , JUL 2 , NJUL , M, N , MS , NS , NOROW, NCOMP , RHS , COST 
COMMON  EXP2,CURV  , POINT, EXP1,IWW,ICOL 
COMMON  X,FZ,RZ,Z,FX,FPZ,FSZ,RX,6 
I F( SENSE  SWITCfi  3)750,412 

750  !|9.JP4S,JUL,JUL2,NJUL,M,N,MS,NS,NOROW,NCOMP,RHS,CURV 
READ  181,EXP2, COST, POINT, EXP1,IWW,IC0L 
00  32  1-1, N 
READ  181.X(I),FZ(I) 

DO  32  J-2,M 

32  READ  181,RZ(I,J) 

DO  33  1-1, M 

READ  181.2(1), FX(I),FPZ(I),FSZ(I) 

DO  33  J-1  ,N 

33  READ  181 ,RX( I , J) 

ML-M+2 

DO  100  J-1 .ML 
100  READ  181, C(J) 

C  ,  *  *  X  SEC  DIFFERENCES  *  * 

1,12  fseftui  -  -  - 

DO  88  IZ-1 ,M 
1-1 

JS-JPLS 

JF-JS+JX 

DO  487  J-JS, JF 

IF(RX( IZ, I )-20.)501 ,487,501 

501  PUNCH  162, J 
PUNCH  177,J,RHS 
I  ROW-1 
ISTRT-1 

GO  TO  (503, 504), NCOMP 
503  FUPZ— 1 . 

FUSZ-O. 

ncomp)ALL  PR0FL2  ^<‘)."<'z).«pi,fupz,fusz,exp2,po,nt,t,tp.ts 

*  504  PWR-RX(IZ,I) 

CALL  SEDIF1  (T,TP.TS,X( I ),J, IROW,PWR) 

IROW-IROW+2 

CALL  SEDIF2  (X( l),T,TP,TS,PWR,J, IROW) 

„  IP( 1-3)75,75,502 

502  IROW-IROW+2 
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IF(K-  >73,75,75 
73  DIF-X( l)-X(K) 

CALL  SEDIF2  (DIF,T,TP,TS,PWR,J, IROW) 
503  IROW-IROW+2 

75  IF(IZ-1)87,87.76 

76  00  79  NPWR-1,3 
ISTRT-ISTRT+NJUL 
IROW-ISTRT 

PWR-RX(IZ,I)*Z(IZ)**NPWR 

CALL  SEDIF1  (T,TP,TS,X( l),J, IROW.PWR) 

| ROW* I ROW+2 

CALL  SE0IF2  (X( I ) ,T,TP,TS,PWR,J, IROW) 
IROW-IROW+2 
DO  78  K-3,1.2 
l F(K— I ) 77*79*79 

77  DIF-X(I)-X(K) 

CALL  SEDIF2  (DIF, T,TP,TS,PWR,J, IROW) 

78  IROW-IROW+2 

79  CONTINUE 
IF(IZ-JUL)87,87,5t1 

511  DO  86  L-JUL, IZ, JUL2 
IF(L-IZ)512, 87,87 

512  ISTRT-ISTRT+NJUL 

I ROW- I STRT  %  , 

PWR-RX(IZ,I)*(Z(IZ)-Z(L))**3 
CALL  SEDIF1  (T,TP,TS,X( l),J, IROW.PWR) 
IROW-IROW+2  „  . 

CALL  SEDIF2  (X( I ) ,T,TP,TS,PWR,J, IROW) 
IROW-IROW+2 
DO  85  K-3,1.2 
I F(K- 1)84, 86. 86 
84  DIF-X(I)-X(K) 

CALL  SE0IF2  (DIF,T,TP,TS,PWR,J, IROW) 


85  IROW-IROW+2 

86  CONTINUE 

87  ICOL-ICOL+1 

PUNCH  171 , ICOL,CURV 
PUNCH  152, ICOL,J, ICOL, J 
487  1-1+2 
NCOMP-1 

88  JPLS-JF+1 

Z  SEC  DIFFERENCES  *  * 

DO  110  1-1 »N 


GO  TO  (550.551 )*NCOHP 

550  IF(RX(IZ,I)-20.)556,411,556 
556  PWR-RZ(I.IZ) 

FZD-Z(IZ) 

FOZ-FX(IZ) 

FUPZ-FPZ(IZ) 

FUSZ-FSZ(IZ) 

562  CALL  PR0FL2(X(I),F0Z,EXP1,FUPZ,FUSZ,EXP2, POINT, T, TP, TS, 
NCOHP) 

GO  TO  563 

551  PWR-RZ(I.IZ) 

FZD-Z(IZ) 

563  PUNCH  162, J 
PUNCH  177,J,CURV 

GO  TO  (564,565), NCOMP 
565  ISTRT-1+NJUL 
GO  TO  566 

564  ISTRT-1 

CALL  SE0IF4  (TS,PWR, J, ISTRT) 

567  I F( 1-1 >569,568,569 

569  IROW-ISTRT+2 

DO  570  IXPWR-1,3 

PROD-PWR*X(l)**IXPWR 

CALL  SED I F4( TS , PROD , J , I  ROW) 

570  IROW-IROW+2 

IF( 1-3)568,568,571 

571  DO  572  K-3,1,2 

I F(K-I ) 604, 568, 604 

604  PR0D-PWR*(X(I)-X(K))**3 
CALL  SEDIF4(TS,PR0D,J, I  ROW) 

572  IROW-IROW+2 

568  I STRT- I STRT+NJUL 

CALL  SEDIF5(TP,TS,FZD,PWR,J, ISTRT) 

573  IF( 1-1)575,566,575 

575  IROW-ISTRT+2 

DO  576  IXPWR-1.3 

PROD-PWR*X( I )** I XPWR 

CALL  SEDIF5  (TP,TS,FZD,PROD,J, IROW) 

576  IROW-IROW+2 

I F( 1-3)566,566,577 

577  DO  578  K-3,1,2 

I F(K-I >605,566,605 

605  PROD-PWR*(x( l)-X(K))**3 

CALL  SEDIF5  ( TP, TS,FZD, PROD, J, IROW) 

578  IROW-IROW+2 
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( 


( 


IK 


582 


580 

585 

586 


587 


566  ISTRT-ISTRT+NJUL 

CALL  SEDIF6  (T,TP,TS,FZD,J,PWR, ISTRT) 
IF( 1-1)581,580,581 
IROW-ISTRT+2 
DO  582  IXPWR-1,3 
PROD«PWR*X( I  )**  I XPWR 
CALL  SEDIF6  (T,TP, TS,FZD, J.PROD, I  ROW) 

I ROW" I ROW+2 
I F( 1-3)580,580,583 

583  DO  584  K-3,1.2 

I F(K- I >606,580,606 
606  PR0D«PWR*(X(I)-X(K))**3 

CALL  SEDIF6  (T,TP,TS,FZD,J,PROD, IROW) 

584  I ROW- I ROW+2 
*  ISTRT-ISTRT+NJUL 

CALL  SEDIF2  (FZD.T,TP,TS,PWR, J, ISTRT) 
IF< 1-1)586,590,586 
IROVMSTRT+2 
DO  587  IXPWR-1,3 
PROD-PWR*X( I )*» I XPWR 
CALL  SEDIF2  (FZ0,T,TP,TS,PR0D, J, IROW) 
I ROW- I ROW+2 
„  IF( 1-3)590,590.588 
588  DO  589  K-3,1.2 
I F(K- I >597,590,597 

597  PR0D-PWR*(X(I)-X(K))**3 
CALL  SEDIF2  (FZ0,T,TP,TS,PR0D, J, IROW) 
I ROW- I ROW+2 
I  F( I Z—JUL)591 ,111 ,591 
ISTRT-ISTRT+NJUL 
DO  592  K-JUL. IZ.JUL2 
IF(K-IZ)593. 111,593 

593  DIF-FZD-Z(Kj 

CALL  SEDIF2  (DIF,T,TP,TS,PWR,J, ISTRT) 
IF( 1-1)607,592,607 
IROW-ISTRT+2 
DO  608  IXPWR-1,3 
PROD-PWR*X(l)**IXPWR 
CALL  SEDIF2  (DIF, T, TP, TS, PROD, J, IROW) 
I ROW- I ROW+2 
I F( 1-3)592,592.598 

598  DO  599  KJ-3,1.2  „ 

I F(KJ-I >600,592,600 
600  PR0D-PWR*(X(I)-X(KJ))**3 

CALL  SEDIF2  (DIF, T, TP, TS, PROD, J, IROW) 


589 

590 

591 


1 

607 


608 
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599  I  ROW- 1 ROW+2 
592  I STRT- I STRT+NJUL 
111  ICOL-ICOL+1 

PUNCH  171  JCOL.CURV 
PUNCH  152  JCOL.J.  ICOL.J 


411  J-J+1 

110  NCOMP-1 
PUNCH  174 
PUNCH  176 
PRINT  178 
PAUSE  751 
GO  TO  750 

152  F0RMAT(7X,1HCJ4,2H  R, I4/7X, 1HC, I4.2H  R, 

171  F0RMAT(6X,2H  CJ4.6H  OFSET.F12.5) 

173  FORMAT (1H*. 5 X.47H 

174  F0RMAT(7HFIRST  B) 

175  F0RMATM2X.  6H  R  1.F12.5) 

176  F0RMAT(3HE0F) 

177  F0RMAT(12X,2H  R.I4.1H  .F11.5) 

178  FORMAT (28HENTER  NEW  OATA  -  PRESS  START) 

ISO  FORMAT  (10I5.2F15.9) 

162  F0RMAT(12X,2H  R,l4) 

181  FORMAT  (4F15.9J4.I4) 

END 


I4.8H 


1.0) 
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B-290 


( 


SUBROUTINE  StO  I  FI  (TD,TPD,TSD,XD,JO,  IROWD,PWRD) 
CALL  SEDIF4  (TSD,PWRD, JO, IROWO) 

IROWO-IROWD+2 

CALL  SEDIF5  (TPD.TSD.XD.PWRD.JD, IROWO) 
IROWO-IROWO+2 

CALL  SEDIF6  (TO, TPD,TSD,XD,JD,PWRD, IROWO) 

RETURN 

END 


SUBROUTINE  SEDIF2  (D I FD.TD.TPD.TSD.PWRD, JD, IROWO) 

WRM- (3.*DI FD*TD+3 .*TPD*D I FD**2+ ( TSD*D I FD**3 ) / 2 . ) *PWRD 
CALL  SIGN  (JO, IROWO, WRM) 

RETURN 

ENO 


SUBROUTINE  SEDIF3  ( I XD.FXD.PRD, JSD, IROWO) 
EXD-IXD 

WRH-(EXD*FXD**( I XD-1 )-FXD**l XO)*PRD 
CALL  SIGN  (JSD, IROWO, WRM) 

RETURN 

END 


SUBROUTINE  SEDIF4  (TSD.PWRD, JO, IROWO) 
WRMTSD*PWRD)/2. 

CALL  SIGN  (JO, IROWDjWRM) 

RETURN 

END 


SUBROUTINE  SED I F5  (TPD.TSD,XD,PWRD, JO, IROWO) 

WRM»((TSD*XD)/2.+TPD)*PWRD 

CALL  SIGN  (JO, IROWO, WRM) 

RETURN 

ENO 

( 
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SUBROUTINE  SEDIF6  (TO, TPD,TSD, XD,JD,PWRD, IROVrt)) 
WRH-TD 

WRM-WAM+(TSD*XD**2)/2. 

WRH»WRH+2 ,*TPD*XD 

WRM-WRM*PWRD 

CALL  SIGN  (JO, IROWD,WRM) 

RETURN 

END 


SUBROUTINE  SIGN  (IW,JW,WAM) 

JWPL-JW+1 
I F(WAM) 370,373,371 

370  WAH— WAM 

PUNCH  169, JW, IW.WAM, JWPL, IW, WAM 
RETURN 

371  PUNCH  167» JW, IW, WAM, JWPL, IW.WAM 
373  RETURN 

167  F0RMAT(6X,2H  C,I4,2H  R, I4,1H-,F11.5/6X,2H  C.I4.2H  R.I4.F12.5 
169  F0RHAT(6X,2H  C,I4,2H  R, I4,F12,5/6X,2H  C,I4,2H  R, 14. 1H-,F1 1 .5 
END 


1.0. 0-2 
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TOOD  DLG  5X5  FNO  PROFILE  X-Z  DBL  SPL-I2/3/62 
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Section  IV 


SMOG  4 


The  SMOG  4  program  produces  a  matrix  for  fairing  a  two-dimensional 
curve,  either  double  or  single  splined,  with  the  primal  formulation 
and  sum  of  the  deviations  fit.  The  program  is  written  in  FORTRAN  II 
for  the  IBM- 1620  computer. 


Fortran  Symbol  Definitions 


N 

- 

The  number  of  data  offsets  along  the  curve 

CO  STB 

- 

The  cost  associated  with  slacks 

COSTN 

- 

The  cost  on  the  non-slack  variables 

X(I) 

- 

The  x  coordinate  of  the  offset 

Y(I) 

- 

The  y  coordinate  of  the  offset  if  read  in  feet  and 
decimals. 

AA 

The  number  of  feet  in  the  offset  if  read  in  feet, 
inches,  and  eighths 

BB 

- 

The  number  of  inches  in  the  offset  if  read  in  feet, 
inches,  and  eighths 

CC 

- 

The  number  of  eighths  in  the  offset  if  read  in  feet, 
inches,  and  eighths 

SX 

* 

Scaling  factor  for  normalising  the  x  coordinates. 
It  equals  x^  -  x^  if  set  by  the  program. 

TX 

* 

Translation  factor  for  the  x  coordinates. 

It  equals  Xq  if  set  by  the  program. 
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Input  Card  Formats 

The  card  formats  are  arranged  as  follows: 


Contents  Card  Mo. 


Column 

1-5 

1 

Contains 

any  alphameric  heading 

Format 

15 

F10.4 

F10.4 

2 

Variable 

N 

CO  STB 

00STN 

Alternate 

Cards : 

Format 

F10.4 

F5.1 

F5.1 

F5.1 

Variable 

X(I) 

4  AA 

4  BB 

4  CC 

Next 

or 

N  Cards 

Format 

F10.4 

F10.4 

Variable 

X(I) 

Y(I) 

Format 

F10.4 

F10.4 

Variable 

sx 

TX 

Last  Card 

Output 

The  output  Is  In  a  form  suitable  for  the  1620  Users  Group,  Linear 
Programming  Routine  No.  10.0.006. 


Sense  Switches 


Switch  1  ON 
Switch  1  OFF 
Switch  2  ON 

Switch  2  OFF 
Switch  3  ON 
Switch  3  OFF 


The  curve  will  be  double  splined 

The  curve  will  be  single  splined 

The  scaling  factors  will  be  read  from  the  card 
reader 

The  program  will  set  the  scaling  factors 
The  second  differences  will  be  printed 
The  second  differences  will  not  be  printed 


'  % 


u 
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< 


SMOG  4 

TYPEWRITER  LISTING  OF  SECOND  DIFFERENCES  OF  SAMPLE  DATA 


2ND  DIFF  -  WUNES 
STA  ACT  VALUE  SIGN 

2  .000169  -1.00 

3  .000238  -1 .00 

4  .000068  -1 .00 

5  .000068  -1.00 


( 


( 
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*  *  SAMPLE  INPUT  *  * 


01037 

•  5  999.9999  .001 

87.500000  30.5729 
105.000000  30.307300 

122.500000  30.09375 
140.000000  29.953120 

157.500000  29.833330 

OUTPUT  MATRIX  FROM  SAMPLE  DATA 


D1037 

006  008 

1 

' 

001  001 

1 .0000 

001  002 

— 

1 .0000 

001  003 

1 .0000 

001  004 

— 

1 .0000 

001  005 

1.0000 

001  006 

— 

1.0000 

001  009 

1 .0000 

002  001 

2.0000 

002  002 

— 

2.0000 

002  003 

4.0000 

002  004 

— 

4.0000 

002  005 

8.0000 

002  006 

- 

8.0000 

002  010 

1 .0000 

003  001 

3.0000 

003  002 

- 

3.0000 

003  003 

9.0000 

003  004 

— 

9.0000 

003  005 

27.0000 

003  006 

— 

27.0000 

003  007 

1.0000 

003  008 

— 

1 .0000 

003  011 

1 .0000 

004  001 

4.0000 

004  002 

— 

4.0000 

004  003 

16.0000 

004  004 

mm 

16.0000 

004  005 

64.0000 

004  006 

- 

64.0000 

004  007 

8.0000 

1. 0.0-2 
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004  008  -  8.0000 
004  012  1.0000 
005  003  -  1.0000 
005  004  1 .0000 
005  005  “  6.0000 
005  006  6.0000 
005  013  1.0000 
006  003  -  1.0000 
006  004  1 .0000 
006  005  -  12.0000 
006  006  12.0000 
006  007  -  6.0000 
006  008  6.0000 
006  014  1.0000 
00J  001  101100000J 
OOJ  002  102100000J 
OOJ  003  103100000J 
OOJ  004  104100000J 
OOJ  005  105100000J 
OOJ  006  106100000J 
OOJ  007  107100000J 
OOJ  008  108100000J 
OOJ  009  109199999R 
OOJ  010  110199999R 
OOJ  011  1 1 1 199999R 
OOJ  012  112199999R 
OOJ  013  1 1 3 199999R 
OOJ  014  114199999R 
001  OOJ  101299999R 
002  OOJ  102299999R 
003  OOJ  103299999R 
004  OOJ  104299999R 
005  OOJ  105200000- 
006  OOJ  106200000- 
001  000  -  .2656 
002  000  -  .4791 
003  000  -  .6197 
004  000  -  .7395 
000 


( 
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FLOW  DIAGRAM  -  SMOG  4 
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COSTS 
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ID-2 

TDCOL'O 


CALL  SMOG 
lN,ioa>L,Req)n 


VM.*VAL-S  +  EXP 


W' 

v  y 


o 


r  > 

(  MATRIX 

v  y 
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FORTRAN  LISTING  -  SMOG  4 


★0810 

C 

1 


201 

202 

2 

3 

4 

5 

6 


8 


10 

11 

12 

13 


14 

15 
26 


126 


SMOG  FOR  1620  LP  4/2/63 
DIMENSION  X(30),Y(30),RX(30) 

READ  600 
PUNCH  600 

READ  150, N ,C0STB , COSTN 

JWD-1 

NL-N+1 

DO  2  1-1, N 

| F ( SENSE  SWITCH  4)201,202 
READ  160.X(I),AA.BB.CC  ' 

Y(|)-AA+(BB+(CC/8.))/12. 

GO  TO  2  ,  x 

READ  152,X( I ) ,Y( I ) 

CONTINUE  x 

I F (SENSE  SWITCH  2)3,4 
READ  152,SX,TX 
GO  TO  5  ,  x 

SX-X(2)-X(1) 

TX-X(1)  x,  , 

I F(SENSE  SWITCH  3)6,7 

PRINT  157  .  * 

*  *  SECOND  DIFFERENCES  *  * 

NS-N-1 

DO  13  1*2, NS 

A-2./( (X( 1+1 )-X( 1-1) ) ) 

A^izA/ixT ill  )-X(  I )  )-(Y(  I  )-Y«  1-1  »/€X(  I  )-X<  1-1 ))  ) 
I F( A) 10.9,8 


i r iv,?, 

RX(  I  )—  1. 

GO  TO  11 
RX(l)-0. 

GO  TO  11 

?f|sENSE  SWITCH  3)12,13 
PRINT  158, I ,A,RX( I ) 
CONTINUE 
RX(N)-RX(NS) 

IF (SENSE— SWITCH  3)14,15 
PRINT  158, I ,A,RX(N) 

*  *  NORMALIZING 

DO -26  1-1 .N 

X(  I  )-(X( I )-TX)/SX 

*  V  MATRIX 

IF(SENSE  SWITCH  1)126,127 

NDS-3 

NSD-2 

NC0L-3+N 
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127 

128 


31 

32 

133 

134 

33 


34 

137 


37 


72 

73 

74 

75 

87 


NR0W-N-1+(N-1)/2 
GO  TO  128 
NDS-2 
NSD-1 

NC0L-(N+1)*2 

NROW-2*(N-1) 

MCOL-NCOL 
I  *2 

CALL  SIZE  (NROW,MCOL) 

DO  37  I  ROW-1, NS 
JCOL-1 

DO  31  NPR-1.3 
TRM-X( I )**NPR 

CALL  MATRIX( IROW, JCOL,TRM, JWD) 

JCOL-JCOL+2 

IF(l-NOS)137,137,32 

DO  34  IT-NOS.I 

IF(I-IT)34, 34.133  , 

IF (SENSE  SWITCH  1)134,33 

IF(lT-2*(IT/2))33.34,33 

TRM-(X(I)-X(IT))**3 

CALL  MATRI X( IROW, JCOL,TRM, JWD) 

JCOL-JCOL+2 

CONTINUE 

NCOL-NCOL+1 

JWD— 1 

TRM-1  • 

CALL  MATRI X( IROW,NCOL,TRM,JWO) 

JWD-1 
1-1  +  1 

*  *  X  SECOND  DIFFERENCES 

DO  88  l-NDS,N,NSO 
JCOL-3 

CALL  HATRI X( IROW, JCOL,RX( I ) , JWD) 

JCOL-JCOL+2 

TRM-3.*X( I )*RX( I ) 

CALL  MATRI X( IROW, JCOL,TRH, JWD) 

JCOL-JCOL+2 

I F  ( I  — NDS )  75 » 75  •  72 

DO  74  K-NDS.I.NSD 

IF(K-1 >73.75,75  %%  #  % 

TRM-3.*(X\ I )-X(K) )*RX( I ) 

CALL  MATRI X( IROW, JCOL,TRM, JWD) 

JCOL-JCOL+2 

CONTINUE 

*  *  SLACKS  *  * 

NCOL-NCOL+1 
JWD— 1 


★ 


★ 
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O 


< 


< 


TRM-1. 

CALL  MATRI X( IROW,NCOL,TRM,JWD) 
IROW-IROW+1 
88  JWO-1 

C  *  *  COSTS 

00  91  J-1,MC0L 
91  CALL  C0SN0N(J,C0STN) 

I  ROW- 1  ROW- 1 
MCOL-MCOL+ 1  ROW 


DO  300  l«J,MCOL 

300  CALL  COSNON(l , COSTS) 

DO  92  J-1 , 1  ROW 
IF(J-N)92,95,92 

95  COSTB-O, 

92  CALL  COSBAS( J ,COSTB) 

★  *  R  H  S 

00  94  1-1, NS 
Y(I+1)»Y(I+1)-Y(1) 

94  CALL  RHS( I ,Y( 1+1)) 

PUNCH  161 
PAUSE 
GO  TO  1 

150  FORMAT ( I5,2F10.6) 

152  F0RMAT(2F10,4) 

157  FORMAT ( 10X, 1 7H2ND  DIFF  -  WLINES// 

C  8X,22HSTA  ACT  VALUE  SIGN/) 

158  FORMAT (5X, I5,F12.6.F8,2) 

160  F0RMAT(F10.5.3F5.1) 

161  FORMAT (3HOOO) 

600  FORMAT (5H  ) 

601  F0RMAT(19X,I4,I4,2H  1) 

END 


*0810 

SUBROUTINE  C0SBAS(M,C0STB) 

MM-MflOO 

ID-1 

IDCOL—  1 

MC0ST-C0STB*1000. 

I V  AL -MM* 10000000+ MCOST+ 2000000 
CALL  SMOG  (M, IDCOL, IVAL, ID) 
RETURN 
END 


*0810 

SUBROUTINE  COSNON(M,COST) 

MM-MflOO 

ID-1 

IDROW—  1 

NCOST-COST  *1000. 


* 


* 


* 
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*0810 


i  VAL-MM*10000000+NCOST+1000000 
CALL  SHOG  (IDROW.M,  IVAL.IO) 
RETURN 
END 


SUBROUTINE  RHS(N,REQ) 

ID-2 
IDCOL-O 

CALL  SHOG  (N,IDCOL,REQ,ID) 
RETURN 
END 

*0810 


SUBROUTINE  MATRI X(N,M,VAL, JND) 
IF(JND)509,508,508 

508  ID-3 

CALL  SMOG  (N,M,VAL,ID) 

RETURN 

509  ID-2 

CALL  SHOG  (N,M,VAL, ID) 

RETURN 

END 
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Appendix  C 


OPERATING  PROCEDURES  FOR  LINEAR  PROGRAMING 


There  are  several  considerations  in  the  fairing  of  a  surface  using  linear 
programming  which  may  have  a  bearing  on  the  final  result  attained.  The 
following  is  a  short  discussion  of  these. 

SENSE  OF  COORDINATE  SYSTEM 

The  coordinate  system  employed  in  the  fairing  procedure  is  very  general. 
The  origin  of  each  surface  faired  is  at  the  intersection  of  the  first 
waterline  and  the  first  station  in  the  surface.  Thus  it  may  be  at  any 
one  of  the  four  corners  of  the  surface,  depending  on  which  station  and 

which  waterline  is  chosen  as  "first." 

r 

n 

For  surfaces  without  profiles  it  is  generally  more  convenient  to  choose 
the  origin  to  be  at  the  corner  with  the  least  value  for  the  offset. 

Since  this  offset  becomes  the  constant  for  the  equation  and  is  subtracted 
from  the  other  offsets,  this  allows  the  remaining  portion  of  the  other 
offsets  to  be  positive.  On  most  ships  the  origin  then  becomes  the  lowest 
waterline  on  the  surface  and  the  station  furthest  forward  if  the  surface 
is  forward  on  the  ship,  or  aft  if  the  surface  is  aft  on  the  ship. 

For  surfaces  with  profiles,  the  origin  of  the  surface  must  be  the  same 
as  the  origin  of  the  profile  equation.  Therefore,  for  the  bow  the  origin 
of  both  the  surface  and  profile  is  usually  made  the  intersection  of  the 
forward  perpendicular  and  the  baseline;  aft,  the  origin  is  the  inter¬ 
section  of  the  aft  perpendicular  and  the  baseline. 

The  positive  sense  of  the  coordinates  is  always  along  the  surface.  If 
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the  origin  is  at  the  lower-right  corner,  for  instance,  the  positive  x 
direction  is  to  the  left,  and  the  positive  z  direction  is  up. 

BASIC  FEASIBLE  SOLUTION 

Mcst  linear  programming  codes,  such  as  LP-90  and  the  code  for  the 
IBM-1620,  have  two  inodes  of  operation.  That  is,  when  they  first  start 
to  solve  a  problem  the  Phase  I  Algorithm  is  used  to  find  a  first  feasible 
solution,  without  regard  to  optimality;  and  after  finding  this.  Phase  II 
is  used  to  find  an  optimum  solution. 

In  larger  problems  the  search  for  the  first  basic  feasible  soluticn  is 
time  consuming  and,  due  to  the  increase  in  the  number  of  calculations, 
encourages  the  buildup  of  round-off  error.  It  has  been  found  advan¬ 
tageous  then,  in  these  larger  problems,  to  provide  the  program  with  a 
basic  feasible  soluticn  so  the  optimization  can  begin  at  once. 

This  is  done  by  producing  in  the  matrix  an  identity  matrix  of  rank  equal 
to  the  number  of  rows  in  the  matrix.  The  identity  matrix  is  identified 
to  the  program  for  use  as  the  first  basic  feasible  solution. 

Each  of  the  SMOG  programs  has  a  routine  for  generating  this  basis  and 
identifying  it  to  the  linear  program.  In  SMOG  3  and  SMOG  4,  this  routine 
is  self-contained.  For  SMOG  1  and  SMOG  2,  the  routine  of  Appendix  K  is 
used. 


INVERSION 

During  the  solution  of  a  linear  programing  problem,  difficulties  may 
arise  from  three  sources :  first,  from  the  increase  in  number  of  non-zero 
elements  in  the  matrix,  which  causes  the  iteration  times  to  become  greater. 
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the  origin  is  at  the  lcwer-rlght  corner,  for  instance,  the  positive  x 
direction  is  to  the  left,  and  the  positive  z  direction  is  up. 

BASIC  FEASIBLE  SOLUTION 

Mcst  linear  programming  codes,  such  as  IP-90  and  the  code  for  the 
IBM-1620,  have  two  modes  of  operation.  That  is,  when  they  first  start 
to  solve  a  problem  the  Phase  I  Algorithm  is  used  to  find  a  first  feasible 
solution,  without  regard  to  optimality:  and  after  finding  this,  Phase  II 
is  used  to  find  an  optimum  solution. 

In  larger  problems  the  search  for  the  first  basic  feasible  solution  is 
time  consuming  and,  due  to  the  increase  in  the  number  of  calculations, 
encourages  the  buildup  of  round-off  error.  It  has  been  found  advan¬ 
tageous  then,  in  these  larger  problems,  to  provide  the  program  with  a 
basic  feasible  solution  so  the  optimization  can  begin  at  once. 

This  is  done  by  producing  in  the  matrix  an  identity  matrix  of  rank  equal 
to  the  number  of  rows  in  the  matrix.  The  identity  matrix  is  identified 
to  the  program  for  use  as  the  first  basic  feasible  solution. 

Each  of  the  SMOG  programs  has  a  routine  for  generating  this  basis  and 
identifying  it  to  the  linear  program.  In  SMOG  3  and  SMOG  _4,  this  routine 
is  self-contained.  For  SMOG  ]  and  SMOG  2,  the  routine  cf  Appendix  K  is 
used. 


INVERSION 

During  the  solution  of  a  linear  programming  problem,  difficulties  may 
arise  from  three  sources  :  first,  from  the  increase  in  number  of  non-zero 
elements  in  the  matrix,  which  causes  the  iteration  times  to  become  greater. 
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Second,  because  of  a  buildup  of  round-off  error,  some  of  the  numbers  may 
lose  significance.  Third,  from  the  buildup  of  the  number  of  transforma¬ 
tions  in  the  product  form  of  the  inversion,  the  time  per  iteration 
increases . 

Most  linear  programming  routinas  which  are  meant  to  handle  large  problems, 
such  as  If -90,  have  proviaion  for  inverting  the  basis  matrix.  This 
inversion,  carried  out  at  intervals  during  the  solution  of  the  problem, 
will  be  extremely  helpful  in  reducing  all  of  the  difficulties  mentioned 
above . 

The  best  method  for  choosing  the  inversion  Interval  seams  to  be  that  of 
simply  keeping  a  running  check  on  the  time  per  Iteration.  The  iteration 
time  will  be  found  to  increase  slowly  for  several  Iterations,  then  to 
increase  more  rapidly.  At  this  point  it  is  time  to  Invert.  On  surfaces 
of  eighty  to  one-hundred  offsets  this  usually  occurs  at  about  every 
thirty  iterations.  It  is  not  efficient  to  Invert  in  any  case  at 
intervals  of  less  than  twenty  to  twenty-five  iterations,  because  the 
inversion  time  cancels  any  savings  in  iteration  time. 


